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Abstract 

We outline in detail the properties of generalized parton distributions (GPDs), 
which contain new information on the structure of hadrons and which enter the de- 
scription of hard exclusive reactions. We highlight the physics content of the GPDs 
and discuss the quark GPDs in the large Nc limit and within the context of the chiral 
quark-soliton model. Guided by this physics, we then present a general parametriza- 
tion for these GPDs. Subsequently we discuss how these GPDs enter in a wide variety 
of hard electroproduction processes and how they can be accessed from them. We 
consider in detail deeply virtual Compton scattering and the hard electroproduction 
of mesons. We identify a list of key observables which are sensitive to the various 
hadron structure aspects contained in the GPDs and which can be addressed by 
present and future experiments. 
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1 INTRODUCTION 

The famous Lagrangian of Quantum Chromodynamics (QCD): 



contains, in principle, all phenomena of hadronic and nuclear physics ranging from the 
physics of pions to the properties of heavy nuclei. The main difficulty and challenge in the 
derivation of strong interaction phenomena from the QCD Lagrangian (|^) is that the theory 
is formulated in terms of quark and gluon degrees of freedom whereas phenomenologically 
one deals with hadronic degrees of freedom. Understanding of how colorless hadrons 
are built out of colored degrees of freedom (quarks and gluons) would allow us to make 
predictions for strong interaction phenomena directly from the Lagrangian (|l|) . The physics 
of hadronization of quarks and gluons is governed by such phenomena as confinement and 
spontaneous chiral symmetry breaking which are in turn related to the non-trivial structure 
of the QCD vacuum. It implies that the studies of hadronization processes provide us with 
valuable information on the fundamental questions of the vacuum structure of non-abelian 
gauge theories. 

The quark and gluon structure of hadrons can be best revealed with the help of 
weakly interacting probes, such as (provided by Nature) photons and W , Z bosons. One 
needs probes which are weakly coupled to quarks in order to "select" a well defined QCD 
operator expressed in terms of quark and gluon degrees of freedom of the Lagrangian (Q). 
By measuring the reaction of a hadron to such a probe, one measures the matrix element 
of the well-defined quark-gluon operator over the hadron state revealing the quark-gluon 
structure of the hadron. 



Historically the famous experiment of Otto Stern et al. | Fri33 measuring the anoma- 
lous magnetic moment of the proton revealed for the first time the non-trivial (now we 
say quark-gluon) structure of the proton. In modern language we can say that in this ex- 
periment with the help of the low energy photon probe (weakly interacting to quarks) one 
selects the QCD operator Cqq^^q, then one couples it to the proton, (p '| eqqj'^q\p), 
and extracts the structural information encoded in the values of hadron form factors. The 
experimental observation of a large Pauli form factor i^2(0) of the proton |Fri33| has 



created the physics of hadrons as strongly interacting many body systems. 

Are we limited to explore the structure of hadrons by QCD operators created by 
photons and W, Z bosons? Can we find other weak coupling mechanism to select more so- 
phisticated QCD operators to explore the structure of hadrons? Actually weak interactions 
as such are an inherent property of QCD due to the phenomenon of asymptotic freedom 
meaning that at short distances the interactions between quarks and gluons become weak. 
This implies that if one manages to create a small size configuration of quarks and gluons 
it can be used as a new probe of hadronic structure. The possibility to create small size 
configurations of quarks and gluons is provided by hard reactions, such as deep inelastic 
scattering (DIS), semi-inclusive DIS, Drell-Yan processes, hard exclusive reactions, etc. 

^ The irony is that this form factor carries the name of W. Pauli who tried to talk out O. Stern of his 
experiment, because for all theoreticians at that time th ere were no doubts that the proton is a point like 



particle. See the nice historical overview in Ref. | DrcOO 
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The common important feature of hard reactions is the possibiHty to separate clearly 
the perturbative and nonperturbative stages of the interactions, this is the so-called factor- 
ization property. Qualitatively speaking, the presence of a hard probe allows us to create 
small size quark, antiquark and gluon configurations whose interactions are described by 
means of perturbation theory due to the asymptotic freedom of QCD. The non-perturbative 
stage of such a reaction describes how a given hadron reacts to this configuration, or how 
this probe is transformed into hadrons. 

In the present work we deal with hard exclusive reactions of the type: 



7*((Z) + T{p) ^ 7(q ') + T'{p ') , 72(g) + T{p) ^ M{q ') + T'{p ') , 



(2) 



in which a photon 7* with high energy and large virtuality —q^ — > Q scatters off 
the hadronic target T and produces a meson M (or a low mass mesonic cluster) or a 
real photon 7, and a low-mass hadronic state T' (where the mass is small compared to 
Q^). The all order factorization theorems for such kind of reactions have been proven in 
Ref. |Col97| (meson production) Refs. |Ji98a, Col99| , Rad98| (real photon production) . 
These theorems assert that the amplitude of the reactions (^ can be written in the form 
(we show this for the case of hard meson production |Col97|): 



XB ; m) Hi-j [Q^xi/xB , Q^, z, m) (I^j {z, fJ.) 



-|-power-suppressed corrections, 



(3) 



where fi/p is a generahzed parton distribution (GPD), |Bar82[ pit88| , |Mul94| , |Ji97b| , |Rad9^ , 
C0I97] (see details in Sec. xi is the fraction of the target momentum carried by the 
interacting parton, (f>j is the distribution amplitude (DA) of the meson, and Hij is a 
hard-scattering coefficient, computable as a power series in the strong coupling constant 
as{Q). The amplitude depends also on the Bjorken variable xb = Q^/'2,{pq), and on the 
momentum transfer squared t = [p — p which is assumed to be much smaller than 
the hard scale Q^. In Eq. (|^) the GPDs fi/p and the meson DAs 4>-j contain the non- 
perturbative physics and Hij the perturbative one. The proof of cancellation of the soft 
gluon interactions in the processes (||) is intimately related to the fact that the final meson 
arises from a quark-antiquark (gluon) pair generated by the hard interaction. Thus the 
pair starts as a small size configuration and only substantially later grows to a normal 
hadronic size, i.e. into a meson. 

Qualitatively one can say that the reactions (^) allow one to perform a "microsurgery" 
of a nucleon by removing in a controlled way a quark of one flavor and spin and implanting 
instead another quark (in general with a different flavor and spin). It is illustrated in Fig. |l| 
for the case of the deeply virtual Compton scattering (DVCS) and in Fig. || 
for hard exclusive meson production (HMP). The lower blobs in both flgures correspond to 
GPDs. It is important that according to the QCD factorization theorem in both types of 
processes the same universal distribution functions enter. This allows us to relate various 
hard processes to each other. Additionally the detected final state can be used as a filter 
for the spin, flavor, C-parity, etc. of the removed and implanted quarks. 

In the processes as in Eq. (^ the short distance stage of the reaction described by Hij 
in Eq. (^ corresponds to the interaction of a parton with a highly virtual photon. This 
stage is described by perturbative QCD. The corresponding hard scattering coefficients 
have been computed to the NLO order in Refs. [|Ji984 ^el98b|, [Man98bt for DVCS and 
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Figure 1: "Handbag" diagrams for DVCS. 




recently in Ref. [Bel01c| for the hard pion production. Also the perturbative evolution of 
generalized parton distributions is elaborated to the NLO order [ Bel99 , BelOOd BelOOe . 
This shows that the perturbative QCD calculations of the hard scattering coefficients 
and evolution kernels are under theoretical control. In this work we shall not review the 
corresponding perturbative calculations, for a review see Refs. |Ji98b| , iRadOlbj , and for 
our analysis we will always stick to the LO expressions. 

We can say that the hard scattering part ("handbag part" of Fig. |^ and the part 
denoted as Tfj in Fig. ||) of the process "creates" a well-defined QCD operator which is 
"placed" into the target nucleon and the outgoing meson. We concentrate in present work 
mostly on studies of such non-perturbative objects entering the factorization theorem (j^). 
They can be generically described by the following matrix elements: 



-I ab 



(4) 



where the operators are on the light-cone, i.e. = 0. Depending on the initial and the 
final states A and B, the following nomenclature of the non-perturbative matrix elements 
exists: 
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A = B = one particle state - parton distributions in the hadron A. These objects 
enter the description of hard inclusive and semi-inclusive reactions. 

B = vacuum, A = one particle state - distribution amplitudes or light-cone wave 
functions of a hadron A. They enter the description of, say, hadronic form factors at 
large momentum transfer as well as hard exclusive production of mesons. 

Both A and B are one particle states with different momenta, i.e. p{A) ^ p{B) - 
generalized parton distributions. They enter the description of hard exclusive pro- 
duction of mesons and deeply virtual Compton scattering. 

B — vacuum, A= many particle state - generalized distribution amplitudes. They 
enter the description of hard exclusive multi-meson production or transition form 
factors between multi- meson states at large momentum transfer. 

A = one particle state, B 



many particle state, e.g. B = irN, KN, etc. - can be 
called many body generalized parton distributions. They enter the description of 
semi-exclusive production of mesons and semi-exclusive DVCS. 

• A = B = many particle state - can be called interference parton distributions. An 
interesting object which has not yet been considered in the literature. 

Although all these objects have different physical meaning, their unifying feature is that 
they describe the response of hadronic states to quark-gluon configurations with small 
transverse size, which have an extension along the light-cone direction. This can also be 
rephrased by saying that one probes the transition from one hadronic state to another by 
low-energy extended objects- "QCD string operators" (Q). From such a point of view, the 
studies of hard exclusive reactions allow us to extend considerably the number of low-energy 
fundamental probes that can be used to study hadrons. An advantage of hard exclusive 
reactions is that in this case the number of probes with different quantum numbers is very 
large and also that these reactions probe non-diagonal transitions in flavor and spin spaces. 
To great extent, the properties of these new probes have not been studied, especially such 
aspects as chiral perturbation theory for these probe^, the large limit, the influence 
of spontaneous breaking of the chiral symmetry, etc. In the present work we attempt to 
address these questions. 
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Table 1: Combinations of parton distributions entering the amplitude for the production of a 
meson cluster with quantum numbers indicated in the first row (J^*^) and with zero isospin. In 
the second row we give typical meson representatives for these quantum numbers. 
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As an illustration of the wider scope of hard exclusive reactions for studies of hadronic 
structure we give here a qualitative discussion to what type of usual parton distributions 
these reactions are sensitive. We shall see that the hard exclusive reactions allow us (leaving 
aside truly non-forward effects) to probe flavor and C-parity combinations of usual parton 
distributions which are not accessible in inclusive measurements. The amplitude of a hard 
exclusive meson production on the nucleon depends on the nucleon generalized parton 
distributions. In the limiting case of zero "skewedness" (i.e. when the momentum transfer 
A ^ in Figs. |^,^) the latter are reduced to usual parton distributions in the nucleon. 
In this sense it holds that the hard exclusive meson production is, amon g othe r things, 
sensitive to the usual parton distributions in the nucleon. In Tables |Fra01 | we give 
examples of hard exclusive production processes of mesons with various quantum numbers 
on a proton target and list the corresponding parton distributions to which the amplitudes 
are sensitive. In these tables the non empty entries mean that the given combination of 
parton distributions can be probed in the production of a mesonic system with the given 
quantum numbers. The physics content of the GPDs is of course not reduced to the 





Q-+ 


0++ 


1 — 


1+- 


1++ 


2++ 






ao 




bi 


ai 




2{u -u) + {d-d) 




X 








X 


2{u + u) + {d + d)+'^G 






X 








2{Au - Au) + {Ad - Ad) 


X 








X 




2{Au + Au) + {Ad + Ad) 








X 







Table 2: Combinations of parton distributions entering the amplitude for the production am- 
plitude of a meson cluster with the quantum numbers indicated in the first row ( J^*^) and with 
isospin one and zero charge. In the second row we give typical meson representatives for these 
quantum numbers. 
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Table 3: Combinations of parton distributions entering the amplitude for the production am- 
plitude of a meson cluster with the quantum numbers indicated in the first row {J^^ , instead of 
C— parity we indicate the G-parity) and with charge +1. In the second row we give typical meson 
representatives for these quantum numbers. 

forward quark distributions. In fact, the "truly non-forward" effects play a crucial role in 
the physics of hard exclusive processes (^), i.e. HMP as well as DVCS. These effects also 
contain new information about the structure of mesons and baryons which is not accessible 
by standard hard inclusive reactions. 

The paper is organized as follows. In section]^ we review the formalism of generalized 
parton distributions. We concentrate on such properties of quark GPDs which are related 
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to the structure of the nucleon leaving aside a discussion of the perturbative evolution 



Section || is devoted to 
In this chapter we also 
The basic 



which has been extensively reviewed in Refs. |Ji98b, RadOlt]. 
studies of quark GPDs in the limit of a large number of colors 
review and present new calculations of GPDs in the chiral quark-soliton model 
structures of GPDs discussed in the sections |^ and ^ are then used in chapter |^ in order to 
construct a phenomenological parametrization of the GPDs. These parametrizations are 
then used to make predictions for a wide variety of hard electroproduction processes. Since 
our main objective is to investigate how to reveal the basic features of the nucleon structure 
in these reactions we restrict ourselves in the calculation of the hard interaction kernels to 
the leading perturbati ve order for ob s ervable s. Studies of the NLO effects on observables 
can be found in Refs. [ Bel98a , BelOOa Bel01c[ ^ In particular, we perform detailed studies 
for DVCS observables in Sec. 4.3 at twist-2 and twist-3 level. We demonstrate how various 
DVCS observables are sensitive to parameters of G PDs , e.g. to the quark total angular 
momentum contribution to the nucleon spin. In Sec. 4^ we study hard meson production 
discussing in detail for each specific process the information about nucleon structure one 
accesses. In all calculations we restrict ourselves to the region of not too small Xb (valence 
region), where one is mainly sensitive to the quark GPDs. Discussions of various aspects 



Bro94 


, Fra96 


Rcf. 


Fra98b 



Fra99a|. We note that in the hard exclusive processes at small xb the GPDs can be 
related to the usual parton distributions in a model independent way. 

The GPDs have also been discussed within the context of wide angle Compton scatter- 
ing I Rad94 Die99 1 where the GPDs are related to soft overlap contributions to elastic and 
generalized form factors (for a review of such processes see Ref. |Rad01b |). Additionally 
the GPDs enter the description of the diffractive photoproduction of heavy quarkonia, see 
e.g. Refs. [ |Rys97| , |Fra984 |Fra99c|, [Mar99t, as weU as in diffractive 7 2jets ]Gol9^ | and in 
TT 2jets processes |Fra9S, FraOOb, BraOl, pher01 |. We do not review these applications 
in the present work. 

Finally, we summarize our main findings in section |^ and list a number of key exper- 
iments to access the GPDs and to extract their physics content. 



^In these papers it was found that for some of observables the NLO corrections can be rather large. 
This implies that more detail studies of the NLO effects arc needed. 
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2 GENERALIZED PARTON DISTRIBUTIONS (gpds) 



In this section we consider basic properties of the generalized parton distributions (GPDs). 
We shall pay special attention to such properties of GPDs which are related to the non- 
perturbative structure of hadrons and the QCD vacuum. We shall not discus s the per - 
turbative evolution of the GPDs, as this issue is excellently reviewed in Ji98b , RadOlbj. 
Our aim will be to cov er topics which are complementary to the ones discussed in those 
reviews fJiQSU |Rad01b|] . 



2.1 Definitions, link with ordinary parton distributions and nu- 
cleon form factors 

The factorization theorem for hard exclusive reactions allows us to describe the wide 
class of such reactions in terms of universal generalized parton distributions. This nu- 
cleon structure information can be parametrized, at leading twist-2 level, in terms of 
four (quark chirality conserving) generalized structure functions. These functions are the 
GPDs de noted by H, H, E, E (we do not consider chirally odd GPDs, which are discussed 
in Rcfs. [ Hoo98 , Die01b |) which depend upon three variables : x, £, and t. The light-cone 
momentum ^ fraction x is defined (see Figs. ^ and ^ by fc+ = xP^ , where k is the quark 
loop momentum and P is the average nucleon momentum {P — {p + p where p{p ') 
are the initial (final) nucleon four- momenta respectively). The skewedness variable ^ is 
defined by A"*" = — 2^^"*", where A = p ' — p is the overall momentum transfer in the pro- 
cess, and where 2^ — > xb/{^ — Xb/^) in the Bjorken limit. Furthermore, the third variable 
entering the GPDs is given by the Mandelstam invariant t = A^, being the total squared 
momentum transfer to the nucleon. |^ In a frame where the virtual photon momentum 
and the average nucleon momentum P^ are collinear along the z-axis and in opposite 
direction, one can parametrize the non-perturbative object in the lower blobs of Figs. ^ 
and I as (we follow Ref. ]ji97b| ) |: 



27 



v+=y±=o 



H^{x,i,t) N{p)j+N{p) 
+ E%x,tt) N{p')ic7+^ 



A. 
2m M 



Nip) 



+ (757")a/3 [H'^{x,tt) Nip h+-/5N{p) 

A+ 



+ E\x,^,t) 7V(p)7E 



2mN 



N{p) 



(5) 



where ijj is the quark field of flavor q, N the nucleon spinor and niN the nucleon mass. The 
Ihs of Eq. (|^) can be interpreted as a Fourier integral along the light-cone distance y~ of 

^using the definition = l/\/2{a^ it a^) for the Hght-cone components 

^In what follows we shall omit the t-dependence of GPDs if the corresponding quantity is assumed to 
be defined at i = 0. 

^In all non-local expressions we always assume the gauge link Pexp(jg J daj'^A^) ensuring the color 
gauge invariance of expressions. 
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a quark-quark correlation function, representing the process where a quark is taken out of 
the initial nucleon (having momentum p) at the space-time point y/2, and is put back in 
the final nucleon (having momentum p ') at the space-time point —y/2. This process takes 
place at equal light-cone time {y~^ = 0) and at zero transverse separation {y± = 0) between 
the quarks. The resulting one-dimensional Fourier integral along the light-cone distance t/~ 
is with respect to the quark light-cone momentum xP~^ . The rhs of Eq. (|^) parametrizes 
this non-perturbative object in terms of four GPDs, according to whether they correspond 
to a vector operator (7~)q/3 or an axial-vector operator (757~)o!/3 at the quark level. The 
vector operator corresponds at the nucleon side to a vector transition (parametrized by 
the function H'^, for a quark of flavor q) and a tensor transition (parametrized by the 
function E"^). The axial- vector operator corresponds at the nucleon side to an axial- vector 
transition (function Sf) and a pseudoscalar transition (function E'^). 

In Fig. |l|, the variable x runs from -1 to 1. Therefore, the momentum fractions {x + £, 
or X — ^) of the active quarks can either be positive or negative. Since positive ( negative) 



momentum fractions correspond to quarks (antiquarks), it has been noted in [Elad96a 
that in this way, one can identify two regions for the GPDs : when x > £^ both partons 
represent quarks, whereas for x < — ^ both partons represent antiquarks. In these regions, 
the GPDs are the generalizations of the usual parton distributions from DIS. Actually, in 
the forward direction, the GPDs H and H reduce to the quark density distribution q{x) 
and quark helicity distribution Aq{x) respectively, obtained from DIS : 

^'(^,0,0) = I (6) 



^'(-.0-0) = { ^^''li ^ ; (7) 



-q{—x), X < 0. 

A(j(a;), a; > 0. 

Aq{-x), x < 0. 

The functions E and E are not measurable through DIS because the associated tensors in 
Eq. (H) vanish in the forward limit (A 0). Therefore, E and E are new leading twist 
functions, which are accessible through the hard exclusive electroproduction reactions, 
discussed in the following. 

In the region — ^ < x < ^, one parton connected to the lower blob in Fig. |l] represents 
a quark and the other one an antiquark. In this region, the GPDs behave like a meson 
distribution amplitude and contain completely new information about nucleon structure, 
because the region — ^ < a; < ^ is absent in DIS, which corresponds to the limit ^ — > 0. 

Besides coinciding with the quark distributions at vanishing momentum transfer, 
the generalized parton distributions have interesting links with other nucleon structure 
quantities. The first moments of the GPDs are related to the elastic form factors of the 
nucleon through model independent sum rules. By integrating Eq. (H) over x, one obtains 



for any ^ the following relations for a particular quark flavor |Ji97b 



dxH''{x,^,t) = Flit), (8) 

1 

+1 

dxE%x,S,,t) = F^it), (9) 

1 

+1 

dxH''{x,^,t) = g\{t), (10) 

-1 
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+1 



(11) 



where F^(t) represents the elastic Dirac form factor for the quark flavor q in the nucleon. 
When referring to the quark form factors in the following, we understand them in our 
notation to be for the proton, e.g. ^"(t) = F^^^{t). In this notation, the u-quark form 
factor is normalized at t = as i^i'(O) = 2 so as to yield the normalization of 2 for the 
w-quark distribution in the proton, whereas the d-quark form factor is normalized at i = 
as F^{Q) = 1 so as to yield the normalization of 1 for the d-quark distribution in the 
proton. These elastic form factors for one quark flavor in the proton, are then related to 
the physical nucleon form factors (restricting oneself to the u, d and s quark flavors), using 
SU{2) isospin symmetry, as : 



u/p 



f: 



d/p 



2Ff 
2Fi" 



Ff 



F? 



(12) 



where Ff and f " are the proton and neutron electromagnetic form factors respectively, 
with Ff (0) = 1 and F"(0) = 0. In Eq. (|l^) Ff is the strangeness form factor of the 
nucleon. Relations similar to Eq. ( p^ ) hold for the Pauli form factors . For the axial 
vector form factors one uses the isospin decomposition : 



9a 



^9a + ^9a 



9i 



-^9 A + ^9 A 



(13) 



where 3^(3^4) are the isovector (isoscalar) axial form factors of the nucleon respectively. 
Similar relations exist for Ha- The isovector axial form factor qa is known from experiment, 
with (7yi(0) ~ 1.267. The induced pseudoscalar form factor Ha contains an important pion 
pole contribution, through the partial conservation of the axial current (PCAC), as will 
be discussed in more details below. 

An interesting connection of GPDs to light-cone wave functions of the nucleon was 
considered in Refs. |Dic01a, BroOl| . In these papers an exact representation of generalized 
parton distributions for unpolarised and polarized quarks and gluons as an overlap of the 
nucleon light-cone wave functions has been constructed. Using such an overlap representa- 
tion of the GPDs model calculations have been presented in Ref. [ChoOl, TibOl, BurkOlj. 



2.2 Polynomial conditions, D-term 

One of the non-trivial properties of the generalized parton distributions is the polynomiality 
of their Mellin moments which follows from the Lorcntz invariance of nucleon matrix 
elements |Ji97b|. Indeed the {N + l)-th Mellin moment of GPDs corresponds to the 
nucleon matrix element of the twist-2, spin-(iV-f 1) local operator. Lorentz invariance then 
dictates that the Mellin moments of GPDs should be polynomials maximally of the order 
N - 



1, 



the polynomiality property implies that [Ji97b 



(14) 



^ The g eneral method for counting of generalized form factors of twist-2 operators can be found in 
Ref. IjiOll 
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dx x^ Ei{x,0 = 



Note that the corresponding polynomials contain only even powers of the skcwcdness 
parameter ^. This follows from the time reversal invariance, see Ref. Man98a , Ji98b |. 
This fact implies that the highest power of ^ is iV + 1 for odd N (singlet GPDs ) and N 
for even N (nonsinglet GPDs). Furthermore due to the fact that the nucleon has spin 1/2, 
the coefficients in fro nt of the high est power of ^ for the singlet functions H'^ and E"^ are 
related to each other |Ji97t, Ji98b|: 



q{N) _ _,q{N) 
f^N+l — "'N+1 



(15) 



The polynomiality conditions ( p^ strongly restrict the class of functions of two variables 
H'^{x,C) and E^{x,£). For example the conditions ( p^ imply that GPDs should satisfy 
the following integral constrains: 



dx 
1 X 



H'^{x,^ + xz)-H'^{x,^) 





dx 




J-l X 



E'^{x,£,+xz)-E'^{x,£,) 



n=0 



(16) 



Note that the skewedness parameter ^ enters the Ihs of this equation, whereas the rhs of 
the equation is ^-independent. Therefore this ^-independence of the above integrals is a 
criterion of whether the functions 77'^(a;,^), E'^{x,C) satisfy the polynomiality conditions 
( [l4|) . Simultaneously these integrals are generating functions for the highest coefficients 
n^^i- In addition, the condition ( p^ shows that there are nontrivial functional relations 
between the functions H'^{x,£^) and E'^{x,£^). 

An elegant possibility to implement the polynomiality conditions (|l^) for the GPDs 
is to use the double distributions |Mul94, Rad96a, Elad97 . A detailed discussion of the 



double distributions has been given in the review of Ref. |Rad011:]. In this case the gener- 
alized distributions are obtained as a one-dimensional section of the two-variable double 



distributions F^(/?, a), if (/?, a) (Recently in Refs. BelOOc , TerOl the inversion formula 
has been discussed): 



/•I ri-m 

Hdd{x,0^ dp da S{x- f3-aO Fi{l3,a) , 

J-l J-1+\I3\ 

and an analogous formula for the GPD E''{x,£^): 

/•I /-i-l/^l 

El^{x,0= dp / da 5{x- 13- aS) K^iP^a) . 
Obviously, the double distribution function F''{P,a) should satisfy the condition: 



da F'^{x, a) — q{x) 



(17) 



(18) 



(19) 



l+\x\ 
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in order to reproduce the forward limit (g) for the GPD H''{x,£,). It is easy to check that 
the GPDs obtained by reduction from the double distributions satisfy the polynomiality 



but always lead to /i^^]-) 



-N+1 



0, i.e. the highest power of ^ for 



conditions 

the singlet GPDs is absent. In other words the parametrization of the singlet GPDs in 
terms of double distributions is not complete. It can be completed by adding the so-called 
D-term to Eq. (0) |Pol99b| : 



dp 



1-1/31 



da S{x-l3-a^) F«(/3,a) +i 



1 + 1/31 
1 .1-1/31 

dp / daS{x-P-a^) K%P,a)-9 
-1 J-i+\0\ 



„21 



£19 



.(20) 



Here D^(z) is an odd fimction (as it contributes only to the singlet GPDs) having a support 
— 1 < z < 1. In the Mellin moments, the D-term generates the highest power of ^: 



- 9W 



dz D^z) . 



(21) 



Note that for both GPDs H'^{x,£^) and E'^{x,£^) the absolute value of the D-term is the 
same, it contributes to both functions with opposite sign. The latter feature follows from 
the relation (p^). We shall see in the section 3.2 that estimates of the D-tcrm in the chiral 
quark-soliton model gives 13" (z) « D'^{z). Therefore in the calculations below we shall 
assume that: 



D^iz) 



Nf 



D{z). 



(22) 



where Nf is the number of active flavors and D{z) — D'^{z) is the flavor singlet D-term. 

The D-term evolves with the change of the renormalization scale according to the 
ERBL evolution equation | Efr80 , Lep79| . Hence it is useful to decompose the D-term in a 
Gegenbauer series (eigenfunctions of the LO ERBL evolution equation) : 



D{z) 
D\z) 



(1-z^) 
3 



di Cl'\z) + d3 Cl^^iz) + dg C^^^iz) + 
df + d^C',^'{z) + dfc!/\z) + . 



■-{l-zr 



(23) 
(24) 



where D^{z) represents the gluon D-term. Because the D-term is a singlet quantity there- 
fore the quark D-term D{z) is mixed under evolution with the gluon D-term D^{z). Asymp- 
totically (for a renormalization scale /.t oo) we obtain: 



Das{z) 

Disiz) 



Nf 
Nf + 4Cf 

ACf 3 
Nf + 4Cf 4 



(l-z^) Cl'\z), 



(l-z^)^ 



(25) 
(26) 
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where Cf — 2N • Furthermore, the scale independent constant d is a nonperturbative 
parameter characterizing the D-term at a low normalization point, d = di + df with di 
from Eq. ( p3| ) and df being the first Gegenbauer coefficient of the gluon D-term (p^). We 
see that the D-term survives in the limit /i — > oo and therefore the complete fornirTof the 
singlet quark GPDs at an asymptotically large scale ii^ is the following: 



Nf 



Nf + iCf 



5 



(27) 



The asymptotic form of the corresponding gluon GPD is the following: 



4Cf 



Nf + ACf 



4e 



1-^ 
e 



Analogously, one can obtain the asymptotic form of the singlet GPD E{x,£^): 



(28) 



Nf 



Nf + ACf 



3/2 I X 



(29) 



The asymptotic form of the corresponding gluon GPD is the following: 



4Cf 



M 



Nf + iCF 4 



e 1- 



e 



1 



(30) 



Furthermore, note that asymptotically the GPD E{x,£,) is completely determined by the 
D-term. Note that all expressions for the asymptotic singlet GPDs depend on the scale 
independent (conserved) constant d. From this point of view this constant is as funda- 
mental as other conserved characteristics of the nucleon, such as the total momentum or 
total angular momentum. 

Up to now we considered the polynomiality properties of the GPDs H{x,S,) and 
E{x,£,). For the quark helicity dependent GPDs H{x,£,) and E{x,£,) the polynomiality 
conditions are very similar to those for H(x,^) and E{x,£)^ see Eq. ([l^). The only differ- 
ence is that in the case of the quark helicity dependent GPDs the highest powers of the 
polynomial in ^ is iV — 1 for the singlet case and N for the nonsinglet one. This implies 
that the D-term is absent for the GPDs H{x,£) and E{x,£). 



2.3 Angular momentum sum rule 

The second Mellin moments of the quark helicity independent GPDs are given by the 
nucleon form factors of the symmetric energy momentum tensor. The quark part of the 
symmetric energy momentum tensor is related to the quark angular momentum operator 
by: 

Jl ^ e'^^ J d^x x' T°'= . (31) 

*In the literature, to the best of our knowledge, only an incomplete form of the asymptotic singlet 
GPDs was presented. 
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This relation implies that the forward nucleon matrix element of the angular momentum 
operator can be related to the form factor of the symmetric energy momentum tensor. 
This results in the sum rule relating the second McUin moment of the GPDs to the angular 



momentum carried by the quarks in the nucleon |Ji97t]: 



1 

dx X {Hiix, + E^x, 0) = 2J« . (32) 

-1 

In Eq. ( |32| ) J"^ is the fraction of the nucleon angular momentum carried by a quark of 
the flavor q (i.e. the sum of spin and orbital angular momentum). The closest analogy 
for this sum rule is the relation between the magnetic moment of the nucleon and the 
eletromagnetic (e.m.) form factors: 



f,N^Fr{0) + F.nO). (33) 

In this case the magnetic moment is given by the "forward" form factor (0) plus the 
"non- forward" form factor ^2(0)^. In analogy to Eq. ( |3l| ) the relation between the mag- 
netic moment operator and the e.m. current: /i' cx e'-''^ J d^x x^ Jem contains explicitly 
the coordinate operator x^ . A detailed derivation and discussion of the angular momentum 
sum rule d32) can be found in Refs. |ji97b| , [li98b|] . 



Let us discuss here the role played by the D-term in the angular momentum sum 
rule. Unfortunately it is very hard to find an observable in which the GPDs H'^{x,£,) and 
E'^{x,^) enter as a sum. Also we should keep in mind that in a hard exclusive process, 
we have kinematically ^ 7^ and that the GPDs H{x,£,), E{x,^) enter observables with 
different kinematical factors. Therefore, we shall discuss the angular momentum sum rule 
(p3) separately for the functions H and E. We rewrite Ji's sum rule (p3) in the following 



equivalent way [ |Ji97a 



f dxxY] HHx, = -^^d^e, 

J^^dxxY,E''{x,0 = (2jQ-M^)-^die- (34) 

Here di is the first Gegenbauer coefficient in the expansion of the D-term (p3|), J'^ = J'' 
and M2 is the momentum fraction carried by the quarks and antiquarks in the nucleon: 



q ' 

We see that the sum rule (|3^) is sensitive only to the combination 2 J'5 — . What 
can we say about this combination? Asymptotically one has that lim^^oo 2 J'^ — M2 = 



|Ji96|. At some finite scale, the value of 2 J'' — M| which one would like to extract from 



the data should be compared with the contribution of the D-term on the rhs of the sum 
rule (p4) . We shall see below that estimates in the chiral quark-soliton model for the first 



15 



Gegenbauer coefficient of tlie D-term give tlie value di « —4, wliicli leads at ^ = 0.15 to 
a value of 0.07 for the second term in Eq. (p^. Therefore at accessible values of ^ the 
extraction of the angular momentum carried by the quarks from the observables requires 
(among other things) an accurate knowledge of the D-term. The size of the D-term can be 
determined from the observables which are sensitive to the real part of the amplitude of the 
corresponding process. An example of such an observable is the DVCS charge asymmetry 



(accessed by reversing the charge of the lepton beam), see a detailed discussion in Sec. 4.3.4 



2.4 Chiral properties of GPDs 

Here we shall argue that the spontaneously broken chiral symmetry of QCD plays an 
important role in determining the properties of the generalized parton distributions. We 
shall illustrate this on examples of the D-term and the GPD E{x, ^, t). 



2.4.1 Spontaneously broken chiral symmetry and the D-term 

First we will discuss how the physics of spontaneously broken chiral symmetry plays an 
important role in determining the size and the sign of the D-term. In particular, this role is 
seen clearly in the case of the GPDs in the pion. In this case the value of the coefficient di 
in the parametrization of the pion D-tcrm ( p3| ) can be computed in a model independent 
way and it is strictly nonzero. To compute the pion D-term we use the soft-pion theorem 
for the singlet GPD in the pion derived in |Pol99a|. This soft-pion theorem has been 
derived using the fact that the pion is a (pseudo)Goldstone boson of the spontaneously 
broken chiral symmetry. The theorem states that the singlet GPD in the pion vanishes for 
^ = ±1 and momentum transfer squared t = (corresponding to to a pion with vanishing 
four momentum): 



^ilW(x,e = ±l) = 0. 



(36) 



Hence (see also [ Pol99b| ), 



dx : 



(37) 



Evaluating Eq. (^ at ^ = 0, determines : 



Jo 



x)+q^''Hx)] 



(38) 



being the fraction of the momentum carried by the quarks and antiquarks in the pion. As 
the highest power in ^ in Eq. ( ^7|) (i.e. the term in ^^) originates solely from the D-term, 
one easily obtains the following expression for pion D-term : 



i?W(z) = -^^^2 



(l-z2) 



(39) 



We see therefore that the first Gegenbauer coefficient of the pion D-term is negative and 
strictly nonzero. Its value is fixed by the chiral relations in terms of the momentum fraction 
carried by the partons in the pion. 
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The nucleon D-term is not fixed by general principles. However one may expect that 
the contribution of the pion cloud of the nucleon can be significant. Indeed in the chiral 
quark-soliton model, which emphasizes the role of broken chiral symmetry in the nucleon 
structure, the D-term is large and has the sign of the pion D-term. Quantitatively for the 



coefficients ^1,^3, ds (see Eq. 23), the estimate which is based on the calculation of GPDs 
in the chiral quark soliton model |Pet98| at a low normalization point /x ~ 0.6 GeV, gives 
[ KivOlbl (see also Sec. pTl] ) : 



di « -4.0, 



-1.2, 



-0.4, 



(40) 



and higher moments (denoted by the ellipses in Eq. (^3|)) are small. Notice the negative 
sign of the Gegenbauer coefficients for the D-term in Eq. (^), as obtained in the chiral 
quark-soliton model. The coincidence of the signs in the nucleon and pion D-terms hints 
that the D-term in the nucleon is intimately related to the spontaneous breaking of the 
chiral symmetry. The chiral contributions of the long range two-pion exchange shown in 
Fig. ^ to the GPDs of the nucleon are sizeable giving large and negative contribution to 
the leading Gegenbauer coefficient of the nucleon D-term di. The sign is related to the 
negative sign of the pion D-term which is fixed by the soft pion theorem (^) . 



Figure 3: Two-pion exchange graphs giving the chiral contributions to the nucleon form factors 
of the energy momentum tensor. The dashed lines represent the pion, and the solid lines represent 
the nucleon. The black blob is the operator of the energy momentum tensor. 



2.4.2 Pion pole contribution to the GPD E{x,^,t) 



The role of the spontaneously broken chiral symmetry in the structure of GPDs is seen par- 
ticularly clearly in the case of the quark helicity dependent GPD E{x, ^, t). We remind that 
this GPD satisfies the sum rule of Eq. ([ll|), in terms of the pseudoscalar nucleon nucleon 
form factor /i^(t). It is well known (see e.g. Ref. [Adl68|) that due to the spontaneously 
broken chiral symmetry this form factor at small t is dominated by the contribution of the 
pion pole of the form: 



where gA ~ 1.267 is the nucleon isovector axial charge and is the Pauli matrix in the 
flavor space. The presence of the chiral singularity on the rhs of the sum rule (^) implies 
that one should also expect the presence of the chiral singularity in the GPD E'^{x,£^,t) 
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at small momentum transfer squared t |Fra98a|. The form of this singularity has been 
specified in Refs. [Man99a, Fra99a, PenOOa to be: 



hm E'^ix £t) = -r' ^^^"^^ 



1 ^ 



(42) 



where is the pion distribution amplitude entering e.g. description of the pion e.m. 



form factor at large momentum transfer and the hard reaction 7*7 — s- tt . In Ref. | PenOOa 
in the framework of the chiral quark-soliton model Eq. (^) and the deviations from it 
have been computed. The corresponding results are discussed J.n Sec. |3.2| . 

The presence of the pion pole singularity in the GPD E leads, in particular, to a 
strong dependence of the differential cross section of, for example, hard 7r+ production on 
the transverse polarization of the target | Fra99a | . The dependence of the exclusive cross 



section on the transverse polarization of the proton target has the following dominant 
structure ^ 



a CX 



15, |sin$ 



— t + m'i 



(43) 



where sin $ is the azimuthal angle between lepton plane and the plane spanned by 3- vectors 
of the virtual photon and produced meson. The specific t-dependence of the cross section 
is due to the contribution of the chiral singularity (pion pole) to the "truly non-forward" 
GPD E whose presence is dictated by the chiral dynamics of QCD. 



2.5 Twist-3 GPDs 



Recently the DVCS amplitude has been computed in Refs. |AniOO, PenOOb, BelOOb, RadOC, 
RadOlbj including the terms of the order 0{1/Q). The inclusion of such terms is manda- 
tory to ensure the electromagnetic gauge invariance of the DVCS amplitude to the order 
V^/Q. At the order 1/Q the DVCS amplitude depends on a set of new generalized parton 
distributions. Let us introduce generic generalized "vector" T^(x,£) and "axial" T^{x,£^) 
distribution^: 



ZTT 



dy- e"^^y'{p'm-y/2)ri:{y/2)\p) 



T^\x,^,t) 



p+ r°° -+ _ _ 

— / dy- e"^ y-{pyj{-y/2)rj,i^{y/2)\p) 



(44) 



y+=j/j_=o 



If in the above equations the index fi is projected onto the "plus" light-cone direction (i.e. 
H = +) we reproduce the twist-2 GPDs, see definition (||). The case of /i = — , according 

^Detailed expressions and quantitative estimates are shown in Sec. [4.5.2 



-^In this section wc do not specify the flavor index as all formulae here are valid for each flavor separately. 
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to the general twist counting rules of Ref. | Jaf91| corresponds to the twist-4 GPDs which 
contribute to the DVCS amplitude at the order of 0(1/Q^). The twist-3 GPDs correspond 
to /i =_L. For a detailed discussion of the DVCS observables to the twist-3 accuracy see 



Sec. 4.3.3. Here we first discuss the general properties of the twist-3 GPDs and then will 
show that they can be related to the twist-2 GPDs in the so-called Wandzura-Wilczek 
approximation. 

2.5.1 General properties of the twist-3 GPDs 

Several properties of the twist-3 GPDs can be derived without invoking any approximation. 
Forward limit : In the forward limit the "vector" GPD can be parametrized as follows: 

Uin 2[g(a;)p^ + m%fi{x)n^] , (45) 

where q{x) is the unpolarized forward quark (antiquark for x < 0) distribution of twist- 
2, f4{x) is the twist-4 quark distribution, see [[Jaf91 |. Also we introduced the light-cone 



vectors and n^, the former has only "plus" non-zero component j3+ = P+ and the 
latter only "minus" component — 1/P+. As it could be expected, the twist-3 part 
of J^^ disappears in the forward limit, because unpolarized parton densities of twist-3 are 
absent. 

The forward limit for the twist-3 part of the "axial" GPD Tj^J^ is nontrivial: 



^mT^A^x,0 = 2S^ATqix). 



Here ATq{x) is the transverse spin quark distribution which can be measured in DIS on a 
transversely polarized target. 

Sum rules : Performing the integral over x of the distributions JF^(a;,^) and J-^{x,£^) given 
by Eq. jii] ) we obtain on the rhs the nuclcon matrix elements of the local vector and axial 
currents. This results in the sum rules: 



' dx T^{x,^,t) = N{p ')^,N{p)F,{t) + N{p ')'.^f^N{p)F^{t) , 
1 ~ - A - 

^dx T,Xx,i,t) = N{p')^^-i^N{p)gA{t) + ^N{p')^5N{p)hA{t). (46) 
These sum rules can be viewed as a generalization of the Burkhard-Cottingham sum rule 



for the transverse polarized distribution Arg(a;) |BurC7C|. Indeed, in the forward limit 
the sum rule ( p6| ) for twist-3 GPD JF^^^ corresponds directly to the Burkhard-Cottingham 
sum rule 



/ dxATq{x) ~ / dxAq{x). 
Jo Jo 
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The second Mellin moment of GPDs T^^{xX) and JF^(a;,^) can be also computed in 



terms of moments of twist-2 GPDs without invoking any approximation |PenOOb, |Kiv01c 



Indeed taking the x-moment of Eq. ([44[) we obtain on the rhs the local operator of the type 
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■07^(75) V '^ip which with the help of QCD equations of motion can be reduced to either 
the local operators of twist-2 or their total derivatives. Let us illustrate how this works on 
the example of the twist-3 part of the "vector" GPD ^^^^(a;,^). From the definition 
we obtain: 



1 ^ 



(47) 



Here Vq, is the gauge covariant derivative in the fundamental representation. We rewrite 
the operator in Eq. (E^) identically as: 



V'7± V +0 



1 
2 
1 

+ 2 



./3 Y7 a. 



(48) 



The first bracket is the operator of the twist-2 (it is symmetric and traceless) and can 
be expressed in terms of the Mellin moments of the twist-2 GPDs H, E. The operator 
in the second bracket can be reduced to the total derivative of the axial current with the 
help of QCD equation of motion for a massless quark: y'0 = 0. Indeed, starting with the 
identities of the type "0707/3 VV" = one derives the following identity: 



0(x) 



la V/3 -7/3 Va 



0(a;) = iEafjpa d° 



tlj{x)Yl5i^ix) 



(49) 



which can be used to reduce the second bracket in Eq. (^7|) to the total derivative of the 
axial current. The latter is obviously the operator of twist-2. 

Puting everything together we obtain the following sum rules JPenOObl , |Kiv01a| : 



dxxT^Jx,^,t) = N{p')^'^N{p) 



dx X {H{x, t) + E{x, t)} + gA{t) 



+ ^^(P')^(^^)^ J ^dxx E{x,C,t). 



We see that the GPD JT^^ is sensitive to the combination 



(50) 



dx X {H{x, t) + E{x, t)) + gA{t) 



which is related to the spin structure of the nucleon in the forward limit: 
lim [ dxx{H'i{x,tt)+E'i{x,£„t)}+9A{t) = 2J'i + Aq, 

A^O J _i 



(51) 
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where J' is for a given flavor a fraction of the total angular momentum of the nucleon 
carried by quarks and Ag is the corresponding fraction of the spin. To derive Eq. we 
made use of angular momentum sum rule (32). We see that the twist-3 DVCS observables 
can be used to probe the spin structure of the nucleon. In principle, this part of the 
amplitude can be extracted from the data through the angular, spin and Q dependence 
of the differential cross section Die97, BelOOa, BclOlb]. As will be discussed further in 
Sec. 4.3.4 and Sec. 4.3.5, this may provide us with additional possibility to probe J'^ in 



the nucleon. 



In the case of the axial transverse GPD J^^^ wc obtain in the same way the following 
result for the corresponding second moment: 



dx X T^^{x,^,t) - ^(P ')7^75iV(p) - / dxxH{x,0+e{Fi{t)+F2{t)} 



other structures . 



(52) 



For simplicity we do not write all possible spinor structures which appear on the rhs. The 
sum rule ( |5^ ) can be viewed as the generalization of the Efremov-Leader-Teryaev sum rule 
[Efr97| for the polarized structure function ATq{x). The reason for such sum rules is the 
absence of the genuine twist-3 local gauge invariant operators of the form •ipGip (G is the 
gluon field strength) having dimension four. The third and higher Mellin moments are not 
generically related to the moments of twist-2 GPDs H, E and i?, they get additional 
contributions from the genuine twist-3 local gauge invariant operators of the form rpGip. 
For instance the third Mellin moment of the GPDs JF^ and JF'' has the following structure: 



dx x^ J^^(J^^) = 



1 



r(p'lW(75)5G+'^V^|p) 



3(P+)2 
moments of H, E,H,E . 



(53) 



Here G'^" = |e 



iivafiQ^^ is a dual gluon field strength tensor. Generically the value of the 
matrix element of the genuine twist-3 V'G'0 operator in the above equation is comparable 
with the size of matrix elements of twist-2 operators. The genuine twist-3 ma trix element in 



Eq. (53) in the forward kinematics has been estimated using QCD sum rules [Bra87 



3te95 |, 



in the instanton model of the QCD vacuum | Balla98 | and recently in lattice QCD | Goc01 |. 

In the instanton model of the QCD vacuum the matrix elements of genuine twist-3 
operators are parametrically (and numerically) suppressed relative to the corresponding 
matrix elements of twist-2 operators [ Balla9S| , DresOOaj. This suppression is due to the 
small packing fraction of instantons in the QCD vacuum which is related to the ratio of the 
average size of the instantons to their average relative distance: ~ \- The dependence 
of the twist-3 to twist-2 parts of forward matrix elements on this small parameter was 
obtained in Refs. Balla98|, DresOOa| as: 



twist — 3 p"* f 
twist - 2 ^ W ^^Vb? 



(54) 



Several indications for the smallness of the above ratio have been provided recently by mea- 



surements of the transverse spin structure functions at SLAG [Mit99, BosOC|. We see that, 
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at least for the forward matrix elements, the genuine twist-3 operators can be neglected 
relative to the twist-2 ones. If one does so, one can express the twist-3 GPDs in terms 
of twist-2 ones. Such approximation is called Wandzura-Wilczek (WW) approximation 
| Wan77| . 

The twist-3 GPDs in the Wandzura-Wilczek approximation are discussed in the next 
section. 



2.5.2 Wandzura-Wilczek approximation for the twist-3 GPDs 

The aim of this section is the calculation of the matrix elements (Q) in the Wandzura- 
Wilczek (WW) approximation. The operators entering the definition of the GPDs and 
J-^ of Eq. ( ^ ) can be decomposed in a gauge invariant way into two parts, one of which 
is fixed by the twist-2 GPDs iJ, E, H and E, and another - the "genuine twist-3" part - 
associated with operators of the type ipGijj. The latter operators are interaction dependent 
and measure quark-gluon correlations in the target hadron. The WW approximation 
consists in the assumption that the nucleon matrix elements of the interaction dependent 
operators are small and can be neglected. 



To obtain the WW relation we make use of the operator identitiesf^ derived in |BelOOb 
KivOlc, RadOO, Kiv0l4 Rad01a| on the basis of technique developed in Ref. |Bal8£|: 



3l. 

{Vx) 



'ijj{x):^ip{—x) 



(xd) 



" (xV) 
^e^^^fipX^j^dP 



2 7o 



da 



|ga(xx>) ^ e-"(^^)} i/.(a2;);/<i/;(-aa;) 

a |e-"(^^) - e"(^^)} V;(ax)9<75V'(-aa;) + 



(55) 



where a = 1 — a and ellipses stand for the contributions of cither twist-3 quark gluon 
operators or twist-4 operators which both are neglected in the WW approximation. An 
analogous expression for the operator %l){x)^ ^^^'ip{—x) can be obtained from Eq. ( |55| ) by 
the replacement ^ ^75 on rhs of this equation. The symbol T) denotes the derivative 
with respect to the total translation: 



d - 

Pa {tl){tx)T[tx, -tx\il;{-tx)] EE — {V'(te -f y)T[tx + y, -tx + yji^i^tx + y)} 



(56) 



with a generic Dirac matrix structure F and [x, y] = Pexp[igJ^dt (x — ?/)pA^(te+ (1 — t)y)]. 
Note that in the matrix elements the total derivative can be easily converted into the 
momentum transfer: 



(p '\Vf,ii;{tx)T[tx, -tx]il;{-tx)\p) ^ i{p ' - p)^{p '\'iJj{tx)T[tx, -tx]'iPi-tx)\p) (57) 
and therefore in the matrix elements one can associate V with the momentum transfer A. 



The identity (55) allows to express the matrix elements of the non-local "vector" 
and "axial" operators {Ihs of (^5|)) in terms of matrix elements of symmetric operators 
generated by non-local operators V'(x)^(75)?/'(— a;) standing on the rhs of Eq. ( |55| ). Since 



Alternative derivation based on generalized rotational invariance can be found in Ref. [AniOl 



22 



eventually (after performing differential operations on the rhs of Eq. (|55|)) the coordinate 
x'^ is put on the light-cone the symmetric operators on the rhs are reduced to the twist-2 
operators and their total derivativesm. This allows us to obtain the WW-relations for 



twist-3 nucleon GPDs [BclOOb, KivOl 



tive£^. 
[ivOla] : 



:N{p ')Nip)Eix,0 - ^NiP 'Wip){H + E){x,0 



2^ 



+ / duG^{u,OW+ix,u,0 + i^±^.k duG''{u,OW-{x,u,0,m 



T^^{x,0 



N{p ')j,N{p)E{x,0 - ^^{V ')fij5N{p)H{x,0 



2miy 
1 



2^ 



+ / duG^{u,OW+{x,u,0 + ii±f.k duG''iu,OW-ix,u,^).{59) 



The following notations are used : 
G^iu^O 



N{p 'h'[N{p)iH + E){u,0 + TTT^Nip ')N{p) 



du dS, 



-^N{p ')fiN{p) 



du dS, 



{H + E){u,0 



E{u,0 
(60) 



Nip ')7^75A^(p)i?(",0 + ^Nip 'h^Nip) 

IrriN 



du d£, 



E{u,0 



A'' 

-^Nip ')fr/5N{p) 



'u—+£,— 
du dS, 



Hiu,0^ 



(61) 



The functions W±{x, u, £,) can be called Wandzura-Wilczek kernels. They were introduced 
in Ref. [BelOOt, Kiv01a| and are defined as : 



/ ,\ 1 f^/ ^\ ^(u > x) ^, ^.9(u<x) 
W±{x,u,0 = 7:\Six>0- -^-Oix<0- ^ 



2 1 u + S, 



ix<-0- 



[u < x) 



(62) 



We also introduce the metric and totally antisymmetric tensors in the two dimensional 
transverse plane (eoi23 = +1) : 



a~f3 



■pn p 



(63) 



^■^ Note that in the pioneering approach of Ref. piuOO | to the WW relations for GPDs the total derivatives 
were neglected. In this paper the reader can find a detailed discussion of the mathematical aspects of the 
twist decomposition of light-cone operators. 
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In Sec. 4.3.4 and Sec. 4.3.5 we shall use the WW-relations (58 5^) in order to estimate 
the contributions of the twist-3 effects to various DVCS observables. 

In the rest of this section, we discuss properties of the twist-3 GPDs peculiar for the 
WW approximation. Let us start with a discussion of the properties of the WW-kernels 
( |6^ ) . We introduce the following notations for the action of the WW kernels on a function 



W±(^f[x,^]= duW±ix,u,Ofiu,0, (64) 



-1 



with W± given by Eq. ( |6^ ) . We shall call "WW transform" the resulting functions W± (8) 
f[x,£] and we shall call the "WW transformation" the action of the WW kernels. 

Limiting cases. We consider two limiting cases of the WW transformation: the forward 
limit ^ ^ and the 'meson' limit ^ — > 1. In the forward limit we easily obtain: 

\imW+®f[x,^] = 0(:r>O) /(u,e = O)-0(x<O) r — /(u,C = 0), 

f^o U J_i u 

\miW-® f[x,^] = 0. (65) 
We can see that the action of the W+ in the forward limit reproduces the Wandzura- 



Wilczek relation for the spin structure function |Wan77|. The term with 9{x > 0) 
corresponds to the quark distributions and the term with 6{x < 0) to the antiquark 
distributions. The W- kernel disappears in the forward limit, so that this kernel is a 
'genuine non-forward' object. 

In the limit ^ — > 1 the generalized parton distributions have properties of meson 
distribution amplitudes. In this limit the WW transforms have the form (we use the 
notation f{u,£, = 1) = ^^(m)): 

„r ^1 1 { du , ^ du 

iml¥±®/ x,C = - / ^{u ± T— ^ « 

C^i 2 {J_i 1-u 1 + u 

which corresponds to the WW relations for the meson distribution amplitudes derived in 



|Ball96, Ball9g 



We can see that WW transforms of generalized parton distributions interpolate be- 
tween WW relations for parton distributions and for the meson distribution amplitude. 
The general form of the WW kernels (^2|) also allows us to derive WW relations for dis- 
tribution amplitudes of a meson of arbitrary spin. 

Mellin moments : One can easily derive the Mellin moments of the WW transform. The 
result has the form: 



^dxx'' W±(g>f[x,^] = — 



1 



du 



,N+1 



cN+l 



,N+1 



± 



i-0 



N+1 



u + £, 



/(u, 0.(66) 



From this simple exercise we can see an important property of the WW transformation, 
namely, if the function f{u,^) satisfies the polynomiality condition, i.e.: 
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du 



.N 



f{u, ^) = polynomial in ^ of the order + 1 



(67) 



its WW transform also satisfies the same polynomiality condition. 



Discontinuities . In Refs. [ KivOlc , RadOC , RadOlaj it was demonstrated that the twist-3 
skewed parton distributions in the WW approximation exhibit discontinuities at the points 
X = This feature is related to the properties of the WW kernels. Let us compute the 
discontinuities of a WW transform at the points a; = 



lim \W± 



hm \W± 



1 



du 
du 



/(^,e) 



(68) 



Here vp means an integral in the sense of valeur principal. We see that for a very wide 
class of functions /(u, ^), the discontinuity of the corresponding WW transforms is nonzero. 
This feature of the WW transformation may lead to the violation of the factorization for 
the twist-3 DVCS amplitude. 

Using general properties of the WW transformation (pq) and Eqs. (p8|,p9|) one obtains 



that J^^^ and J^J^'^ have discontinuities at the points x = ±^ . But using a certain sym- 



metry of the Eqs. (58) and (p9|), one can find that some combinations of the distributions 
T^"^ and T^"^ are free of discontinuities. For example, using Eq. ( |68| ) one can see that 
the combination: 



has no discontinuity at x 



^. On the other hand, the 'dual' combination: 



WW 



(69) 



(70) 



is free of the discontinuity at a; = — ^. The cancellation of discontinuities in these particular 
combinations of the GP Ds ens ures the factorization of the twist-3 DVCS amplitude on the 
nucleon, as given in Sec. 4.3.2 , 

Pion pole contribution: The generalized parton distribution E at small t is dominated by 
chiral contribution of the pion pole, see Eq. (^). In the WW relation (58) the twist-2 
GPD E enters only in the combination: 



1 



d 



d 



du dS, 



E{u,0- 



(71) 



One can easily see that the contribution of the pion pole ( ^ ) nullifies under the action 
of the differential operator in Eq. (£l]) . The only place where the pion pole contribution 
survives is the GPD J^^, see Eq. (p9|). In this way we obtain the following simple results 
for the contribution of the pion pole to the twist-3 GPDs in the WW approximation: 



^WW,piQn Polc^^ 



0, 

A 

2™ 



N 



-Nip')^,Nip)- 



1 



'AT 



(72) 
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We see an interesting result that the pion pole contribution to the twist-3 GPDs is expressed 
in terms of twist-2 pion distributions amplitude. This observation is in nice agreement with 
the fact that in the WW approximation there exists no twist-3 pion distribution amplitude 



associated with vector- and axial- vector operators |Ball96 



D-term contribution: Now we discuss the D-term contribution to the twist-3 GPDs. If we 

|20|) for the GPDs H and E into the WW relations 



substitute the D-term contributions ^ 
of Eqs. dUg]) we see that the result is 



WW D-torm 



WW D-torm 



-N{p ')N{p)D ' 



2^ rriN 



0{\x\<0- 



(73) 



In some sense the D-term is an eigenfunction of the WW transformation. 



2.6 GPDs for ^ A, E transitions : SU(3) relations and beyond 

The studies of hard exclusive processes with strangeness production give an access to the 



flavor non-diagonal generalized parton distributions [Fra99b|. These distributions corre- 
spond to the process where a quark of one flavor is taken out of the initial nucleon at 
the space-time point y/2, and then a quark with another flavor (say a strange quark) 
is put back exciting a hyperon at the space-time point —y/2. The strangeness changing 
distributions ior N ^ Y transitions {Y = S, A) can be defined as (cf. Eq. (^)): 



2tt 



dy e 



ixP+y / I 



{p\Y\-S0{-y/2)q^{yl2)\p,N) 



+(757 )aP 



H 



{x,tt) Y{p)^+N{p) 



E 



{x,^,t)Y{p)ia 



A. 



2mN 



N{p) 



H 



N~*Y 



(a;,^,t)y(p)7+75iV(p) 



E 



N^Y 



{x,^,t)Y{p')j5:^N{p) 



(74) 



where sp{—y/2) is the field operator of a strange quark, q{y/2) is the field operator of a 
nonstrange quark (with q = u or d depending on the charge of the hyperon Y) . 

The flavor non-diagonal GPDs provide a new tool to study the non-perturbative 
structure of various nucleon-hyperon transitions. Below we shall relate strangeness chang- 
ing ^ F GPDs to ^ iV GPDs with help of flavor SU{3) relations. But before doing 
this let us discuss some interesting features of the N Y GPDs which are related to the 
SU{3) symmetry breaking effects. 

Let us start with a discussion of the form factor sum rules for strangeness changing 
N ^ Y GPDs. Integrating the Ihs of Eq. ( [Z^ over x one gets the local sq operator 
sandwiched between a nucleon and a hyperon state. Such matrix elements are parametrized 
in terms of vector and axial N Y transition form factors entering the description of 
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semi-leptonic decays of hyperons (see e.g. | Gar85{ ) p^ : 



{p\Y\srq\p,N)=Y{p') 



2mN ^ 2mN 



Nip), 



{p\Y\sn5q\p,N)=Y{p') 



2mN 



2m 



N 



75iV(p). (75) 



We see that for the strangeness changing transitions two additional form factors f^^^ (t) 
and 92 {t) related to the second class currents (52 is also often referred to as weak 
electricity ) appear. Both these form-factors are proportional to flavor S't/(3) breaking 
(~ TOs) effects and they are absent for transitions without strangeness change. 

The form factor sum rules for strangeness changing GPDs have the following form: 



~ ^ 2mN ^ ' ' 

/.^-^W+e/3^-^(i), 



dx H'^-^"" {x,^,t) 



9i 



dx E^-^^{x,^,t) 



33^-^W + ^5f-^(0. 



(76) 



One sees that in contrast to the corresponding sum rules for the GPDs without strangeness 
changing (see Eqs. (|§-0)) the sum rules ( [76| ) get additional contributions from the SU{3) 
symmetry breaking form factors /^(t) and 32(^)1 associated with the second class currents. 
Note that these additional contributions also violate the symmetry of GPDs under the 
transformation ^ — ^ [ Man98a| . The reason is that for the derivation of this symmetry 
the time inversion transformation has been used. Under this transformation the initial 
and final baryon states are exchanged. In the case when this transition occurs between 
baryons belonging to the same symmetry multiplet the T-invariance gives the symmetry 
property under ^ — > — ^. However for the case of transitions between baryons from dilferent 
multiplets the T-invariance is not restrictive (it gives only th e realn ess of the GPDs). 

Studies of semiletponic hyperon decays showed (see e.g. [ Gar85| ) that the assumption 
of the SU{3) flavor symmetry works very well for form factors 51, /i, and f2- Unfortunately 
the semileptonic decays of hyperons are not sensitive to /a, (72 and 53. On general grounds 
we can state that the form factors fs and (72 are proportional to the explicit SU (3) breaking, 
i.e. they are proportional to the strange quark mass, therefore their effect can be non- 
negligible. For the form factor (73 (t) we can apply the generalized PC AC relations to obtain 
the contribution of the kaon pole: 



,3---(t)«^^^, (77) 

TOjf — t 

^^We change the definitions of the form factors /2,92 and 33 by a factor of two to have correspondence 
with the standard definition of the nuclcon form factors. 
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where qkny are the KNY coupling constants, and Jk ~ 159 MeV is the kaon decay 
constant. Comparing the i— dependence of this equation with the one of the the pion pole 
contribution to pseudoscalar nucleon form factor (|4^) we notice a big difference due to 
the large difference between the kaon and pion masses. This shows that SU (3) symmetry 
breaking effects are sizeable for g^^^it). 

From the sum rules (^6|) we see that the SU (3) symmetry breaking effects are sup- 
pressed for the GPDs H, E, and H, at least they are suppressed in the sum rules by either 
^ (like for H and E) or my — tun (like for H). On contrary, the sum rule for GPD 
E^~*^ {x, ^, t) shows that the flavor symmetry breaking effects can be very strong for this 
GPD. Indeed the first term on the rhs of the sum rule is proportional to g^~~*^{t) for 
which one expects large symmetry breaking effects (see discussion above). Additionally 
the weak electricity (symmetry breaking) form factor enters the corresponding sum rule 
with the factor of 1/^. This, in particular, implies that in the limit ^ — » the function is 
divergent^, and this divergency is dominated by weak electricity. This observation opens 
an exciting possibility to study the weak electricity of nucleon-hyperon transitions (oth- 
erwise hardly accessible) in hard exclusive strange meson production. In this respect, the 
transverse spin azimuthal asymmetry in the production of kaons on the nu cleon is a very 
promising observable as it is sensitive to the function E^^^ , see Sec. 4.5.4 and qualitative 
discussion below. _ 

Given the large expected SU{3) symmetry breaking effects for the GPD E^^^ we 
model this function by the contribution of the kaon pole plus the weak electricity contri- 
bution: 



E 



N->Y 



m 



1. 



K 

N^Y 
WE 



t 



i*„.(i 



(78) 



where ^k{z) is the kaon distribution amplitude and E^l^^ {x, t) is the "weak electricity" 
GPD normalized to the corresponding weak electricity form factor 5^~*^'(i): 



dx E'lYj^ I 



x,^,t) 



(79) 



This GPD is proportional to explicit SU{3) symmetry breaking effects (~ mg). If for the 
weak electricity GPD we use the following rough (asymptotic) ansatz E^~^^ = 3/4(1 — 
2^^/'?^)/^ 92 ^^{^) transverse spin asymmetry in production of the kaons has the 
dominant contribution (c/. Eq. (E^)): 



N^YK 



cx 



|5j_|sin$ ^/U, 



t 



gKNYrUNfK 



9t 



'{t) 



-t + ra\ 2^ 



H 



(C,e,i). 



^*One can introduce an "ovecomplete" set of GPDs adding to the rhs of Eq. one additional "weak 
electricity" GPD _B^;g'^(T+'' A^75/2mjv . In this way one avoids working with GPDs non normaUzable at 
f = 0. 
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This expression illustrates that at small ^ the measurement of hard exclusive kaon produc- 
tion can give us (among other exciting things!) an access to the symmetry breaking weak 
electricity form factor 5^^^ (t) . The contribution of the weak electricity is enhanced at 
small xb- For example for the j*p — > AK~^ processes at —t = 0.5 GeV^ and xb ~ 0.05 



the contribution of the weak-electricity part in (80) is of the same order as the contribu- 
tion of the kaon pole. We take for this estimate the value of 32^^(0) « 0.1 obtained in 
Ref. |Kim98| in the framework of the chiral quark-soliton model. This simple estimate 



shows that the studies of the second class currents are feasible in hard exclusive reactions 
with production of strangeness. 

For the strangeness changing GPDs jjN^y ^ j^n^y j^n~>y ^^-^ apply the 
flavor SU (3) relations as it seems that for these GPDs the effects of the sym metry br eaking 
are small. The SU{i) relations for these functions have the following form [ Fra99a |: 



H^^^{x,^,t) 

HP^''\x,U) 
HP-^\x,U) 



1 

1 

71 



(2i/" -H"" -H' 



(80) 



and similar expressions for GPDs E^^^ and H^^^ . For transitions from the neutron 
we can analogously derive: 



H 



[xX.t) 



H 



1 

1 



{2H'^ - i/" - H') 



'ix,^,t) = 



V2 



(81) 



The SU (3) relations ( pO| , pll ) give us a possibility to access strange quark distributions 
in the nuclcon. Indeed, using these relations we can derive the following obvious relations 
[we show only the example of H] : 



lim HP^^(x,^,t) 



1 
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{2Au{x) - Ad{x) ~ As{x)) , 



^mHP-^{x,tt) = -^{Ad{x)-As{x)) 
Um^ HP^^^ (x, i) = - {Ad{x) - As{x)) . 



(82) 



This illustrates that measurements of strangeness changing GPDs can give us important 
information about strange (anti-)quark polarized distributions. 



2.7 GPDs for TV ^ A transition 

In the hard exclusive processes with an excitation of the A one accesses the parton distri- 



butions for ^ A transition |Fra98a, FraOOa|. A knowledge of such distributions would 
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give us the detailed information about the physics of iV ^ A transition, in particular it 
would allow us to understand the physics of this transition at the level of partons. 

The ^ A quark distributions can be defined through nondiagonal matrix elements 
of products of quark fields at light-cone separation. We give the definitions of — > A 
generalized quark distributions for p —f A"*" transition. The GPDs of other transitions, e.g. 
p A++, p A", etc. can be obtained with the help of the following isospin factors: 

b^A++]^y|[p^A+], [p^A"] = y|b^A+], (83) 

and analogous relations for transitions from the neutron. Keeping in mind this simple 
isospin relations, we introduce for the p A+ transition the following three GPDs (cf. 
Eq. I with if = An'')): 

^e^^^A,p'm-Xn/2)fiT^ ^(An/2)|iV,p) = (84) 

Here is a light-cone vector, 

= 0, n-P = 1, P^^^{p + p'f (85) 

A is the four-momentum transfer, 

A=p'-p, (86) 

tp^{p ') is the A Rarita-Schwinger spinor, N{p) is the nucleon Dirac spinor. The N ^ A 
quark distributions IIm,e,c{x, A^) are regarded as functions of the variable cc, the square 
of the four-momentum transfer, i, and its longitudinal component 

^ = -\{n-A). (87) 

The covariants /C^f,'^'*^ are magnetic dipole, electric quadrupole and Coulomb quadrupole 



covariants |Jon73|: 

AL«„ = -i^ 77 tt; :z£au\a-r A 



_ 3(mA+mjv) 
-'^'^ " "*2mw((mA+m^)2-i)''3MAa^ 



where we introduced the notation Z{t) = [(jtia + jtijv)^ — <][(mA — niisi)'^ — t]. We choose 
such definition of the iV ^ A quark distributions in order to have the following sum rules: 



J dx HM,E,cix,S,,t) = 2 Glj E ci*) 



(89) 
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where G%j ^ (^{t) are the standard magnetic dipole, electric quadrupole and Coloumb 
quadrupole transition form factors respectively |Jon73|. 

Analogously to Eq. ( p^ ) we introduce the following quark helicity dependent p — > A+ 
generalized quark distributions: 



^e^^- {A,p'\^i-Xn/2)frr\^^{Xn/2)\N, p) = 



C,{x, t)np + C2{x, 1 1) + Cz{x, t)^{npl^ - A^^} 



"N 



rriN 



Nip) 



(90) 



Actually in this definition not all the functions Ci are independent, for example the function 
C4 can be absorbed in a redefinition of the functions Ci and C2 . We choose the "overcom- 
plete" definition ( pO| ) to have simple (^-independent) sum rules for the first Mellin moment 
of generalized distributions Cf. 



' dx C^{x,i,t) = 2 C^{t) , C dx C2{x,i,t) = 2 Ci{t) , 

-1 J ~i 
1 .1 

dxC3{x,i,t)^2C^{t), / dxCi{x,£,,t)^2Ct{t). (91) 

1 J -I 



Here C/^ (<) are axial A'' ^ A transition form factors in the notations of Adler | Adl75|| 



(A+,p'|^(0)7^75^(0)b,p) 



''N 



m, 



N 



rriN 
N[p) 



(92) 



In Sec. I we shall derive large Nc relations between the A A generalized quark distri- 
butions and the usual N ^ N ones. 

For the A^ — > A GPDs one can repeat almost word for word the discussion of such 
issues as polynomiality, pion pole contribution, etc. Howev er in applications for hard 
exclusive processes with excitatio n of A (see Sec. 14 and Sec. 4.5.4 ) we shall use only the 
large A^, relations derived in Sec. |3.1.l| . Therefore here we restrict ourselves to discussion 
of these relations. They have the following form: 



HM{x,^,t) 

Ciix,tt) 



-^[E^{x,tt)-E''{x,U)] 



H^ix,^,t)-H\x,^,t) 



^ \E-ix,U)-E''ix,U) 



(93) 
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All other N A GPDs are zero at leading order of 1/Nc expansion. Let us also note 
that if we take the first moment of the Eqs. (^3|) and take into account the sum rules for 
N ^ N GPDs (p[]ll|) we reproduce the well known result of soliton models for the ^ A 
transition form factors. Integrating Hm, given by the large Nc relation of Eq. ^ over x 
at t = gives the relation between magnetic moments: 



fJ-NA 



(94) 



which is satisfied rather well by the experimental data on magnetic moments. As we 
discussed in Sec. 2.4.2, the GPD E contains a contribution of the pion pole. From the 
large Nc relation for C2 in Eq. ( p3| ) we conclude that the A^ ^ A GPD C'2{x,^,t) also 
contain the pion pole contribution of the form: 



lim C2{x,tt) - 



9Am% 
mi — t 



1 ^ 



(95) 



The A?^ — > A distribution Hm{x, S,,t) in the large Nc limit contains information about 
the isovector part of the angular momentum of the nucleon carried by the quarks - J" — J"^. 



Indeed with help of angular momentum sum rule ( 34 ) we easily obtain 



lim [ dx X Hm (a;, ^ , t) = ^ [2 ( J" - J'^) - + M^] 



(96) 



where, as usually: 

r-l 



M| = dx X [q{x) + q{x)] . 
Jo 

Note that the contribution of the D-term does not appear in the N ^ A transition, because 
the D-term is a flavor singlet quantity. 

Obviously one can easily introduce GPDs for other A^ N* {N* is a nucleon res- 
onance) transitions which enter description of hard exclusive processes with excitation of 
baryonic resonances. Such kind of reactions can be viewed as a new tool for the study of 
baryon spectroscopy (spectroscopy by light-cone operators) . 

2.8 Further generalizations: GPDs for N — > ttN transitions 

The transitions of the nucleon to a nN system with low invariant mass can also be studied 
by light-cone probes encoded in the A^ — > ttN GPDs^^ which can be introduced as follows 
(we only show the parametrization of the "vector" bilocal quark operator): 



^ / dy-e"^^y' {7T{k)N{p 'm^y/2)^+i^{y/2)\N{p)) 
Zn 



1 r 
J 



H^^ix, t, a, U.m^N) Nip h+j^Nip) 



^ Previously TV — > ttTV GPDs were briefly discussed in Ref . [ Pol9^ 



32 



+ Hi^'^ix, t, a, m^jv) N{p')i5N{p)\ . (97) 

Here « 92.4 MeV is the pion decay constant and we introduced the foUowing kinematical 
variables: 

P ^ i-^ , (98) 

and the momentum transfer 

A = p' + k~p. (99) 

n is a Ught-cone vector, i.e. = 0, which is normahzed such that n- P = 1. The N ttN 
GPDs H^^ are the functions of the invariants on which the matrix elements depend. 
These are the skewedness: 

e--^, (100) 

the square of the momentum transfer t = A^, the square of the nucleon momentum 
difference t^^ ^ [p ' — p)^, the invariant mass of the final TriV system m^^v and the light- 
cone fraction a, which characterizes the longitudinal fraction of the final pion: 

n ■ k , , 

a=^. (101) 

These GPDs can be used to describe the hard "semi-exclusive" processes of the type: 



J* + N (ttN) , J* + N ^ M+ (ttN) 



(102) 



in which the target nucleon dissociates into a low mass ttN system. The description of 
such reactions is important because in the high energy experiments the resolution in the 
missing mass of the recoil baryonic system can be rather low. In such cases the data for 
e.g. 7*iV always include the 7*7V j{TrN) 'contamination' due to a soft pion. 

Another motivation for studies of — > ttN GPDs is that in the threshold region, 
where m^^r — > m-jj+miq, a — > TOjr/(TOjr+m jv), etc, the i:N GPDs can be expressed in terms 
of nucleon and pion GPDs [ MosOl , StrOlj. Additionally one can develop the systematic 
chiral perturbation theory for the N — » ttN GPDs. Indeed the first Mellin moment of 
the A^ — > ttN GPDs ( p7| ) corresponds to the amplitude of the pion electroproduction off 
the nucleon. Near the thre shold this process has been described by chiral perturbation 
theory, for the review see [Ber94|. The second Mellin moment of the A^ — + nN GPDs 



corresponds to pion production by the operator of the energy momentum tensor, this can 
be rephrased as the pion production by a graviton. This simple example shows that the 
studies of the hard reactions of the type (102) open a new, large field for applications of 
chiral perturbation theory. Measurements of these GPDs would give a new sensitive tool 
for studies of chiral symmetry breaking mechanism in QCD. 

The A^ — > ttN GPDs in the resonance region, i.e. for m-^^ ~ TnR, give us a new tool 
for studying of the resonance spectroscopy, i.e. the spectroscopy with light-cone probes. 
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3 GPDs IN THE LARGE iV, LIMIT 



The large Nc limit in QCD |'tHo74, Wit79[ proved to be an useful tool to analyze the 



properties of mesons and baryons. The idea is that our world with the number of colors 
Nc = 3 is not qualitatively different from an imaginary world with large number of colors. 
For the c omprehe nsive review of applications of the large Nc expansion in hadronic physics 
see Ref. Mano98 |. In this section we study the GPDs of the nuclcon in the large Nc limit, 
first discussing general properties of the GPDs in this limit and then presenting calculations 
of the GPDs in the specific dynamical realization of the low-energy large Nc QCD - chiral 
quark-soliton model. 



3.1 Generalized parton distributions at large 
erties 



general prop- 



In order to consider the generalized parton distributions in the large Nc limit we start from 
the quark correlator {B2\4'{x2)'ip{xi)\Bi) . The Hartree picture of the larg e Nc n uclcon 
described in | Wit 79 leads to the following form of this correlation function | PobOC |: 



X / dRrB,iR)MR) 



Rf2fiPs2fi^sif[ ( a;? - 



2' 



(103) 



Here ^"s2/2,si/i (a^i — 2:2,xi,X2) is the correlation function {'ijj{x2)ip{xi)) corresponding to 
the static solution of the Hartree equation, Si are Dirac spinor indices and fi are SU (2) 
isospin indices. The x'^ — y'^ time dependence expresses the fact that we deal with a static 
solution. The solution -^32/2.81/1 (^^^ ^ xj, X2, Xi) also violates the SU{2) flavor and the 
3-space translation invariance. Therefore in Eq. (103) we used plane waves e*P "^ and the 
SU{2) r otator wave functions 4>b(R) expressed in terms of Wigner finite rotation functions 
[Adk8|: 



(104) 



in order to construct baryon states with given spin, isospin and momentum. The expression 
in the square brackets in Eq. ( |103| ) (where one puts (x — y)'^ =0) can be called a parton 
distribution of a soliton. 

Since the exact form of the effective Hartree hamiltonian corresponding to the large 
Nc limit generically is not known we do not know the function -^52/2,51/1 (2^2 ~ 2;",X2,Xi) 
eitheip]. However, following |Wit79| we assume that this solution has the spin-flavor 
symmetry: 

S,,,,{R)Rf,J,Fs,j,^,,j,(^X°^y°,Ospace{R)^,OspaU^^ (105) 



Fs2f2,sih{x -y ,x,y) 



(106) 



Here we use the notation Sg^gi (R) for spin rotations and Ospace{R) for space rotations, 
which are parametrized by a SU{2) matrix R. 

^®In the next section we shall use the chiral quark-soliton model to compute this function 
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Using the expression ( |103| ), one can obtain any quark parton distribution function 
in the large Nc hmit. Using the spin-flavor symmetry (106) of the soliton distribution 
(the expression in square brackets of Eq. ( |103| )) we can write explicitly the dependence on 
the orientation angles R of particular ■tfjip operators entering the definition of GPDs, see 
Eq. (||). In the leading order of 1/Nc we have: 



2tt 



7^ / dy-e' 



cP+i 



"(sol P'|V;/'(-2//2)7+^/(y/2)|sol P) 



y+=y±=o 

Ho{x,U)5ff' - , ^i(^'^'^) lTr(fltM(/-^')flr,-) • eg.feAfe , 

771 TV ^ 



(107) 



where R is SU{2) matrix which defines the orientation of the soliton in the spin-flavor 
space, M^f^ ^ is the matrix, of which all elements are zero except the element in row / and 
column /' which is unity. Furthermore in Eq. 107, _ffo,i(2;, ^, t) are the universal GPDs 
of the soliton. Analogously we can write for the quark helicity dependent GPDs of the 
soliton: 



— / d2/-e'-^^^"(solP'|Vi/'(-2//2)7+75V'/(y/2)|solP) 

ZTT 



A3 A* 



'■N 



'-TT:{R'^M^ff'^RTi) . 



(108) 



Now to get the GPDs for a baryon with particular quantum numbers (or for transition 
GPDs) we have to sandwich equations (107 108) between rotational states (104) of the 
soliton which corresponds to baryons from the octet and/or decuplet, i.e. 



(5' = r',5^,T^|...|5 = r,^3,T3) - j 



dR^^s'rf'iR) ... 4::f^{R). (109) 



We note here that in the derivation of Eqs. (107, 108) we used that in the large Nc limit 
the baryon masses rise linearly with Nc, with mass splitting m/\ — ttln ^ 1/Nc. Therefore 
in the first approximation we can neglect nucleon-A mass difference and consider them as 
nonrelativistic particles if ~t <C rn'j^. It is useful to work in the Breit frame, where: 

p' = -p = A/2^N°, 
1 



= ^(P'° + P°, 0, 0, 0) « (mAT, 0,0,0) ~ (7V„ 0,0,0) 



(110) 



we introduce a light-like vector n which is normalized by n • P = 1. In the large Nc limit, 
we can choose it in the form 



1^==— (1,0,0,-1) 



1 



(1,0,0,-1) 



which leads to 



t^N° and ^ = -i(n.A)«-^^l/Ar, 
2 ZiriN 



(111) 



(112) 



These relations fix the large Nc counting rules for baryon masses and kinematical variables 
in the problem. 
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Let us first apply the expressions for GPDs of the soliton (107 108) to the case of 
N ^ N GPDs. The projection of the rhs of Eqs. (107,108) onto particular nucleon states 
can be done with help of the standard integral: 



dR 



1 



1 



where C!j^2-j2m2 Clebsch-Gordan coefficient. Performing the projection we observe 

that in the leading order of the 1 /Nc expansion only the following flavor combinations of 
nucleon GPDs are nonzero: 



and 



i7«(x, e, t) + H\x, 1 1), H^ix, ^, t) - H^ix, ^, t) . 



E^{x, t) - E\x, 1 1), E^{x, t) - E^'ix, t) 



(114) 



(115) 



Gomplementary combinations are zero in the leading order of the l/Nc expansion. Further, 
taking into account the Nc counting rules for kinematical variables and for the nucleon 
form factors (see e.g. Ref. [Das95, Mano98| ) we obtain the following counting rules for the 
nucleon GPDs: 



and 



H"ix,tt)+H''ix,^,t) 
H^ix,tt)-H''ix,tt) 
H^{x,i,t)-H\x,i,t) 
H^{x,C,t)+H'^{x,C,t) 



E^x,^,t)-E''{x,tt) 
E^{x,^,t)-E\x,^,t) 
E^ix,^,t)+E''ix,tt) 
E^x,^,t)+E^{x,^,t) 



N^f{NcX,NcC,t), 
' Ncf{NcX,Nctt), 



N!f{NcX,Nc^,t). 
Ntf{NcX,Nci,t) 
N^fiNcX,Nctt) 
N^fiNcX,Nc^,t) 



(116) 



(117) 



In the above equations /(a;i,a;2,i) is a generic function stable in the large Nc limit. Also 
note that for the "truly non-forward" functions E, E the overall powers of Nc are higher 
than for the GPDs H,H because of the kinematical factor 1/mjv ^ and ^/mjv ~ 

l/Nc in the definition of the GPDs E and E respectively (see Eq. (||)). 

Taking the forward limit ^ 0,t ^ of the large- iVc counting rules (116) we 



reproduce the corresponding c ounting rules for polarized and unpolarized nucleon parton 
distributions derived in Refs. [ Dia96 , Dia97 |. 

Let us analyse the large Nc behavio ur o f particular contributions to GPDs. We start 
with the D-term contribut ion, see Sec. 2.2, The large- A'^c counting for the D-tcrm can 
be extracted from either ( |ll6[) or from (117), because the D-term contributes to both 
functions H and E in the same way (up to the sign). For the isosinglet component of the 
D-term we obtain: 



D^{z) + D'^iz) ^ Nlf{z) . 



(118) 
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When we turn to the isvector D-term we see that that the counting rules for iJ" — H"^ and 
i?" — E'^ are different the latter being the dominant one in the large Nc limit. It implies 
that for the counting of the large Nc behaviour of the isovector D-term D"(z) — D'^{z) one 
has to use Eq. (116). One obtains: 



D'-{z)-D''{z)^NJ{z). 



(119) 



Comparing Eqs. (|ll§| ) and (|119| ), one sees that the isovector D-term is suppressed relative 
to the isoscalar one in the large Nc limit. 

Next let us check that the p ion pole contribution of Eq. (^) to i?" — E"^ corresponds 
to the general counti ng ru le ( |117[ ). Indeed, taking into account that ~ Nc and rrif^ ^ Nc 
we easily repro duce (117) from Eq. (|4^). 

In section 3.2 we shall discuss the results obtained in the chiral quark-soliton model 
for the "large" GPDs (114,115). Before doing this we derive the large Nc relations between 
iV ^ and iV A GPDs introduced in Sec. O. 



3.1.1 GPDs for octet to decuplet transition in the large Nc limit 

In order to derive the large Nc relations between the iV ^ A GPDs and usual ones we 
shall use the fact that in the large Nc limit the nucleon an d A stat es are simply different 
rotational excitations of the same classical obje ct - solit on [ Fra98a |. 

To obtain the relations we sandwich Eqs. (107, lOS) between rotational states of the 
soliton which corresponds to baryons from octet and decuplet. Doing this we obtain the 
expressions for N ^ N and N A GPDs in terms of the universal functions i?o.i and 
Hq^i, i.e. the GPDs of the sohton introduced by Eqs. ( p^JlO^ ). Since both N ^ N and 
— > A GPDs are expressed in terms of the same GPDs of a soliton Hq i and iJo,i we 
obtain the following relations [The N ^ A GPDs are defined by Eqs. (^ji 



V3 



[E^{x,i,t)~E''{x,^,t)] 



Ci(x,^,i) = x/3[ff"(x,e,i)-i/''(.T,e,i) 

V3 



E\x,i,t)~E'{x,^,t) 



(120) 



All other A'' ^ A GPDs are zero at lead ing order of \/Nc expansion. The properties of 
these relations were discussed in Sec. 



3.2 GPDs in the chiral-quark soliton model 

Model calculations of the nucleon GPDs are important firstly to understand how non- 
perturbative mechanisms generate various structures in GPDs. In this respect the chiral 



models should be very useful because as we saw in Sec. 2.4 the spontaneously broken chiral 



symmetry plays a prominent role in the physics of GPDs. Secondly, model calculations 
allow us to develop an intuition about the structure of the GPDs which is important for 
the construction of phenomenological parameterizations for GPDs. 
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Model calculations of nucleon GPDs were first performed in Ref. | Ji97a| where GPDs 
were computed in the MIT bag model. Unfortunately in the bag model the chiral symmetry 
is not respected and also the non field theoretic nature of this model leads to problems 
with positivity of antiquark distributions. Calculations of the nucleon GPDs in the chiral 
quark-soliton model |Dia8S| do not show these deficiencies and have been performed in 
Refs. |Pct98, PcnOOa . The chiral quark-soliton model is field theoretic model of the nucleon 
which stresses the dominant role of spontaneous chiral symmetry breaking in the dynamics 
of the nucleon bound state. The review of foundations of this model and its applications 
to calculations of the parton distributions can be found in [ Dia98 , DiaOC |. The review of 
applications to calculations of nucleon form factors and many static observables can be 
found in Refs. Chr96| , Alk96 |. We refer the reader to above reviews for details. Here we 
remind only some basics of this model. 

The chiral quark-soliton model is essentially based on the 1 /Nc expansion. Although 
in reality the number of colors Nc — the academic limit of large Nc is known to be 
a useful guideline. At la rge Nc the nuc leon is heavy and can be viewed as a classical 
soliton of the pion field [Wit?!;; Adk83|. In the previous section we derived from this 
large- iVc semiclassical picture of the nucleon general large Nc counting rules for various 
nucleon GPDs (see Eqs. (116,117)) and relations between N ^ N and N A GPDs (see 
Eq. ( |l2Cl| )). In order to compute the nucleon GPDs we have to specify the dynamics of 
the low energy QCD. Here we work with the effective chiral action given by the functional 
integral over quarks in the background pion field |Dia83, Dia86|: 



exp(iS'eft[7r(a;)]) = / DijjDi; exp i / d^xip{if)- MU'^'-)^; 



U = exp[i7r''(x)T''] 



= exp[i7r''(x)T''75] = 



1 + 75 



U 



1 



75 



J/t. (121) 



2 2 

Here ip is the quark field, M the effective quark mass, which is due to the spontaneous 
breakdown of chiral symmetry (generally speaking, it is momentum dependent), and U is 
the SU{2) chiral pion field. The effective chiral action given by Eq. (121) is known to con- 
tain automatically the Wess-Zumino term and the four-derivative Gasser-Leutwyler terms, 
with correct co efficien ts. Equation (121) has been derived from the instanton model of the 
QCD vacuum [Dia86|, which provides a natural mechanism of chiral symmetry breaking 



and enables one to express the dynamical mass M and the ultraviolet cutoff intrinsic in 
(121) through the Aqcd parameter. The ultraviolet regularization of the effective theory 
is provided by the specific momentum dependence of the mass, M(p^), which drops to zero 
for momenta of order of the inverse instanton size in the instanton vacuum, 1/p 600 MeV. 
We also note the the inverse instanton size sets the renormalization scale for computed 
parton distributions at /i = 0.6 GeV. 



An immediate application of the effective chiral theory (121) is the chiral quark- 
According to this model nucleons can be viewed 



soliton model of baryons of Ref. |Dial 



as Nc "valence" quarks bound by a self-consistent pion field (the "soliton" ) whose energy 
coincides with the aggregate energy of the quarks of the negative-energy Dirac continuum. 
Similarly to the Skyrme model the large Nc limit is needed as a parameter to justify the 
use of the mean-field approximation; however, the 1/A^c-corrections can be — and, in some 
cases, have been — computed ]Chr9(: | . 
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Let us remind the reader how the nucleon is described in the effective low-energy 
theory Eq. (121). Integrating out the quarks in (121) one finds the effective chiral action, 



^offK(a;)] = -iV,Splogi^(C/), D{U) = ido-H{U), (122) 

where H{U) is the one-particle Dirac hamiltonian, 

H{U) ^ -ij^j'^dk + Mj^U'" , (123) 

and Sp . . . denotes the functional trace. For a given time-independent pion field U = 
exp(i7r''(x)T°) one can determine the spectrum of the Dirac hamiltonian, 

i/$„ - £;„$„. (124) 

It contains the upper and lower Dirac continua (distorted by the presence of the external 
pion field), and, in principle, also discrete bound-state level(s), if the pion field is strong 
enough. If the pion field has unity winding number, there is exactly one bound-state level 
which travels all the way from the upper to the lower Dirac continuum as one increases 



the spatial size of the pion field from zero to infinity |Dia88|. We denote the energy of the 
discrete level as i?iev, —M< Ei^^ < M. One has to occupy this level to get a non-zero 
baryon number state. Since the pion field is color blind, one can put quarks on that 
level in the antisymmetric state in color. 

The limit of large Nc allows us to use the mean- field approximation to find the nucleon 
mass. To get the nucleon mass one has to add NcEiev and the energy of the pion field. 



Since the effective chiral lagrangianis given by the determinant (122) the energy of the pion 
field coincides exactly with the aggregate energy of the lower Dirac continuum, the free 
continuum subtracted. The self-consistent pion field is thus found from the minimization 



of the functional [Dial 



niN = min J Ei.^U] + ^ {E,,[U] - E^^^) I . (125) 

I -E„<0 J 

From symmetry considerations one looks for the minimum in a hedgehog ansatz: 

C/c(x) = exp[i7r''(x)T''] = exp [m''T"F(r)] , r = |x|, n = -, (126) 



r 



where P{r) is called the profile of the soliton. 

The minimum of the energy ( |125| ) is degenerate with respect to translations of the 
soliton in space and to rotations of the soliton field in ordinary and isospin space. For 
the hedgehog field ( |126| ) the two rotations are equivalent. The projection on a nucleon 
state with given spin (S^) and isospin (T3) components is obtained by integrating over all 
spin-isospin rotations, R, see Eq. ( |109D . Analogously, the projection on a nucleon state 
with given momentum P is obtained by integrating over all shifts, X, of the soliton, 

(P'|...|P) = y"d3Xe'(P'-P)-^ ... (127) 

Recently it was shown that the chiral quark-soliton model of the nucleon provides a 



framework for a successful calculation not only of nucleon form factors [Chr96| but also 
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of the nucleoli parton distributions, both unpolarized and polarized, at a low normaliza- 
tion point pia96|, pia97|, |Pob96|, |Wei97|, pam98|, |Wak9^, |Wak99|, |Pob99|. Because both 



cases are limiting cases of the GPDs, the chiral quark-soliton model is a promising tool 
to calculate them. Since the chiral quark-soliton model is a quantum field-theoretical 
description of the nucleon, with explicit quark degrees of freedom, it allows an unambigu- 
ous identification of the quark as well as antiquark distributions in the nucleon. We were 
able to demonstrate that all general pro perties of the q uark and antiquark distributions 
(positivity, sum rules, Soffer inequalities | Pob96 FobOO , etc.). 



All GPDs can be expressed in terms of a matrix element of a non-local quark bi- 
linear operator in the nucleon states with definite 4-m omcnta P, P' and spin and isospin 
components, see general large Nc expression Eq. (103). Since in the chiral quark-soliton 
model the dynamics of the mean-field is specified, one can write the explicit expression for 
correlator of the type of Eq. (|l03| ). [We write explicitly all the flavor (/, g — 1,2) and the 
Dirac (i, j = 1, ...,4) indices for clarity]: 

(P', 5^, T^l ^,f,{x^, x) r' (2/°, y) |P, ^3, T3) = 
2mNNcJ d^Xe'^^'-'^"^-^ J dR <j)l,^r^,{R)<j)s,TAR) 



nj'i 



,,(x-X) <i>f:-'(y-X)(i?t)/' 



(128) 



The functions $„ are eigenstates of energy of the Dirac hamiltonian (123) in the 
external (self-consistent) pion field t/c, they can be found numerically by diagonalization 
of the hamiltonian. In this way we compute the function F in general large- A'^c Eq. (103). 



3.2.1 Results for H and E 

Having the expression ( |128| ) for nucleon matrix elements of quark bilocal operators in the 
chiral quark-soliton model we can now express various nucleon GPDs as the sum over 
occupied quark orbitals in the external pion field Uc- We start with GPDs H and E. One 
immediately observes that, in accordance with general large- iVc counting rules ( 
in the leading Nc order only the flavor singlet part of H{x,£^,t) and the flavor-nonsinglet 
part of E{x,^,t) are non-zero. They are given, respectively, by 

Y^m{x,t^') = ^ I dz' f d'X exp[zA.X] 

X ^ exp{zzO[(a; + e)™JV-S„]}*UX)(l+7V)«'n(X-z*'e3), (129) 

occup. 

E^^{x,^,t)^E'^{x,^,t)=-'^^^ J dz^J d'X exp[zA.X] 

X exp{^z°[(a;-f 0™JV--B„]}$t(X)T*=(l-f7V)*n(X-z°e3). (130) 

occup. 

The isovector part of i?(x,^,A^) and the isosinglet part of E{x,£,, A'^) appear only in 
the next-to-leading order of the 1/A^c-expansion, i.e., after taking into account the finite 
angular velocity of the soliton rotation. 
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First we would like to demonstrate that the two limiting cases of the GPDs — usual 
parton distributions and elastic form factors — are correctly reproduced within the chiral 
quark-soliton model. Taking in Eq.(129) the forward limit, A — > 0, one recovers the 



formula for the usual singlet (anti-) quark distributions in the chiral quark-soliton model 
which was obtained in Ref. [ Dia96[ . Thus the forward limit, Eq. (^, is reproduced. On 
the other hand, integrating Eqs. (129 130) over — 1 < 2: < 1 one obtains (up to corrections 
parametrically small in \/Nc) the expressions for the electromagnetic formfactors of the 
nucleon derived in Ref. 



i Diag 



I'^dx Y,H''{x,^,t)^N,J d'Xexp[tA-X] ^ $t(X)$„(X) 



(T=0) 



it), (131) 



occup. 



E 

occup. 



$t(X)7"7V<fn(X)=^^f-'^(t) 



(132) 



The electromagnetic formfactors computed in the chiral quark soliton model on the basis of 
these formulae compare rather well with the experimentally measured ones up to momenta 
of order ^ 1 GeV^ [ phr96[ |. One can also check that |Pct9S| , |PenOOa| the GPDs given 
by Eqs. (129, 130) satisfy the polynomiality condition (Oh. 




''min 



Figure 4: The GPD H''{x,^,t) + H'^{x,^,t) as a function of a; at ^ = 0.3 and t = t„ 
-0.35 GeV^ Dashed curve: contribution from valence level. Dashed-dotted curve:contribution of 
the Dirac continuum. Solid curve: the total distribution (sum of the dashed and dashed-dotted 
curves). The vertical lines mark the crossover points x = 



Eqs. (12£ ,13C ) express the GPDs as a sum over quark single-particle levels in the soli- 
ton field. This sum runs over all occupied levels, including both the discrete bound-state 
level and the negative polarized Dirac continuum. We remind the reader that in the case of 
usual parton distributions it was demonstrated that in order to ensure the positivity of the 
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antiquark distributions it is essential to take into account the contributions of all occupied 
levels of the Dirac Hamiltonian | Dia9(; ] . The so-called "valence level approximation" ad- 



vocated in Refs. |Wei97, Gam98| for structure functions in the chiral quark-soliton model 



leads to unacceptable negative antiquark distributions. We shall see below that also in the 
off-forward case the contribution of the Dirac continuum drastically changes the shape of 
the distribution function, leading to characteristic crossovers of H{x,£,,A) at \x\ = ^, in 
addition this contribution generates the large and negative D-term. 

Now we present results for calculations of the contributions of the discrete bound- 
state level and the negative Dirac continuum to Eq s. ( 12S|.13C| ). Technical details and 
methods of calculations can be found in Refs. |Pct98, PenOOa . 
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Figure 5: Ratio Eq. (134) as a function of x showing the deviation from the factorization ansatz 
for the t-dependence of the GPD H for the values of ^ = (filled circles) and ^ = 0.2 (open 
circles) with to = -4m%£,'^/{l - ,f^) and ti = to - 0.5 GeV^. 

In Fig. H a typical shape of the GPD H^{x,^,t) + i?''(x,^,t) as a function of the 



variable x is shown at ^ = 0.3 and t 



tmin, with 



AmU'l{l-e)- 



For ^ = 0.3 one has tmin —0.35 GeV^. On this figure we show separately the contribu- 
tions of the valence level and of the Dirac continuum as well as their sum. We see that 
the Dirac continuum contribution is essential especially in the region — ^ < a: < also in 
the region x < — ^, corresponding to minus antiquark distributions, the Dirac continuum 
contribution ensures the positivity of antiquarks. We also computed several McUin mo- 
ments at various values of ^ and t (in the model calculations the values of ^ is limited by 
< —t/ (4m^ ~ i)) and with the help of Eq. ( |2l] ) we determined the value of the first few 
coefficients in a Gegenbauer expansion of the D-tcrm (|2^). The result extrapolated to the 
value of t = is: 



q—u,d,s.. 



A.QCl^^{z)- 1.2 Cl'\z) 



3/2, 



0.4 Cl'\z) 



(133) 
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We shall use this result for the D-term in Sec. |j for estimates of various observables. 

Calculations of GPDs in the chiral quark-soliton model give a possibility to study 
their dependence. In modeling the ^-dependence of GPDs usually the factorization 
ansatz: H{x,^,t) — H{x, £^)Fi{t) is adopted. In order to illustrate numerically the non- 
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Figure 6: Ratio Eq. (134) as a function of x for a fixed value of ^ = and for to = and for 
ti = —0.25 (open circles), —0.5 (filled circles), —0.75 (triangles) GeV^. 

factorazibility of the t-dependence of H{x,£_,t) we show in Fig. || the following ratio : 



R = 



H{x,£„ti) J^^dx H{x,^,to) ' 



(134) 



[H — + H'^] as a function of x taking to = t^i^ and ti — t^in — 0.5 GeV'^, and ^ = 
(filled circles), ^ = 0.2 (open circles). If the i-dependence of the GPDs fact ors out then 
R = 1 independently of to, ti and x, From Fig. ^ we see that the ratio (134) can deviate 
considerably from unity and it depends in a non-trivial way on both variables x and ^. 

To illustrate the dependence of the ratio ( |134| ) on the value of ti we plot in Fig ^ this 
ratio as the function of a; with <^ — Oandio = for several values of ti = —0.25, 0.5, 0.75 GeV^. 
The results of the chiral quark-soliton model on the t— dependence of GPDs shown in Fig. ^ 
and Fig. ^ illustrate that the assumption of factorization of the dependence for GPDs 
should be taken with caution and that more work on t— dependence of phenomenological 
parametrizations of GPDs should be done. 

It is interesting to note that the features of the t-dependence of GPDs shown on Fig. ^ 
can be qualitatively reproduced by the following ansatz for H'^{x,£^ = 0,t): 



J?'(x,e = 0,t) = 477 g(^), (135) 
x" ^ 

where a' can be interpreted as the slope of a Regge trajectory. This Regge theory motivated 
ansatz reproduces all qualitative feature of Fig. || obtained in the chiral quark-so liton 
model. As a curiosity we note that a rough fit of the curves on Fig. || by the ansatz (135) 
gives a' « 0.8 GeV~^, amazingly close to the slopes of meson Regge trajectories. 
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Let us now discuss results for the GPD E'^ — E"^ . In Fig. |^ we plot the GPD 
E'^{x,^,t) — E'^{x,£,,t) at ^ = 0.3 and t — ^min- We see that in this case the Dirac con- 
tinuum contribution is also essential and it generates a symmetric sea contribution to the 
GPD E^ — E''-. Computing several Mellin moments of the function E" {x, ^, t) — E'^{x, ^, t) 
we are able to extract the value of the isovector D-term D'^{z) — D'^{z). The leading- TVc 
result for this flavor combination is zero, i.e. 

D''{z)-D'^{z) = 0. 



This is in accordance with general large Nc counting rules for the D-term (119). Therefore 
we conclude that the D-term is independent of quark flavor. 




Figure 7: The GPD E'^{x,^,t)-E'^{x,^,t) as a function of s at ^ = 0.3 and f = Wn = 0.35 GeV^ 
Dashed curve: contribution from valence level. Dashed-dotted curvexontribution of the Dirac 
continuum. Solid curve: the total distribution (sum of the dashed and dashed-dotted curves). 
The vertical lines mark the crossover points x — 

For the construction of the phenomenological model for the GPD E"^ we shall need 
the forward limit of this functions which we denote as e'^ {x) : 

e''{x)= hm^E'^{x,^,t). (136) 

In Fig. ^ we present the result of the model calculation of e"{x) — e'^''{x). We see that it 
consists of the valence contribution and of the sea part peaked around x = 0. With help 
of sum rules for the function e"(a;) — e'^{x): 

J dx [e"(x) - e''(x)] = kP-k", 

j dxx [e"(a;) - e\x)\ = 2 [ J" - J'^] - + M!^ (137) 

where is an anomalous magnetic moment of a nucleon N ^ p,n and J"? and M| are 
the parts of angular and of total momentum carried by a quark, we obtain in the model: 

- k" = 1.0™, -f 2.22,ea = 3.22 vs. 3.7 (exp.), 
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Figure 8: The GPD £"(a;,^ = 0,t = 0) ~ E''{x,^ = 0,t = 0) = e"(x) - e'^{x) as a function of x. 
Dashed curve: contribution from valence level. Solid curve: contribution of the Dirac continuum. 

2 [J" - J''] - + = 0.18™; + Osea = 0.18 . (138) 

Here we show also the individual contributions from the valence and sea parts. We see that 
the value of the anomalous magnetic moment is dominated by the sea contribution and 
its value is close to the experimental number. The model prediction for 2 J"^'' — 2A/2 
is dominated by the valence contribution. 

Since the combination J" — J"^ of angular momenta is dominated by the valence 
part of the function e"(a;) — e''(a;)0 it is instructive to compare the shapes of the valence 
contributions to e'^{x) and to q{x) with each other. The corresponding comparison is shown 
in Fig. 1^ where we plot 



- eialix) Qvalix) + italix) 



versus 



/-I dx {e^^iix) - etiix)) j\ dx (^^(x) + qti{x)) ' 

computed in the chiral quark-soliton model. We see that the shapes of both these valence 



distributions in the model are very close. This observation will motivate us in Sec. 4.2.3 
for the ansatz: 

el^i{x)=A^q,ai{x), (139) 

where A'^ is a ^-independent constant. 

The qualitative features of the GPDs H and E obtained in the chiral quark-soliton 



model will be used in Sec. 4.2 for the construction of phenomenological parametrization 
of these distributions. 

3.2.2 Results for H and E 

Now we turn to the results for the quark helicity dependent GPDs H and E. Expressions 
for these distributions in terms of a summation over occupied quark orbitals in the external 



pion mean field can be found in Ref. [ PenOOa |. Here we discuss only the main results 



^'^As shown in Ref. [ |Pob99| the AfJ — A/j is also dominated by valence level contributions. 
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Figure 9: Comparison of shapes of the valence level contributions to u{x) + d{x) (solid curve) 
and to e^{x) — e'^{x) (dashed curve). Both contributions are normalized to unity. 



According to the general large Nc counting rules derived in Sec. 3.1 at the leading 



order only isovector combinations of GPDs _ff " — H'^ and i?" — E'^ are non-zero. This 
general counting is fulfilled in the chiral quark-soliton model. 

Before presenting the numerical results for the GPDs H and E let us discuss specific 
contribution to these GPDs originating from the long range pion tail of the pion mean- 
field. The behaviour of the mean pion field at large distances is governed by linearized 
equations of motion and PCAC: 



lim U{x) = 1 + — — -^(1 -hm^|f|) — p— exp(-TO^|x|) , 



(140) 



where qa ~ 1.267 is the axial charge of the nucleon, « 92.4 MeV is the pion decay 
constant. Con tribution of this long-range tail to E can be computed analytically (see 
details in Ref. |PenOOa|) with the result: 



Ml/2 



mNF{t) f d'^p 



P. 
(p + A)2 



(2^)4 



{P - A)^ 



{p + A)2 _ M2 ?0 {p- A)2 - M2 iO 



+ u 



(141) 



Here M[p'^] is the momentum dependent constituent quark mass as obtained in the theory 
of th e instanton vacuum | Dia86t| . It provides the UV regularization of the effective theory 
(121). In the Eq. (141) we introduced the following form factor: 



Fi-k') 



3fc3 



d'^x exp(ifc ■x)Tt U{x) - I t' 



(142) 
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The small momentum asymptotic of this form factor can be obtained by substituting into 



Eq. (142) the long-range asymptotic of the mean-field (|140|). The result is: 



lim Fit) . 1^ 
t^ml ml — t 



(143) 



Now the crucial observation [PenOOaj is that the integral overp in Eq. (141) coincides 
exactly (up to trivial renaming of variable) with the expression for the light-cone pion 
distribution amplitude in the instanton model of the QCD vacuum [ Pet99 . Therefore the 



expression (141) for the at t can be written in the compact form: 



This is exactly the general result obtained in Sec. 2.4.2, see Eq. m 



(144) 
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Figure 10: Comparison of pion pole contribution and non-pole part of the isovector GPD E at 
various values of t and ^. The positive curves correspond to pion pole contributions. 



In Ref. [PenOOa| the non-pole contributions to i?" — E''' have been computed in the 
chiral quark-soliton model. The corresponding results are presented in Fig. |l^. We see 
that for wide range of t and ^ the pion pole contribution dominates, the non-pole part 
being negative. The result of the chiral quark-soliton model for iJ"(a;, ^, t) — H'^(x, t) at 
^ = 0.3 and t — tmin is presented in Fig. |ri|. We observe the sizeable contribution of the 
Dirac continuum. 
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-0.2 



Figure 11: The GPD - H'^{x,^,t) as a function of a; at ^ = 0.3 and t = f^in = 

-0.35 GeV^ Dashed curve: contribution from valence level. Dashed-dotted curve:contribution of 
the Dirac continuum. Solid curve: the total distribution (sum of the dashed and dashed-dotted 
curves). The vertical lines mark the crossover points x = 



To illustrate the non-factorizability of the i— dependence also for the GPD H we 
plot in Fig. [ij the ratio analogous to Eq. (134) with the same values of all parameters 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

Figure 12: The same as Fig. | but for the GPD 5"(a;,C,t) - H''{x,^,t). 

replacing GPD i?"(a;,C,t) + H'^{x,£,,t) by H'^{x,^,t) - H'^{x,£^,t). We observe again a 
noticeable deviation from the factorization ansatz for the dependence, though in this 
case the deviation is smaller. 

With these lessons from chiral quark-soliton model calculations we next turn to a 
discussion of the phenomenological applications of the GPDs to various hard exclusive 
reactions. 
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4 GPDs FROM HARD ELECTROPRODUCTION RE- 
ACTIONS 



4.1 Introduction 

In this section, we show how one can access the GPDs through the measurement of hard 
exclusive reactions. As discussed in Sec. in leading twist there are four GPDs for the 
nucleon, i.e. H, E, H and E, which are defined for each quark flavor (u, d, s). These 
GPDs depend upon the different longitudinal momentum fractions a; + ^ (a; — ^) of the 
initial (final) quark and upon the overall momentum transfer t to the nucleon. Therefore, 
these functions contain a wealth of new nucleon structure information. 

In the following, it will be shown how the different GPDs enter in the observables 
for different hard exclusive processes. We will discuss in detail the deeply virtual Comp- 
ton scattering (DVCS) process on the nucleon and hard electroproduction of vector and 
pseudoscalar mesons (HMP) on the nucleon. Besides the hard processes involving nucleon 
GPDs, it will also be discussed how the iV — + A DVCS process and the hard electro- 
production of ttA final states depend upon GPDs for the N ^ A transition, providing 
new information on the quark structure of the A resonance. Furthermore, it will be dis- 
cussed how the electroproduction of kaon-hyperon {KY) final states contains GPDs for 
the N —> Y transitions, yielding information on the quark structure of hyperons. 

It is clear that to extract all those different GPDs from hard electroproduction reac- 
tions and to map out their dependencies on the three variables (x, ^, t) is a quite challeng- 
ing program. As a first step towards this goal, one strategy is to use phenomenological 
parametrizations for the GPDs which include all general constraints and which are flexible 
enough to allow for a fit to different hard electroproduction observables. The aim of this 
section is to present such phenomenological and physically motivated parametrizations of 
the GPDs and to test then their sensitivity on observables. 



4.2 Phenomenological parametrization of GPDs 

In the following, we discuss successively the parametrizations for the GPDs H"^, E'^, H'' 
and E'^ for each quark flavor q (= u,d,s). As wc address here mainly hard exclusive 
reactions in the valence region at relatively small values of —t { < 1 GeV^ ), we start by 
parametrizing the GPDs at t = and discuss firstly their x and ^ dependencies. For each 
GPD, we then come back to their t-dependence. 



4.2.1 Parametrization of the GPD 

For the function H'' (for each fiavor q), the i-independent part H'^{x,^) = = 0) 

is parametrized by a two-component form : 



1 



H^ix,0 = H%j,{x,^) + 6{i-\x\) — D{-), 



(145) 



where H'^jj is the part of the GPD which is obtained as a one-dimensional section of a 
two- variable double distribution (DD) F'^ |Rad99|, imposing a particular dependence on 
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the skewedness ^, as given by Eq. (IT] 



dl3 



1-1/31 



1+1/31 



da (5(x-/3-aO F''{(3,a) . 



(146) 



The D-tcrm contribution D in Eq. ( |145| ) completes the parametrization of GPDs, restoring 
the correct polynomiaHty properties of GPDs |Ji98b, Pol99t| as discussed in Sec. 2.2. It 
has a support only for \x\ < |^|, so that it is 'invisible' in the forward limit. Because the 



D-term is an isoscalar contribution, it adds the same function for each flavor. In Eq. (145), 
Nf = 3 is the number of active flavors (u, d, s). 

For the double distributions, entering Eq. (|146D, we follow Radyushkin's suggestion 



|Rad99| to use the following model 



(147) 



where q{(3) is the forward quark distribution (for the flavor q) and where /i(/3, a) denotes 
a profile function. In the following estim ates, w e parame trize the profile function through 
a one-parameter ansatz, following Refs. | Rad99| , Rad01t| : 



h{j3, a) 



r(26 + 2) [(1- 



22''+ir2(6+l) (1 



\2b+l 



(148) 



In Eq. ( |148| ), the parameter b characterizes the strength of the ^ dependence of the GPD 
H{x,^). The limiting case & ^ oo corresponds to the ^ independent ansatz for the GPD, 
i.e. H{x,£) = q{x), as used in Refs. pui98| , |Vdh98[ . The power b in Eq. (|l4|) is a free 
parameter for the valence contribution (byai) and for the sea/antiquark contribution (fo^eo) 
to the GPD, which can be used in such an approach as fit parameters in the extraction of 
GPDs from h ard ele ctroproduction observables. As an example, the twi st-2 DVCS predic- 
tions of Ref. I Vdh9£] and the twist-3 DVCS predictions of Ref. |Kiv01t] correspond to the 
choice bjjai = ^sea = 1-0. In the following, we discuss the dependence of DVCS observables 
on the shape of the profile function, i.e. on the choice of the parameters byai and bsea- 

In Eq. ( 147 ), we use the phenomenological forward quark distributions q{(3) as mea- 
sured from DIS as input, ensuring the correct forward limit for the GPDs, i.e. 7?'"(x, 0) = 
u{x), H'^{x,0) = d{x), H^{x,0) = s{x). Specifically, the ca lculation s shown in the follow- 
ing were performed using the NLO quark distributions of |Mar98a (more precisely, their 
so-called "central gluon parametrization" ) . We checked explicitely however that the results 
in the valence region practically do not change when using another current parametrization 
for the NLO quark distributions, as those parametrizations are very close to each other in 
the valence region which we consider here. As an illustration, we show in Figs. |l^, |l^ the 
parametrization of Eqs. (145 148) for the double distribution part to the u-quark GPD 
H"{x,^) and d-quark GPD _ff^(x, ^) respectively, for the parameter value byai — bsea — 1- 
One sees how the GPDs approach the forward u- and d-quark distributions as ^ — > 0. 
Note that the negative x-range at ^ = corresponds to the anti-quark distribution, i.e. 
q{-x) = -q{x). 

Fig. |l^ shows, at a fixed value of ^ = 0.3, the dependence of the DD part to the GPD 
on the parameters bvai,bsea- In the limit b^ai = bsea = oo, one finds back the foward 
w-quark distribution. 

We next discuss the parametrization of the ^-dependence of the double distribution 
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Figure 13: Double distribution part to the w-quark GPD ff" at t = for different ^, using the 
ansatz ( |f46| ) as desribed in the text for byai ~ hsea = 1- Upper panel : valence ii-quark GPD, 
lower panel : total (valence and sea) u-quark GPD. The thin lines = 0) show the ordinary 
forward ii-quark distributions (MRST98 parametrization &t — 2 GeV'^). 
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Figure 14: The same as Fig. ^ but for the double distribution part to the d-quark GPD H'' at 
t = 0. 
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Figure 15: Double distribution part to the GPD H^{x,(_,t = 0) at ^ = 0.3 for different values 
of the parameters h^ai, 6sea- The limit b^ai = &sea = oo corresponds to the forward ii-quark 
distributi on. The D-term contribution to , as it comes out from the chiral quark soliton model 
(see Eq. (133)) is shown by the dotted curve, and should be added to the double distribution part 
to obtain the total GPD H". 



part to the GPD H'^{x, ^, t). As for the hard exclusive reactions which we discuss t should 
be a small scale compared to the hard scale of the reaction, we will only be concerned 
here with the small —t<l GeV^ region. The main constraint which one has on the 
f-dependence of the GPDs comes through the sum rules of Eqs. (|^|lo|), through which the 
first moment of the GPDs is given by the elastic nucleon form factors. For the GPD H'^ 
and for a quark of flavor q, the first moment is determined through Eq.(^. in terms of 



the elastic Dirac form factors F^{t) for a quark of fiavor q in the nucleon. In Eq. (12), the 
strange form factor F^, which is not so well known and expected to be small, is set equal 
to zero in the following numerical evaluations. 

The simplest parametrization of the t dependence of the GPD H'^ in the small —t 
region which fulfills the sum rule of Eq. (||), consists of the factorized ansatz for the t- 
dependence : 

H^{x,tt) = ff"(x,0i^i"(t)/2, 

H'{x,i,t) - 0, (149) 

where F]"(t) and F'}{t) are determined through Eq. (|l^) using the empirical parametriza- 
tions for the proton and neutron Dirac form factors. Besides satisfying the sum rule 
constraint, the ansatz of Eq. ( |149D also yields the correct forward limit for u- and d-quark 
distributions, i.e. i?"(x,^ = 0, i = 0) = i?"(a;, ^ = 0) = u{x) and 7?''(x,^ = 0, i = 0) = 
H'^ix,^^ 0) = d{x). 

Note that the sum rule of Eq. (H) only gives a constraint for the valence quark part to 
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the GPD, as the sea-quark contribution drops out of this sum rule. Indeed at ^ = 0, where 
the GPDs reduce to the forward quark distributions, one sees from q{x) — qvai{x) + q(x) 
and <z(— = that the sea-quark and the anti-quark contributions cancel each other 

in the sum rule of Eq. (||). One could therefore use a different t-dependence for the sea- 
quark part to the GPD H'^. However, in lack of more information, we adopt here the same 
t-dependence for both valence and sea-quark parts. 



Such factorized forms for the t-dependence as in Eq. (149) have been used in prac- 
tically all theoretical estimates for hard exclusive reactions at small —t up to now. It 
remains to be investigated how realistic such parametrizations are if one goes away from 
the t = region. In particular, the i-dependence of the differential cross section for 
electroproduction on the proton, in which the GPDs iJ" and H '^ enter (see further), has 
been measured at HERMES at values of xb around 0.1 - 0.15 | Tyt01 |. These cross sec- 
tions clearly display an exponential fall-off with —t in the range Unin <-t< 0.4 GeV^. A 
form factor type ansatz as in Eq. (149), even if it gives the correct normalization at tmin, 
deviates over this t range already noticeably from an exponential fall-off. 

We also saw in Sec. |3.2.1 that the calculations of GPDs in the chiral quark-soliton 
model indicate a deviation from the factorization ansatz (149) for the dependence of 
GPDs, as seen from Figs. ^j^. The chiral quark-soliton model results for the i— dependence 
of GPDs in the small —t region and at a low normalization point can be described by a sim- 



ple Regge theory motivated ansatz (135). If one uses such a simple ansatz for H'^{x, ^ = 0, i) 
then the nucleon form factors follow from the sum rule Eq. (H): 



Flit) 



1 



(150) 



Here qvai{x) is the valence quark distribution at a low normalization point. The ansatz of 
Eq. ( |150[) is only valid for the small —t region, because at larger values of —t the integral 
in Eq. ( |150| ) is dominated by the large x region, for which a Regge ansatz does not hold. 
Furthermore, note that the model of Eq. ( |150D for the nucleon Dirac form factor as well 
as the ansatz of Eq. ( |135[ ) is valid only at a low normalization point. Strictly speaking the 
form of the ansatz of Eq. ( |135| ) should be changed when one goes to a higher normalization 
point. In particular, when increasing the normalization point, the slope a' should decrease 



[Bro94, Fra96|. This should lead for instance to the disappearence of the shrinkage of the 
diffractive cone in hard exclusive processes when increasing 

In F ig. ^ we show the ratio of the proton Dirac form factor Fl{t) obtained from 



Eq. (150), using the MRST98 parametrization for the valence quark distributions at the 



scale /X = 1.4 GeV (which will be relevant in the discussion of observables further on) to 
the dipole form of this form factor, given by: 



dipoio/.N _ 1 - (1 + kP) t/Am% 



= (t) 



1 - t/AT 



(1-^/0.71)2 



(151) 



The result is shown for several values of the slope a' = 0.75, 0.8 and 0.85 GeV^^, which are 
close to the phenomenological slope of Regge trajectories. It is seen from Fig. ^ that the 
simple model of Eq. ( |150D with e.g. a' = 0.8 GeV~^ gives a rather satisfactory description 
of the proton Dirac form factor over the range < — t < 1 GeV^. This finding therefore 



suggests that the Regge theory motivated ansatz of Eq. (135) can be used as a guide for 



53 



p DIPOLE 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

-t (GeV^) 



Figure 16: The ratio of t he p roton Dirac form factor obtained from Eq. ( |l50| ) to the dipole form 
of the Dirac form factor (151) for different values of the slope a' (in GeV~ ) as indicated in the 
figure. 



a more realistic parametrizations of the dependence of the GPDs. From the theoretical 
point of view it can give a clue for understanding the interplay between the Regge and 

partonic pictures of hard processes. 

The ^-dependence of the Regge type ansatz of Eq. (135) for H{x,£^ — 0,t) can be 
restored from Eq. (|l^) by using the following model for the double distribution : 



(152) 



where /i(/3, a) is a profile function as in Eq. (148). 

In Fig. |l^, we compare the factorized model of Eq. ( |149D for the 'envelope' function 
H{^,£,,t) (i.e. along the line x — ^) with the unfactorized model of Eq. (152). One sees 
that the unfactorized model leads to an increasingly larger reduction when going to smaller 
^ (typical for a Regge type ansatz) and to an enhancement at the larger ^. 

Since the theoretical status of the ansatze of the type (135,152) still remains to be 
clarified, we restr ict ourselves in the discussion of most observables, to the simple factoriza- 
tion ansatz ( |l49|) unless explicitely indicated. In some observables (like spin asymmetries 
for hard meson production) the dependence of the GPDs tends to cancel in the ratio of 
cross sections and therefore these observables are not very sensitive to the explicit form 
for the dependence of the GPDs. However, for absolute cross sections and also for the 
DVCS single spin asymmetry, the ansat z for the dependence of the GPDs can influence 
the prediction, as will be shown in Sec. 4.3.5. 



4.2.2 Parametrization of the D-term 

To c omp lete the parametrization of the GPD H, wc next specify the D-term part in 
Eq. (145). The D-term contributes only to the singlet GPDs H{x,£^), i.e. yields the same 
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Figure 17: Comparison of the ^-dependence of the 'envelope' function at two values 

of ~t (in GeV^). The function //"(C,C,0) (full curve) is calculated according to the double 
distribution ansatz of Eq. (146) with b^ai = fesea = 1, using the MRST98 parametrization at a 
scale fi^ — 1.4 GeV'^ for the forward -u-quark distri butio n as input. For the function H'"{^,£,,t) 
at —t = 0.55 GeV^, the factorized model of Eq. (149) (dotted curve), is compared with the 
unfactorized parametrization of Eq. (15i) (dashed curve). 



contribution for each quark flavor. 

As discussed in Sec. |2.2| , the function D{z) = D{z,t = 0) in Eq. (145) is an odd 
function of its argument, i.e. D{—z) — —D{z), and can be expanded in odd Gegenbauer 
polynomials as given by Eq. (p3|). The moments di, da, o?5 in such an expansion have been 
estimated in Eq. (133) using the calculation of GPDs in the chiral quark soliton model 
|Pet98| at a low normalization point fi « 0.6 GeV. In principle, the D-term also has a scale 
dependence but it will be neglected as for the moment the uncertainties in the modeling 
of the D-term are larger than the logarithmic scale dependence, for the scales relevant for 
present and planned experiments. 

The general expansion of the D-term in Eq. (23) can be used in phenomenological fits 
to data, and the moments di, da, ds, ... can be used as fit parameters to be extracted from 
hard electroproduction observables which are sensitive to the D-term. For the following 
estimates, we will use the values of Eq. (133) obtained from the chiral quark soliton model 
calculation. The r esulting D-term is represented in Fig. (|l5| ) for one quark flavor (i.e. by 
dividing Eq. ( [l33|) by N f — 3. Note the negative sign of the Gegenbauer coefficients for 
the D-term in Eq. (133), as obtained in the chiral quark-soliton model, which fixes the sign 
of the D-term relative to the double distribution part to the GPD H'^. 

In Fig. |l^, we display the full x and ^ dependence of the parametrization of Eq. ( |l45| ) 
for the GPD iJ"(x,^,t = 0), i.e. including both the DD part and the D-term. By moving 
along lines of constant ^, one nicely sees how the sensitivity of the GPDs to the D-term 
part increases relative to the double distribution part when increasing the skewedness 
parameter ^. In the end, at ^ = 1, the D-term contribution dominates the GPD H. 
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Figure 18: The x and ^ dependence of the GPD H'"{x,£,,t = 0) for the it-quark distribution, 
including the double distribution part (for the value b^ai ~ bsea ~ 1) and the D-term. 



The <-dependence of the D-term is again not constrained through the sum rule of 
Eq. (H), because the D-term drops out of this sum rule. In absence of more information, 
we will assume the same i-dependence for the D-term as for the DD contribution to iJ^. 
The study of the t-dependence of the I?-term is a subject for future work. 



4.2.3 Parametrization of the GPD E"^ and quark contribution to the proton 
spin 

The parametrization of -E', which corresponds to a nucleon helicity flip amplitude, is 
more difficult as we don't have the DIS constraint for the x-dependence in the forward 
limit. In the amplitudes for hard electroproduction processes at small —t, E'^ is multiplied 
by a momentum transfer and therefore its contribution is suppressed at small t in the 
observables. 

One contribution to the GPD E"^ is however determined through the polynomiality 



condition as explained in Sec. 2.2, which requires that the D-term contributes with opposite 



sign to H and E. This leads to the D-term contribution for : 

El-t.rJ^,U) - -e{(~\x\)l-D(^^,?j , (153) 

which guarantees that the D-term contribution is cancelled in the combination H + E. 
Therefore, the D-term part does not enter in the angular momentum sum rule of Eq. (p2[). 



Similarly to Eq. (145) for H"^, one can parametrize E'^ by adding a double distribution 
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part to the D-term of Eq. (153). This leads to the parametrization 



E''{x,^,t) = El^{x,^,t) ~ e{^~\x\)^D(^^,t] , (154) 



where E'j^j^ is the double distribution part. 

In the forward limit, the DD part reduces to the function e'^{x) = E'^jj{x, ^ = 0, t = 0) 
which is a priori unknown. However one constraint exists on the normalization of Ej^jj 
due to the form factor sum rule Eq. (^) in terms of the Pauli form factor F|(t) for the 
quark of flavor q in the proton. The quark form factors i^j' ^ind F2 are normalized at t = 
through the anomalous magnetic moments of proton and neutron, i.e. -F|'(0) — kP = 1.793 
and i^2"(0) = = -1.913. Defining k" ee F^[0) and k'* ee F^{0), one finds : 

k" = 2kP + k" = 1.673, 

k'^ = + 2 k" = -2.033. (155) 

Taking Eq. in the forward limit (^ = 0, t = 0) leads to the normalization condition 

-1 

dx e«(.T) = . (156) 

L 

As a first guess for the function e'?(a;), which is consistent with this normalization condition, 
one can assume that it has the same x dependence as a valence quark distribution, i.e. 

e"(x) = iu„Q;(a;)K", 
e''(a;) = dyai{x) n"^ , 

e'{x) = 0. (157) 

Indeed, it is seen from Fig. || that the chiral quark-soliton model calculation for the isovec- 
tor function e"(a:) — e'^{x) (which is leading in the large Nc limit) contains such a valence 
quark component. 



Using the ansatz of Eq. (157), one obtains a prediction for the total angular momen- 



tum of the quarks in the nucleon, evaluating Ji's sum rule at ^ = : 



1 



2 



Ji = - dx X [q{x) + e'(x)] . (158) 



-1 



By defining the total fraction of the proton momentum carried by the quarks and anti- 
quarks of flavor q as 



r-l /■! 

ri - 



M| = / dxx [q{x) + q{x)] ^ / dx x [q^ai{x) + 2q{x)] , (159) 
Jo JO 

and the momentum fraction carried by the valence quarks as 



M|"-' = / dxxq^aiix), (160) 
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the ansatz of Eq. (157) for e'^{x) leads to an estimate for the total angular momentum 
carried by the u, d and s-quarks in the proton, J", J'^, and J'^ respectively, as : 



J" = - I M. 

-- 

2 

r = - 



2 2 



(161) 
(162) 
(163) 





Mf = 1 GeV2) 


M| (/x2 = 1 GeV2) 


2 J9 = 1 GeV2) 


u 


0.34 


0.40 


0.69 


d 


0.14 


0.22 


-0.07 


s 





0.03 


0.03 


u + d + s 


0.49 


0.65 


0.65 



Table 4: Esti mate of 2 J' at the scale fj? = 1 GeV^ according to Eqs. (161-163) using the valence 
model of Eq. (157). For the forward par ton distr ibutions, the MRST98 parametrization (so-called 
'central gluon parametrization' of Ref. [VIar98a|) has been used. 



In Table ^, we show the values of the momentum fractions A/, 
scale /i^ = 1 GeV^, using the MRST98 parametrization [Mar98a| for the forward parton 



and M| at the 



distributions. We also show the estimate for J", J"^, and J'* of Eqs. (161-163) at the same 
scale. These estimates lead to a large fraction (69 %) of the total angular momentum of 
the proton carried by the u-quarks and a relatively small contribution carried by the d and 
s-quarks. Using this estimate, the sum of u, d, and s-quark contributions is around 65 % 
at this l ow scal e. This estimate is in good agreement with the recent lattice QCD estimate 
of Ref. UMatOOj for J?, which found the total v alue 2 J *? = 0.60 ±0.14 at a comparable low 
scale, and the QCD sum rule estimate of Ref. | Bal97 |. 

One furthermore sees from Table | that 2 J" -I- 2J'^ + to a rather high 

accuracy. For the sep arate u and d contributions however, one sees from Table ^ that in 
the ansatz of Eq. (157), 2 J" is quite different from M2', and similarly 2 J"* is quite different 



from M2 . This leads us to the observation that the following remarkable relation between 
nucleon anomalous magnetic moments and valence quark distributions in the proton seems 
to hold in general : 

dx X [u{x) — u{x) ] w —2 K / dx X {d[x) -d{x)\ . (164) 











Such a relation allows us to express the ratio of anomalous magnetic moments of proton 
and neutron in terms of the proton momentum fractions M|"°' carried by the valence 
quarks (see dchnition of Eq. (16C)) as follows: 



1 4M° 



2 m; 



(165) 



At this point one could say: 
exiguitas non caperet" 



..cuius rei demonstrationem mirabilem sane detexi. Hanc marginis 
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To check the vahdity of the above relation, wc show in Fig. n9 the scale dependence of 



-0.92 



-0.93 



-0.935 



-0.94 - 



-0.95 









CTEQ5M 






GRV98NLO 




k'/k" 






MRST98 







10 



10 

(GeV^) 



Figure 19: Scale depend ence of the rhs of Eq. (165) for three d ifferent NLO forward parton 
distributio ns : M RST98 [ |Mar98at (full curve), CTEQ5M [ |LaiOO[ (dashed curve), and GRV98 
NLO(MS) |G1u9^ (dashed-dotted curve). Also shown is the Ihs of Eq. (p^), i.e. the value k^/k" 
(dotted curve). 



the rhs of Eq. (165) for three current parametrizations of the NLO parton distributions. 
Fig. |l^ shows that the scale dependence drops out of the rhs of Eq. ( |165| ) , although the 
numerator and denominator separately clearly have a scale dependence. In addition, it is 
seen from Fig. |l9|, for the three parametrizations of the NLO parton distributions, that 
the relation of Eq. ( |165| ) is numerically verified to an accuracy at the one percent level! 
Both of these observations suggest that Eq. (165) holds and that the unpolarized valence u 
and d-quark forward distributions contain a non-trivial information about the anomalous 



magnetic moments of the proton and neutron. The relation of Eq. (164) can for example 
be checked in a theory with weak coupling as well as in chiral perturbation theory. Note 
that in a simple quark model with M^""' — 2 M2"°', one recovers from Eq. (165) the result 
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from SU(&) symmetry, i.e. k'p ^ - k". In quark model language, the relation of Eq 
implies that the breaking of the SU{6) symmetry follows some rule which is encoded in 
the valence quark distributions. 



Because Eq. (164) holds to such good accuracy, when applying it to the estimates 
given in Eqs. ( p^J162| ), it leads to the resuh that 2 J" + 2J'^ « M^ + M^. For the separate 



u and d contributions however, it was already noted from Table ^ that in the ansatz of 
Eq. (|157D, 2J" is quite different from M^, and similarl y 2J '^ is quite different from 



One therefore concludes that the valence ansatz of Eq. (157) for e''{x) leads to the result 
that the sum over all quarks yields the same contribution to the proton spin as to the 
proton momentum, at a scale as low as /i^ ~ 1 GeV ^. 

In Table |, we show the estimate of Eqs. (|l6l[jl63| ) for 2 J" + 2J'^ + 2 J" at different 
scales, using the three different NLO forward parton distributions discussed in Fig. 091 One 
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firstly sees that the three estimates are very close. One furthermore sees that J" + J'^ + J" 
decreases (very slowly) with increasing scale. The scale dependence of the total quark (J'') 
and gluon (J^) contributions to the nucleon spin has been calculated in Ref. [Ji96|. In 
particular, it has been shown that as /i 



J«(oo) 



1 



oo, there exists a fixed point solution 
3^/ 



2 16 + 3iV/ 



(166) 



and J» follows from J'' + js = 1/2. Applying Eq. ( |166| ) for Nf = 3, one obtains that the 
total quark contribution to the nucleon spin in the limit /i — )■ oo amounts to 36 %. 





2 J« 

(m^ = 2.5 GeV2) 


2 J9 

(/x2 = 10 GeV2) 


2 J« 

(m^ = 100 GeV2) 


u + d + 3 
(MRST98) 


0.60 


0.55 


0.51 


u + d -\- s 
(CTEQ5M) 


0.59 


0.55 


0.52 


u + d + s 
(GRV98 NLO) 


0.59 


0.55 


0.51 



Table 5 : Dependenc e of th e estimate for 2J" + 2 J"* + 27^ on tlie scale /i^ using the MRST98 
|Mar98a|], CTEQ5M |LaiOO|, and GRV98 NLO(MS) [|G1u98| forward parton distributions. 



One can further interpret the spin structure of the proton by decomposing the quark 
contribution to the proton spin in an intrinsic part (Ag) and an orbital part {L'^) |Ji98b|, 
for each quark flavor, as : 



J'(m) 



1 



(167) 



The flavor decomposed polarized quark contributions Aw, Ad and Aq to the proton spin 
have been determined by SMC iAde9^ at a scale = lo GeV^, and by HERMES |Ack99|] 
at a scale jj? = 2.5 GeV^ from the measurement of semi-inclusive spin asymmetries. Both 
experiments are in good agreem ent when extrapolated to the same scale /x^ = 2.5 GcV^, 
as has been done in Ref. | Ack99 |. The sum of the quark intrinsic spin contributions to the 
proton spin yields about 30 % at this scale. By using our estimate for J"^, and the measured 
values Aq for the different quark flavors, one can extract from Eq. (167) the corresponding 
orbital angular momentum contributions for the u, d, and s-quarks to the proton 
spin, as is shown in Table One sees from Table ^ that at the scale /i^ — 2.5 GcV^, 
our estimate for the total quark contribution to the proton spin yields about 60%, which 
is nearly entirely due to the w-quarks. The measured intrinsic spin contributions show 
that both Au and Ad are large, and have opposite sign. Our estimate leads then to the 
picture where the quark orbital angular momenta contribute about 30 % to the proton 
spin, of comparable size as the intrinsic spin contribution. For the separate flavors. Table ^ 
shows that for the u-quarks, the dominant contribution to the proton spin comes from the 
intrinsic contribution Au, with only a small amount due to orbital angular momentum. 
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2 J« ^ 2.5 GeV2) 


Aq {fi^ = 2.5 GeV2) 
(HERMES) 


2 LI ^ 2.5 GeV2) 


u 


0.61 


0.57 ± 0.04 


0.04 T 0.04 


d 


-0.05 


-0.25 ± 0.08 


0.20 T 0.08 


s 


0.04 


-0.01 ± 0.05 


0.05 T 0.05 


u + d -\- s 


0.60 


0.30 ± 0.10 


0.30 T 0.10 



Table 6: Estimate of 2L« at the scale /i^ = 2.5 GeV^ from Eq. ( |l67| ). The value of 2 J« is 
obtained according to the model estimate of Eqs. (161-163) using the MRST98 | Mar98a ] parton 
distributions at the same scale. The values for A q (inclu ding both quark and anti-quark contribu- 
tions) are taken f rom the HERMES experiment [AckQE] at this same scale. We added the errors 
from Ref. [Ack99| in quadrature. 



For the d-quarks however, the intrinsic contribution Ad is nearly entirely cancelled by a 
large orbital contribution, yielding a small net value for J'*. 

Having discussed a valence parametrization for the function e'^{x) up to now, we next 
investigate a possible sea quark contribution to e'(x). It is seen from Fig. || that the 
chiral quark soliton model calculation for the leading (at large Nc) isovector GPD i?" — E"^ 
contains such a sea quark component, implying that the valence model of Eq. (157) for 
e'^(x) is incomplete. This sea quark component to e'(x) is symmetric in x, in contrast 
to the sea quark contribution to q{x), which is antisymmetric in x. As the chiral quark 
soliton model calculation indicates that this sea quark component of e{x) is a very narrowly 
peaked function around a; = 0, we propose to parametrize it by a (5 function in x. Adding 



this piece to the valence parametrization of Eqs. (157) for u and rf-quark flavors, leads to 
the following parametrization for e'^{x) : 



e"(x) = A'^u^aiix) + B'^dix), 
e'^ix) = A'-d^aiix) + B^ix), 
e'ix) = 0, 



(168) 



where the par ame ters A", A are related to J", J through the total angular momentum 



sum rule Eq. (|158| ) as 



2 J? 



Ml 



(169) 



and where the parameters B", B"^ follow from the first sum rule Eq. (156) as : 



B" 



2^ 



2 J" - Mi 



2 J" - Af ° 



M. 



dva 



(170) 
(171) 



One sees that the total angular momenta carried by u- and c?-quarks, J" and J"^, enter now 



directly as fit parameters in the parametrization of Eq. (168), in contrast to the valence 



61 



model for e'^(x) of Eq. (157), where the values J" and J'^ are fixed as in Eqs. (161,162). 
Therefore, such a parametrization as in Eq. ( |168D can be used to see the sensitivity of hard 
electroproduction observables on J" and J"*, as will be shown further on. 

The physical interpretation of the sea quark part of e'(a;) in Eqs. ( |16^ ) can be un- 
derstood as being due to "vector meson exchange" because e'^{x) is normahzed to k'. In 
a vector meson dominance picture, the isovector (isoscalar) anomalous magnetic moments 
can be mainly understood as due to exchange of a p (w) meson, through its tensor coupling 
to the nucleon. 



Having specified the model for the forward distribution e'^{x), i.e. Eq. ( 16^ ), w e now 
turn to the ^-dependence of the GPD E^j^ix,^) = E%jj{x,^,t = 0) in Eg. |l54D. We 



generate the ^-dependence of El)p {x, £,) through a double distribution K'^{(3,a) as given 
by Eq. (|l|) similar as in Eq. (|l4§) : 



/•I /-i-l/JI 
Koi^^O = dp da 

J -I "'-1+1/31 



5{x- l3-a(,)K''{l3,a) . 



(172) 



The double distribution K'i{(3,a) is then obtained, in an analogous way as in Eq. (147) 
by multiplying the forward distribution e''(/3) with a profile function as : 



X«(/3,a) =/i(/3,a) e«(/3). 



(173) 



and where the profile function /i(/3, a) is taken as in Eq. (|148| ). Using the parametrization 
of Eq. ( |168D , one obtains then for the GPD £'|j^(x,^) the expression : 



r(26 + 2) 1 



22b+ir2(6 + i) 



1-^ 

e 



(174) 



where the first (second) term in Eq. (174) is the DD part originating from the valence 
(sea) contribution to e"^ respectively in Eq. (168). In the following, we take the power b in 
Eq. (174), which enters in the profile function in Eq. (173), equal to the value b — 1. This 
is motivated because this corresponds to "vector meson exchanges" with an asymptotic 

distribution amplitude. 

In Figs. |20| , pl| , we show the parametrization of Eq. (174) for E^j-, and Efyjy respec- 
tively at a fixed value ^ = 0.3. The two components in Eq. (174) are displayed for different 



values of J" [J"^), which are the free parameters in E^jj (E^jj) respectively. The valence 
contribution is calculated using the forward MRST98 valence quark distributions at a scale 
fi^ — 2 GeV^. It is seen that the sea contribution due to "vector meson exchange" displays 
the profile of an asymptotic distribution amplitude. One also sees from Figs. that 
the two components are correlated as the integral of their sum is fixed to k'. In particular, 
for the value J" « 0.32, the sea contr ibut ion to E" vanishes and the ansatz of Eq. ( |168| ) 
reduces to the valence ansatz of Eq. (157). Similarly, for the value J"^ ~ —0.03, the sea 
contribution to E'^ vanishes. For other values of J" ( J'^) however, the relative contribution 
of both components changes. As both components have a different shape in x, one gets 
a sensitivity on J" and J'' in observables where the GPDs E"^ and E'^ enter, as will be 
shown further on. 

In Fig. we show the ^ dependence of the three contributions of the mod el fo r 
i?". Besides the valence and sea contributions in the double distribution part of Eq. (172) 
which are shown at a fixed value of J" = 0.25, the third contribution originates from the 
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Figure 20: Double distribution part to the u-quark GPD JS" at t = using the ansatz of Eq. ( |l68|) , 
for different values of J" as indicated on the curves. Upper panel : valence part (VAL), lower 
panel ; sea contribution due to vector meson exchange (VM). 
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Figure 21: The same as Fig. G0| but for the double distribution part to the d-quark GPD E"^ at 
t = 0, for different values of J . 
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Figure 22: The three contributions of the model for the u-quark GPD at t = 0, for diff eren t 
values of ^ as indicated on the figure. Upper panel : valence part (VAL ) in the ansatz of Eq. (16S) 
for J" = 0.25; middle panel : sea contribution (second term in Eq. ( hOq)) d ue to vector meson 
exchange (VM) for J" = 0.25; lower panel : D-term contribution of Eq. (|153|) to . 



D-term of Eq. ( ]153| ), which contributes with opposite sign to as compared to H^, and 
therefore drops out of the spin sum rule. Notice also from Fig. ^ that for small ^, the sea 
contribution approaches a i5-function in our parametrization. 

Having specified the parametrization of (a; , ^, t = 0), we next discuss the i-dependence 
of the GPD E'^{x,^,t). This t-dependence of E'^ is constrained through the first sum rule 
of Eq. (^), which is given by the Pauli form factor F2{t). As stated before, we are only 
concerned here with hard exclusive reactions in the small —t region, with —t small com- 
pared to the hard scale Q^. In this small —t region, the simplest parametrization which is 
consistent with this sum rule constraint consists of a factorized ansatz 

E''{x,tt) = E'^ix,0 (i,Jo.7i)2 ^ (175) 

where the t-dependence has been expressed through a dipole form factor, which is known 
to be a good parametrization for the t-dependence of the nucleon magnetic form facto rs. 



In all estimates at small —t which contain i?', the factorized form (in t) of Eq. (175) is 



used. As discussed for the GPD H"^, a more realistic estimate of the i-dependence of the 
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GPD E"^ might also be given by a Regge ansatz analogous to Eqs. (135,152), which is 
postponed to a future work. 

To summarize this section, we constructed a three component model for the GPDs 
E"^ as illustrated in Fig. |2^. The valence component was fixed in terms of valence for- 
ward quark distributions. This assumption is supported by the calculation in the chiral 



quark-soliton model (see Sec. 3.2). For the sea part, we included the two physically im- 



portant contributions due to the D-term and the component associated with the "vector 
meson exchange" . The D-term contribution is required by the polynomiality condition 
(see Sec. 2^), whereas the "vector meson exchange" component is motivated by the vector 
meson dominance model for the nuclcon anomalous magnetic moments. 



4.2.4 Parametrization of the GPD H'' 

To model the GPD H"^, we nee d the corresponding polarized quark distributions. For 
this we follow the work of Ref. | Lea98 |, where a next to leading order QCD analysis of 
inclusive polarized deep- inelastic lepton-nucleon scattering was performed and which yields 
an excellent fit of the data. In this analysis, the input polarized densities (at a scale — 
1 GeV^) are given by : 



Aq{x,tx^) = ri,As S{x,^?), 



(176) 



where a SU{i) symmetric sea has been assumed, i.e. Ag = Au — Ad = As and S 
represents the total sea. On the rhs of Eqs. (176), the normalization factors j4„, A^, 
are determined so that the first moments of the polarized densities are given by % 
respectively. For the valence quark densities, rju and rjd were fixed in Ref. |Lea98| by the 
octet hyperon (3 decay constants which yield : 



riu = 



dx Auv{x) « 0.918 



m = 



r + l 

/ dx Advix) w -0.339 . (1 
"'0 



77) 



The first moment of the polarized sea quark density was determined by the fit of Ref. | Lea98| 
which gives rjq w —0.054. 

Our ansatz for H"^ is to generate it from a double distribution F'^{f3, a, t), analogous 
to Eq. (EtI), for which we take 



F«(/3,a,f) 



= h{(3,a) Auv{P)9l{t)/gm 
= h{p,a)Adv{P)g'i{t)/gi{Q) 



(178) 



where h{l3, a) is the profile function of Eq. ( |148| ). The parameter h in Eq. ( |148[ ) can again 
be chosen as a free parameter in such a parametrization, but we will show all estimates 
with the value 6 = 1 for the GPD H. 

One can check that Eq. (0) is verified by construction and that, using Eq. (|T^ ) to 
evaluate the axial form factors g\ and g\, the sum rule of Eq. ( [l0| ) is satisfied within 10 %. 

The parametrization of the antiquark contribution to H'^ remains to be investigated in 
future work. This might be interesting to study since a remarkable prediction of the chiral 
quark soliton model, noted first in Ref. [Dia96, Dia97|, is the strong flavour asymmetry 
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of polarized antiquaries. Such a feature is missing in other models like, for instance, pion 



cloud models, see discussion in Ref. [DresOObj. 

We stress that the factorization of the t-dependence of GPDs is not expected to hold 
when increasing the value of — i. For example, the calculations of the GPD H in the chiral 



quark soliton model |PenOOa| also show a deviations from the factorization ansatz for the 
f-dependence of GPDs as is shown in Fig. In the calculations presented here, we limit 
ourselves to relatively small values of —t and assume the factorization ansatz. 



4.2.5 Parametrization of the GPD E'^ 

As for the GPD i?', the GPD i?'^, corresponding with the induced pseudoscalar transition 
at the nucleon side in the matrix element of Eq. (|^) , also vanishes in the forward direction 
and does not contribute to DIS. However, it was shown in Ref. |Man99b, Fra99a, PenOOa 
that the pion exchange, which contributes to the region — ^ < a; < ^ of _B, may be non 
negligible at small t due to the proximity of the pion pole at t = m^. As discussed in 
Sec. 3.2.2, it was found [PenOOa in the chiral quark-soliton model calculation that in the 
limit t the function E exactly reduces to this pion pole contribution. Furthermore it 

is seen from Fig. that in the model calculation for E the pion pole part to E dominates 
over a wide range of t and ^ values. ^ ^ 

According to these findings, the main contribution to E is obtained by evaluating E 
assuming it is entirely due to the pion pole. According to this hypothesis, since the pion 
exchange is isovector, one has 



^ — 2 —pole 



(179) 



The t -dependence of ETr-poie{x,S,,t) is fixed by the sum rule of Eq. (|T^ in terms of the 
pseudoscalar form factor hA{t). 

In the region — ^ < a; < the quark and antiquark couple to the pion field of the nu- 
cleon. Therefore, this coupling should be proportional to the pion distribution amplitude. 
For the latter we adopt the asymptotic form. Expressing the quark's longitudinal momen- 
tum fraction z in the pion in the symmetric range — 1 < 2; < 1, the asymptotic distribution 
amplitude <^as is given by ^as{z) = 3/4 (1 — z"^), and is normahzed as J^^ dz ^as{z) = 1. 
The light-cone momentum fractions of the quark and antiquark in the pion are respec- 
tively given by (1 -I- x/S,)/2 and (1 — x/£,)/2. Therefore, E^^-poie is finally modelled in the 
estimates as : 



Et, 



'(-e<.T<0 hA{t)-<i>a 



(180) 



which satisfies the sum rule of Eq. (pi]) . 
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4.3 Deeply virtual Compton scattering (DVCS) 

In order to access the GPDs experimentally, the deeply virtual Compton scattering (DVCS) 



process 



7*{q)+N{p)^^iq')+N{p')., 



(181) 



was proposed as a practical tool in Ref. | lJi97b[ |. In the Bjorken regime, where the photon 
virtuality is large, with xb = l'2.{p ■ q) fixed, the leading order DVCS amphtude 
in the forward direction is given by the handbag diagrams of Fig. |l| [ |Mul94 |Rad964 
Ji97b|. These handbag diagrams express the factorization of the DVCS amplitude in a 
hard scattering part (which is exactly calculable in PQCD) and a soft, non-perturbative 
nucleon structure part, represented by the lower blobs in Fig. |l|. This soft, non-perturbative 
object is parametrized in terms of the GPDs as given by Eq. (||). The factorization of the 
DVCS amplitude in a hard and soft part has been proven in Ref. Ji98a, Rad98, Col99 for 
the leading power in Q and all logarithms. 

To calculate the DVCS amplitude, one starts from its definition as a nucleon matrix 
element of the T-product of two electromagnetic currents : 



= -i 



d'^x e 



i(q-x) 



(p'|T[J,'^^.(^)^c.m.(0)] \V) 



(182) 



where the four-vector index /i (y) refers t o the virtual (real) photon. The DVCS amplitude 
is obtained from the DVCS tensor of Eq. (182) by contracting with the photon polarization 
vectors as : 



T 



(183) 



where ef^{q) ( s^{q') ) are the polarization vectors of the virtual (real) photons respectively. 

In the following, we firstly introduce the kinemati cal v ariables entering the DVCS 
process. We then discuss the DVCS amplitude of Eq. (182) in the leading power in Q 
(twist-2 accuracy), and show its dependence upon the GPDs. In order to be able to extract 
the twist-2 GPDs from DVCS observables at accessible values of the hard scale Q, it is 
necessary to have an estimate of the effect of power suppressed (higher twist) contributions 
to those observables. The first power correction to the DVCS amplitude is of order 0{1/Q), 
hence it is called twist-3. We therefore discuss subsequently, the twist-3 corrections to 
the DVCS amplitude, which have been derived and calculated recently by several groups 
[ AniOOl , iPcnOObj |BelOOb| , |RadO0| , |Kiv0l4 |KivOlH |KivOl4 jRadOlaj iBelOOcj |Bel01b| using 
different approaches. We then explore the experimental opportunities and discuss the 
different DVCS observables. 



4.3.1 Kinematical variables of the DVCS process 



As explained in Sec. 2.1, the calculation of the DVCS amplitude in the Bjorken regime and 
the parametrization of the non-perturbative matrix elements, representing the lower blobs 
in the handbag diagrams, is performed in a frame where the virtual photon momentum 
and the average nucleon momentum are collinear along the z-axis and in opposite 
direction. To calculate the DVCS amplitude, it is therefore convenient to introduce lightlike 
vectors along the positive and negative z-directions as = P+/-\/2(l, 0, 0, 1) and n^^ — 
■ 1/a/2(1, 0, 0, — 1) respectively, satisfying p • n = 1. 
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In this frame, the physical momenta entering the DVCS process (181), where p{p') 
are the momenta of the initial (final) nucleon and q{q') are the momenta of the initial 



(final) photon respectively, have the following decomposition [GuiE 
— 1 ifi'^ 



2e F 
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(20 p^ + {Cfh^) n'^ + A 



-A'; 



(184) 

(185) 
(186) 
(187) 



where A_l is the perpendicular component of the momentum transfer A ( i.e. p ■ Aj 
n ■ Aj_ = ), and where the variables m^, ^ and ^ are given by 



2^ 



rriN 
P-q 



2e = 2e 



A2 

4 ' 

)2 - A2 



XB 



xb ' 
2 



(188) 
(189) 
(190) 



To twist-3 accuracy, Eqs. ( |l8^p^ ) reduce to 

P ^ P, A 



-2CP 



' A 



(191) 



4.3.2 Twist-2 DVCS tensor 



Using the parametrization of Eq. (|^) for the bilocal quark operator, the DVCS tensor in 
leading order in Q, ^dvcs-lo' follows from the handbag diagrams of Fig. |l] as ]Mul94 



Ji97b, Rad96a| 



TTIJ.U 

-"DVCS-LO 



+1 



dx C+{x,i) 



-Ely^,{x,^,t) N{p )zcj^'^'^N{p) 



2mN 



+1 



dx C {x, 



DVCS 



(a;,e,t) N{p)frf5Nip) 



- I A ■ n 

+ ^WCS(2^-^:*) N{p h5^N{p) 



,(192) 



where the symmetrical and antisymmetrical twist-2 tensors were introduced in Eq. (|62 
as : 



-g^ + n^p + n p^ 



(eM^)_ 



f-tivaprfp^ ■ 



(193) 
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In Eq. (192), we also introduced the coefficient functions C^{x,£,), which are defined as 



± 



£, — ie' 



(194) 



In the DVCS process on the proton, the GPDs for each quark flavor enter in the 
combination 



^DVCS 



4 
9 



9 



9 



(195) 



and similarly for H, E and E. The tensor for the DVCS process on the neutron is given 



by Eq. (192) by replacing H 



DVCS 



by 



Trn 

^DVCS 



9 



9 



(196) 



and similarly for H, E and E. Note that in Eqs. ( |195| , |196D the flavor dependent GPDs i/", 
H'^ and in our notation always refer to the corresponding quark flavor in the proton, 
e.g. if" = W^/P as explained in Sec. 4.2. 



One sees from the rhs of the DVCS amplitude of Eq. (|192| ) that the GPDs H, H, E, E 
enter in a convolution integral over the quark momentum fraction x. This is a qualitative 
difference compared to the case of DIS, where one is only sensitive (through the optical 
theorem) to the imaginary part of the forward double virtual Compton amplitude. In the 
case of DIS, the convolution integral collapses and the quark momentum fraction x is fixed 
to the kinematical quantity xb- In contrast, the non- forward DVCS amplitude of Eq. (192) 
displays a richer structure and has both real and imaginary parts. For the imaginary part, 
the convolution integral again collapses, and the value x is fixed to the kinematically de- 
termined skewedness variable ^. Therefore, the imaginary part of the DVCS amplitudes 
is directly proportional to the GPDs evaluated along the line x = (e.g. in Fig.Jl|), and 
measures in this way the 'envelope functions' H{£_, ^, t), E{^, ^, t), H{£,, t), and i?(^, ^, t). 

The x-dependence of the GPDs away from the line a; = ^ is contained in the principal 
value integral of the real part of the DVCS amplitude Eq. (192). As the variable x is not 
a kinematically accessible quantity, the extraction of GPDs away from the line x = f is a 
non-trivial task. The strategy which we propose her e is to start from physically motivated 
parametrizations for those GPDs as outlined in Sec. 12. One can then extract the param- 
eters entering those parametrizations from different observables, which have a different 
sensitivity to the real and imaginary part of the hard electroproduction amplitudes, as 
shown further on. 



4.3.3 Twist-3 DVCS tensor 

The DVCS amplitude of Eq. ( |192| ) displays a scaling behavior as it is independent of 
Q. One easily sees however that this amplitude is not complete when going beyond the 
leading order in Q. Although the twist-2 DVCS amplitude of Eq. (192), is exactly gauge 
invariant with respect to the virtual photon, i.e. -^dvcs-lo ~ electromagnetic gauge 
invariance is violated however by the real photon except in the forward direction. This 
violation of gauge invariance is a higher twist effect compared to the leading order term 
^Dvcs-Lo- Since ql H^ycs-LO ^ -(^i)'' ^dvcs-lc ^n improved DVCS amplitude 
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linear in has been proposed in Refs. | Gui9§| , Vdh9£] to restore gauge invariance (in the 
nonforward direction) in a heuristic way : 



^DVCS ~ ^DVCS-LO 



(197) 



leading to a correction term of higher order in Q to the twist-2 DVCS amplitude, because 
iP-q') = 0{Q^). 

Recently, the DVCS amplitude on the nuc leon to the order 0(1/(5) has been 
derived explicitely in a parton model approach jPenOOb and in a light-cone expansion 
framework |BelOOb|, yielding the result : 



dx 



i-gni- - 



(-9^1- 
(P-q) 



{P ■ q') 



[P ■ q') 



nPTp{x,i) C+(x,C) 



(198) 



(P . q') 



where the functions an d T^, are given by Eqs. ( p8| , p9D . 

In the expression Eq. (198) for the DVCS amplitude to the twist-3 accuracy, the first 
two terms correspond to the scattering of transversely polarized virtual photons. This 
part of the amplitude, containing npT^ and n^T^, depends only on the twist-2 GPDs 



H,E and H,E and was anticipated in Refs. |lGui98| , |Vdh99[ . The third term in Eq. ( |19S| ) 
corresponds to the contribution of the longitudinal polarization of the virtual photon. 
Defining the polarization vector of the virtual photon as 



(199) 



we can easily calculate the DVCS amplitude for longitudinal polarization of the virtual 
photon (L — > T transition), which is purely of twist-3 : 



{sl)^ H^^^ ='^-^ j ^dx ItI C+{x,S,)-zefT^k C-{x,0 



(200) 



It is therefore seen that this term depends only on new 'transverse' GPDs JF^ and JF^, 
which can be related to the twist-2 GPDs E, H and E with help of Wandzura-Wilczek 
relations as discussed in Sec. 2.5, and which are given by Eqs. (^8|^l|). 

The DVCS amplitude of Eq. (198) is electromagnetically gauge invariant, i.e. 



q^H^'' ^{q~A),H^^'' =0. 



(201) 



formally to the accuracy 1/Q^. In order to have 'absolute' transversality of the amplitude 
(i.e. such that Eq. (201) is satisfied exactly) we keep in the expression (|198D terms of the 



^^The DVCS amplitude to twi st-3 accuracy has also been derived for the analogous case of a pion target 
1 Refs. [ |AniOO[ |RadO0| , |Rad01a| . 
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order jQ'^^ by applying the prescription of Eq. (197), i.e. 



5x 



-51 



W7) 



(202) 



for the twist-3 terms in the amplitude. Formally such terms are beyond the twist-3 accuracy 
and they do not form a complete set of l/Q^ contributions, but we prefer to work with 
the DVCS amplitude, satisfying Eq. ( |201| ) exactly. 

A more systematic treatment of corrections of the order t/Q^ for DVCS still remains 
to be done. Recently, the target mass cor rections , which induce correction terms in powers 
of Af^/Q^, have been addressed in Ref. [ BelOlb . 



To calc ulat e the twist-3 DVCS amplitud e, the convolution of the leading order Wilson 
coefficie nts (194 ) with the WW kernels in Eq. ( 198 ) can be performed analytically, as shown 
in Refs. [ KivOlb] , BelOOcj. In particular, it has been shown t hat the strongest (integrable) 
singularity of the resulting integrand i s log arithmic only (see KivOlb for technical details). 

Finally, let us note that in Eq. (198), the quark flavor dependence in the GPDs J^, 
(Tyi) is implicitely understood, i.e. 



E 

q—u,d, 



(203) 



4.3.4 DVCS cross section and charge asymmetry 

In this section, we present results for the DVCS cross sections. We give all results for the 
invariant cross section of the ep — > ep7 reaction, which is differential with respect to Q^, 
XB, and out-of-plane angle <I> (<I> = 0° corresponds to the situation where the real photon 
is emitted in the same half plane as the leptons). The invariant ep ep^ cross section is 
given by : 



da 



1 



XB y 



dQ'^dxBdtd<i> (27r)4 32 



-1/2 



Tbh + TpYcS 



(204) 



where toat is the nucleon mass, y = {p-q)/{p-k), and k is the initial lepton four-momentum. 
In the ep —^ epj reaction, the final photon can be emitted either by the proton or by the 
lepton. The former process is referred to as the fully VCS process (amplitude Tpvcs in 



Eq. (204)), which includes the leptonic current. The process where the photon is emitted 



from the initial or final lepton is referred to as the Bethe-Heitler (BH) process (amplitude 
Tbh in Eq- (204)), and can be calculated exactly. For fu r ther te c hnical d etails of the 
ep — > ep7 reaction and observables, we refer to Refs. lGui98i |Bel0l4 |VdhOOb| 0. 



When calculating DVCS observables, we show all results with the exact expression 
for the BH amplitude, i.e. we do not expand the BH amplitude in powers l/Q- Also in the 
actual calculations, we use t he exact kinematics for the four-momenta of the participating 
particles as given by Eqs. (|l84Hl87| ). In this way we take (partially) kinematical higher 
twists into account. Furthermore, when refering to the twist-2 DVCS results, we include 
those higher twist effects which restore exact transversality of the amplitude as expressed 
through Eq. (202). Similarly, when referring to the twist-3 DVCS results (provided by a 
longitudinally polarized virtual photon), we include the gauge restoring higher twist terms. 



see Ref. | VdhOOa | for a calculation of the QED radiative corrections to the VCS process. 
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Because the BH amplitude contains two lepton electromagnetic couplings in contrast 
to the FVCS process, the interference between BH and FVCS process changes sign when 
comparing the e~^p — s- e+p7 and e^p — s- e~p7 processes. Therefore, in the difference of 
cross sections — dg- , the BH (whose amplitude is purely real) drops out. This difference 



measures the real part of the BH-FVCS interference [Bro72| 



Re T^Hj^FVCS* 



(205) 



and therefore is proportional to the real (principle value integral) part of the DVCS am- 
plitude. In this way, the difference ae+ — a^e- is sensitive to the GPDs away from the line 
a: = ^ (as e.g. in Fig. |l8|) . 

In the following, we firstly show results for DVCS cross sections and DVCS charge 
asymmetries (£e+ — ae-)/{cre+ + cfe-)- 

In Fig. the ^-dependence of the ep — s- ep7 cross section is shown for a lepton 
(either electron or positron) of Eg = 27 GeV (accessible at HERMES). The kinematics 
corresponds to the valence region {xb = 0.3) and to a ratio —t/Q^ — 0.1 (Remark that 
increasing the ratio —tjCf'^ increases the higher twist effects). 

It is firstly seen from Fig. |2^, that in these kinematics, the pure twist-2 DVCS pro- 
cess without D-term contribution (dashed curves) dominates the ep — > ep7 cross section 
compared to the BH process (which is sizeable only around <& = 0°, where it reaches its 
maximal value). In absence of the BH, the twist-2 DVCS cross section would give a $- 
independent cross section which practically saturates the ep ep7 cross section at $ = 
180°, where the BH is vanishingly small. The ^-dependence of the dashed curve in Fig. |2^ 
can be understood as the sum of this constant twist-2 DVCS cross section, the cross section 
for the BH process, and the interference of the BH with the relatively small real part of 



the DVCS amplitude (in the valence kinematics, xb — 0.3, shown in Fig. 23, the ratio of 
real to imaginary part of the DVCS amplitude without D-term, is around 15 %). 

When adding to the twist-2 DVCS amplitude the purely real D-term contribution, the 
resulting cross section for the full twist-2 DVCS process (dashed-dotted curves in Fig. [2^ ) 
shows a noticeable sensitivity to the D-term. At <& = 180° (where the BH 'contamination' 
is very small), the predominantly imaginary DVCS amplitude (in absence of the D-term 
contribution), and the purely real D-term amplitude have only a small interference. In this 
region, the DVCS cross section is enhanced by about 10 %, when using the chiral quark 



soliton model estimate of Eq. (|40| ) for the D-term Eq. (|23|). However, around $ — 0°, 
the purely real D-term amplitude interferes maximally with the BH. This interference is 
destructive for the electron reaction and constructive for the positron reaction. Therefore, 
the effect of the D-term can be very clearly seen in the ^-dependence of the charge asym- 
metry as shown on the lower panel of Fig. ^ Including the D-term contribution, the full 
twist-2 DVCS charge asymmetry changes sign and obtains a rather large value (« 0.15) 
at <& = 0°. On the other hand, at $ = 180°, where the interference with the BH (and 
hence the difference between e" and e+) is small, the charge asymmetry is correspondingly 
small. The pronounced ^-dependence of the charge asymmetry and its value at $ = 0°, 
provides therefore a useful observable to study the D-term contribution to the GPDs, and 



to check the chiral quark soliton model estimate of Eqs. (133). 

We next study the twist-3 effects on the DVCS cross section in the WW approxima- 
tion. It is seen from Fig. ^ (full curves) that they induce an additional (approximate) 
cos $ structure in the cross section. The twist-3 effects would induce an exact cos $ struc- 
ture only when the BH amplitude is approximated by its leading term in an expansion in 
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Figure 23: Invariant cross section for tlie e~p e~p7 reaction (upper panel) and DVCS charge 
asymmetry (lower panel) at Ee = 27 GeV, for the DVCS kinematics as indicated in the figure. 
Dotted curves : BH contribution; dashed curves : BH + twist-2 DVCS (without D-term); dashed- 
dotted curves : BH + twist-2 DVCS (with D-term); full curves ; BH -I- twist-2 + twist-3 DVCS 
(with D-term). The calculations are performed using the values bvai = bsea = 1. 



1/Q, neglecting its additional ^-dependence. In all calculations, we keep however the full 
^-dependence of the BH, due to the lepton propagators, which complicates the interfer- 
ence at the lower Q. One sees from Fig. ^ that the interference of the twist-3 amplitude 
with the BH -I- twist-2 DVCS amplitude is destructive for <i> < 90°, and constructive for 
$ > 90°. Around $ = 0°, the twist-3 effects reduce the full twist-2 DVCS cross section 
(including the D-term) by about 25 %, whereas around $ = 180°, they largely enhance 
the twist-2 DVCS cross section (by about 55 %). 

To see how the size of the twist-3 effects on the ep epj cross section decreases 
when increasing the value of at fixed xb and t, the cross section and charge asymmetry 
are shown in Fig. |^ at the same E^, xb, and t, as in Fig. but at a value = 5 GeV^. 
One firstly sees that going to higher Q^, at the same Ee, xb, and t, enhances the relative 
contribution of the BH process compared to the DVCS process. Consequently, the cross 
section follows much more the $-behavior of the BH process. For the twist-2 cross section, 
one sees again clearly the effect of the D-term which leads, through its interference with 
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the BH amplitude, to a charge asymmetry of opposite sign as compared to the one for the 
twist-2 DVCS process without D-term. 

The twist-3 effects in the kinematics of Fig. |2^, where to ratio t/Q^ is only half the 
value of Fig. |2^, are correspondingly smaller. It is furthermore seen that the twist-3 effects 
induce an (approximate) structure ~ A cos(2$) in the charge asymmetry [A k, - 0.08 in 
the kinematics of Fig. ^). Only for out-of-plane angles $ > 130° one sees a deviation 
from this simple structure induced by the twist-3 amplitude, due to the more complicated 
^-dependence when calculating the BH amplitude exactly. 
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Figure 24: Same as Fig. ||, but for = 5 GeV^ 



Besides exploiting the interference between the DVCS and BH contributions to the 
ep — > ep7 reaction, as e.g. through the charge asymmetry, one can also study the DVCS 
amplitude directly in a kine matical regime where it dominates over the BH amplitude. 
Indeed, it was found in Refs. | Gui9g , Vdh98 | that at some fixed Q^, xb and t, one favors 
the DVCS contribution over the BH by going to higher beam energy. In particular, at 
intermediate values of xb ^ 0.15 — 0.3, one is in such a regim e at 200 GeV beam energy 
at COMPASS, where a DVCS experiment has been proposed |d'Ho99|. 



74 



4.3.5 Electron single spin asymmetry (SSA) 

At intermediate lepton beam energies, one can also extract the imaginary part of the 
interference between the FVCS and BH amplitudes through the ep — > ep7 reaction with a 
polarized lepton beam, by measuring the out-of-planc angular dependence of the produced 
photon [ Kro96| . It was found in Rcfs. | Gui98| , Vdh98 that the resulting electron single 



spin asymmetry (SSA) 



^SSA ^ ^e,.=+l/2-ae,.^-l/2^ ^^^^^ 
CTe,/i=+l/2 + '^e,h=-l/2 



with ae,h the cross section for an electron of helicity h, can be sizeable for HERMES 
{Ee = 27 GeV) and JLab {Eg =4-11 GeV) beam energies. The helicity difference 
{<^e,h=+i/2 ~ <^e.h=~i/2) vanishcs for the BH process and depends linearly on the imaginary 
part of the DVCS amplitude. This helicity difference is therefore directly proportional 
to the GPDs along the line x — and maps out the 'envelope' function _ff(^,^,t), and 
analogously for E, H and E. 

In the following, we discuss the SSA to twist-3 accuracy and show the sensitivity of 
the SSA to different parametrizations of the GPDs in kinematics at HERMES { E^ = 27 
GeV ) and JLab ( E^ = 4.23 GeV ), where this SSA has been measured recently. 

In Fig. we show the ^-dependence of the DVCS cross section and of the single 
spin asymmetry (SSA) for the same values of and t as in Figs. ^ and ^ and for a 
value of xb = 0.15, accessible at HERMES. Note that, due to parity invariance, the SSA 
is odd in We therefore display only half of the $ range, i.e. $ between 0° and 180°. 

When comparing Figs. ^ with Fig. it is firstly seen that the cross sections 
for the ep epj process increases strongly when decreasing xb a.t fixed and fixed 
t, mainly due to the growth of the BH amplitude. Due to this large BH amplitude, the 



relative twist-3 effect on the cross section in Fig. 25 is smaller than the ones in Figs. |2J. 
However, the large BH amplitude leads to a large value for the SSA through its interference 
with the DVCS process. At (pure) twist-2 level, the SSA originates from the interference 
of the imaginary part of the DVCS amplitude and the (real) BH amplitude. In case the 
BH is approximated by its leading term in an expansion in 1/Q, the twist-2 SSA displays 
a pure sin $ structure. Due to the more complicated ^-dependence of the BH at the lower 
values of Q^, this form gets distorted and its maximum displaced. Note that, for practical 
considerations, our "twist-2" DVCS calculations include kinematical higher twist terms as 
well as the gauge restoring terms according to Eq. ( |202| ). Their effect can be seen in the 
slight change in the SSA (of the percent level), due to the DVCS process by itself (i.e. 
when increasing the real part of the DVCS amplitude by adding the D-term contribution, 
the curves for the SSAs are slightly displaced). 

We next discuss the twist-3 effects, calculated in WW approximation, on the SSA 
for the DVCS process. One sees from Fig. ^ that the twist-3 corrections induce an (ap- 
proximate) sin(2$) structure in the SSA. The amplitude of the sin(2$) term is however 
rather small and the twist-3 effects change the SSA by less than 5 % in the kinematics 
corresponding to t/Q^ = 0.1. It was checked that at xb = 0.3 and for a value t/Q^ = 0.1, 
the twist-3 effects on the SSA are of similar size. One therefore observes that although the 
twist-3 effects in WW approximation can provide a sizeable contribution to the real part 
of the amplitude (see Fig. [2^ ), they modify the imaginary part, and hence the SSA, to a 
much lesser extent. 
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Figure 25: Invariant cross section for the e~p e~py reaction (upper panel) and DVCS single 
spin asymmetry (lower panel) at Ee — 27 GeV, for the kinematics as indicated in the figure. 
Curve conventions as in Fig. E3. 



In Figs. 26 and we investigate the sensitivity of the SSA to different parametriza- 
tions of the GPD H in kinematics of a HERMES experiment |Air01| and of a JLab exper- 

) have been performed us- 



The calculations in the previous figures ( Figs. 23 
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uid01[. 


s. 


23 





ing the values h^ai — bsea — 1 for the parameters which enter in the profile function of 
Eq. (|l4|) to reconstruct the GPD H. In Figs. |26| and |7[ the SSA to twist-3 accuracy 
is shown for different values of hvai and bsea- One sees that increasing the value of h^ai 
and bsea decreases the size of the SSA. It is seen from Fig. |l^ that by increasing b^ai 
and bsea^ one approaches the forward quark distribution ( which corresponds to the limit 
bvai — bsea — oo ) . One may use this sensitivity to extract the parameters b^ai and bsea 
from accurate measurements of the SSA. For example, one sees from Fig. ^ that the max- 
imum of the SSA is reduced to about 80 % of its value when increasing the values of b^ai 
and bsea over the indicated range. 
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Figure 26: Sensit ivity of the single spin asymmetry for the e^p e^p^ reaction, in the HERMES 
kinematics [AirOl], to different values of the parameters b^q i and hsea as indicated, which enter 



in the parameterization of the GPDs //" and H ( see Sec. 4.2.1 
twist-3 effects to the DVCS amplitude. 
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Figure 27: Sensiti vity of the single spin asymmetry for the e~p — > 
kinematics |3tep01|, to different values of the parameters by gi and bsea 
in the parameterization of the GPDs //" and i/'* ( see Sec. 
twist-3 effects to the DVCS amplitude. 



4.2.1 



P7 reaction, in the CLAS 
as indicated, which enter 
All calculations include the 
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The previous figures for DVCS observables have all been obtained by neglecting the 
GPD E (except for its D-term contribution of Eq. 



(153)). In Fig. we investigate 
the effect of the GPD E on the DVCS SSA, using the three-component parametrization 

(VM) contribution ) of Eqs 



D-term, valence contribution, and "vector-meson" 



(154 



174). This parametrization allows to see directly the sensitivity of observables to the total 



angular momentum contribution of the w-quark J") and d-quark (J ) to the proton spin. 

In Fig. we show the SSA in kinematics accessible at HERMES for different values 
of J" (for a value J'^ — 0) corresponding to the values shown in Fig. |20[ Because the SSA 
is sensitive to the imaginary part of the DVCS amplitude (apart from a small contribution 
to the SSA due to the DVCS process alone, i.e. in absence of the BH process), it is mainly 
sensitive to the valence component of the GPD E^ (upper plot in Fig. Indeed, the 
D-term and VM contributions to the GPD E give rise to a real amplitude. By comparing 
Fig. (upper plot) and Fig. one sees that neglecting the contribution of the GPD E", 
corresponds approximately to the value J" ~ 0.2. Furthermore, one sees that varying the 
value of J" from J" = 0.1 to J" = 0.4 changes the SSA by about 10% of its value. Since the 
sensitivity of the DVCS SSA to J" (and J'') is modest, it calls for precise measurements as 
well as detailed studies of higher-twist effects and NLO corrections. However, the physical 
importance of J" and J'^ warrants such dedicated efforts. 



e + p— >e + P + Y 

E =27 GeV, (f = 5 GeV\ x„ = 0.30, t = -0.25 GeV^ 



0.35 




20 40 60 80 100 120 140 160 180 

* (deg) 



Figure 28: Sensitivity of the SSA for the e p ~* e p'y reaction, in HERMES kinematics, to 
different values of the quark contributions to the proton spin, J" and J"*, which enter in the 



parametrization of the GPDs _E" and E ( see Sec. 4.2.3 ). The dotted curve corresponds to 
the calculation when neglecting the contribution of the GPDs i?" and The calculations are 
performed for the BH + twist-2 DVCS amplitude. 

In Fig. |2^, we investigate the t-dependence of the DVCS cross section and SSA to 
twist-3 accuracy in CLAS kinematics. One sees from Fig. ^ that both cross section 
and SSA display a qualitative change with increasing values of —t (note that the SSA 
vanishes at i = tmin)- Furthermore, it is seen that when increasing the value of —t at 
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fixed and fixed Xb, the shape of the SSA changes from the "sin<&" dependence and 
adopts an increasingly larger "sin(2$)" twist-3 component. This qualitative change in 
the t-dependence is important to investigate experimentally to have confidence in the 
extraction of GPDs from a fit to SSA data. 
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Figure 29: Effect of the twist-3 corrections on the ^-dependence of the e~p e~p^ invariant 
cross section (upper panel) and SSA (lower panel) in CLAS kinematics, for different values of 
the momentum transfer —t as indicated in the figure. All calculations represent BH -I- twist-2 -|- 
twist-3 DVCS cross sections for bvai ~ bsea ~ 1.0. 

In Fig. we investigate the sensitivity of the DVCS SSA to different models for the 
t-dependence of the GPD H (we perform this comparison for the twist-2 DVCS amplitude). 



We compare the factorized model of Eq. (149) which has been used in all calculations before 
with an unfactorized, Regge inspired model of Eqs. (135 152[ ) (with a' = 0.8 GeV~^), which 
has an exponential i-dependence and exhibits shrinkage at small values of ^, as was shown 
in Fig. One sees from Fig. ^that at —t = 0.25 GeV^, the unfactorized form leads to 
a slightly reduced SSA compared to the factorized one. However, for a value of —t = 0.55 
GeV^ , one already sees a noticeable reduction of the unfactorized model over the factorized 
one, and the SSA gets reduced to about 65 % of its value when using an unfactorized ansatz. 
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Such important sensitivity clearly points out the need to map out experimentally the t- 
dependence of the SSA carefully in order to extract GPDs. It also stresses the need from 
the theoretical side to study the parametrization of the t-dependence of GPDs in more 
detail in the future. 
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Figure 30: Single spin asymmetry for the e p e py reaction in C LAS kinematics, for two 



models for the t-dependence of the GPD H. The factorized model of E g 
is compared to the unfactorized (Regge inspired) model of Eqs. ( 135 

different values of : 0.25 GeV^ (full curves), 0.4 GeV^ (dashed curves), 0.55 GeV^ (dashed- 
dotted curves). All calculations represent BH + twist-2 DVCS cross sections for bvai = bsea = 1-0. 
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4.4 DVCS with ^ A transition 



To derive the ^ A DVCS amplitudes, we use analogously as ioi N ^ N DVCS a frame 
where the virtual photon momentum and the average A^A momentum P'^, defined 
in Eq. (^), are coUinear along the z-axis and in opposite directions. By denoting the 
lightlike vector along the positive and negative z-directions as p'^ and respectively, we 
can decompose the physical momenta of the iV — + A DVCS process as follows : 



1 



ip'^+p'n 



9^ = -(2^ n^ 



P 



1 



A^ = p ^ - = - (2C) p'" 
g'^ = g'' - A^ = -2 - 
where the variables fh%/^, ^ and ^ are given by 
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2e 2e 



2e ) {m\ - m2 
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A^, (210) 



(211) 
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The leading order TV ^ A DVCS tensor follows from the handbag diagrams (Fig. |l| 
where the outgoing nucleon is replaced by a A"*"). By parametrizing the ^ A matrix 
elements of the vector and axial vector bifocal quark operators as in Eqs. ( p4[ ) and ( |90| ) 
respectively, the leading order N ^ A DVCS tensor H'^'^, defined in an analogous way as 
in Eq. (182), is given by : 



H^'" (j*p ^ 7A^ 



- [p^'n" + p-'n^' 



2 1 

3 6 



V7(p) {HM{x,tt) i-IC^i) 
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- [te-^'^^^n^px] 



dx C (x, £,) 
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A-- 
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C4x,^,t) — (A ■ Pnp ~ Ap)} N{p) , (214) 



m 



N 



where the factor 1/6 in Eq. (214) results from the quadratic quark charge combination 
iel-el)/2. 

To give estimates for the N ^ A DVCS amphtudes, we need a model for the N ^ A 
GPDs which appear in Eq. (214). Here we will be guided by the large Nc relations discussed 
in Sec. |^ These relations connect the ^ A GPDs Hm,Ci and C2 to the ^ iV 
isovector GPDs E"" ~ E'^, H"" - H'^ and E'^ ~ E'^ respectively, as expressed by Eq. 
All other (sub-dominant) GPDs, which vanish at leading order in the 1/Nc expansion, are 
set equal to zero in our estimates for the processes involving the N A GPDs. For the 
N N GPDs, we use the phenomenological ^-dependent ansatz as discussed in Sec. 4.2. 

For the N ^ A DVCS in the near forward direction, unlike the N ^ N DVCS case, 
the axial transitions (distributions Ci and C2) are numerically much more important than 
the vector transition (distribution Hm)- This is because the Hm conies with a momentum 
transfer (A) in the tensor /C^^ as seen from Eq. (p8|). In contrast, the GPD Ci, which is 
proportional to the polarized quark distribution in the forward limit (A — s- 0), enters with 
no momentum transfer in Eq. ( |l92| ). Besides Ci, the distribution C2 is numerically most 
important as it contains a pion-pole contribution. The tt*^ pole contribution to the — > A 
DVCS amplitude can be evaluated analytically by using : 



dx 



1 



1 



■ a 



2,77 —pole 



1_ ^/3 
2^ 4 



hA{t) , (215) 



with the pion-pole part of the pseudoscalar form factor h^it) given by Eq. (|4l|). This leads 
then for the tt*' pole contribution to the N ^ A DVCS tensor 



TT—pole 



{-f*p ^ 7A+) 



1 



n^Px 



9A 



V3 



-t + mi 



{p^{p)ApN{p) . (216) 



The N ^ A DVCS process can be accessed experimentally through the ep eA+7 
reaction. In Fig. we show our predictions for the fivefold differential ep eA+7 cross 
sections in the JLab energy range. As for the ep — > epj reaction, there is also for the 
ep eA+7 reaction an analog contribution when the photon is emitted from the electron 
line, which is the Bethe-Heitler (BH) process. In comparison with the BH for the ep epj 
reaction, the BH process in the ep — s- eA+7 reaction is reduced by an order of magnitude in 
the near forward direction. Firstly, this is due to the photon propagator in the BH process, 
which goes like 1/t, which is smaller for the N ^ A case due to the larger value of tmin in 
the unequal mass case compared to the N ^ N transition. Secondly, the electromagnetic 
TV — > A vertex in the BH process contains an additional momentum transfer (magnetic 
N A transition), which results in an additional suppression at the small angles. 

Fig. |3l| shows that if one stays away from the BH peaks by going in the opposite 
halfplane (negative angles on Fig. ^l|) as the electron lines, the ^ A DVCS process 
shows a fast increase with energy and starts to dominate over the BH already around 10 
GeV. 

As mentioned above, the A^ — > A DVCS in the near forward direction is dominated 
by the axial transitions Ci and C2. The vector transitions, which dominate the N ^ N 
DVCS are suppressed here by a momentum transfer in the corresponding tensors. Among 
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the axial transitions, the pion pole contribution in C2 dominates in the valence region. 
This is also seen in Fig. where the Ci contribution to the A DVCS amplitude is 

shown separately (at 12 GeV). 
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Figure 31: Angular distribution for the TV ^ A DVCS cross section at — 2.5 GeV^, xb = 
0.3, and for different beam energies accessible at JLab. A comparison is shown between the BH 
(dashed curves), the A'' — > A DVCS (dashed-dotted curves) and their coherent sum (full curves). 
We also show the contribution of the Ci GPD (dotted curve at 12 GeV). 



Due to dominance of the pion pole in the valence region, our predictions for A'^ — s- A 
DVCS are relatively model independent. Therefore, they can be used to estimate the back- 
ground for the DVCS on the proton as it may be difficult in an actual DVCS experiment to 
pin down the final state unambiguously. If one does not have the resolution to distinguish 
the ep —> ejA'^ — > e^ir^p reaction from the ep — > e^p reaction, an estimate of the former 
process may be helpful to quantify the background from the ep e7A+ contribution. 
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4.5 Hard meson electroproduction (HMP) 



In the previous sections, we discussed how the GPDs enter into different DVCS observables. 
As discussed in Sec. |[ the GPDs reflect the structure of the nucleon independently of the 
reaction which probes the nucleon. In this sense, they are universal quantities and can also 
be accessed, in different flavor combinations, through the hard exclusive electroproduction 



of mesons - tt 
given in Refs 



(see F ig. p[) for which a QCD factorization proof was 
According to Ref. |Col97|, the factorization applies when the vir- 
which corresponds to a small size configuration 
More technically, for a longitudinally polar- 



|Col97|. 

tual photon is longitudinally polarized, 
compared to a transversely polarized photon, 
ized photon the end-point contributions in the meson wave function are power suppressed. 
Furthermore, it was shown that the cross section for a transversely polarized photon is 
suppressed by compared to a longitudinally polarized photon. 

In the hard scattering regime, the leading order diagrams for meson production are 
shown in Fig. ^ (Th part in Fig. |2|) |^ Because the quark helicity is conserved in the 
hard scattering process, the meson acts as a helicity filter. In particular, the leading order 
perturbative QCD predicts [Col97| that the longitudinally polarized vector meson chan- 
nels ujl, 4>l) are sensitive only to the unpolarized GPDs {H and E) whereas the 
pseudo-scalar channels (tt'''^'^, 77, ...) are sensitive only to the polarized GPDs {H and E). 
In comparison to meson electroproduction reactions, we recall that DVCS depends at the 
same time on both the unpolarized {H and E) and polarized {H and E) GPDs. This 
property makes the hard meson electroproduction reactions complementary to the DVCS 
process, as it provides an a dditio nal tool to disentangle the different GPDs. 

It was shown in Ref. [|Die99|l that the leading twist contribution to exclusive electro- 
production of transversely polarized vectors mesons vanishes at all orders in perturbation 
theory. 





Figure 32: Leading order diagrams for the quark exchange contribution to hard meson electro- 
production. 



According to the above discussion, we give predictions only for the meson electro- 
production observables which involve a longitudinal virtual photon. The longitudinal 
j1 + p M + p two-body cross section daL/dt is given by : 



daL 

~dr 



1 



167r(s-m2,)A(s,-Q^m2,) 2 



JEEI-^''(^^^ = 0''^^;MI', (217) 



'i« h' 



For neutral isoscalar mesons with positive C-parity (e.g. /q, /2, ...), besides the quark diagrams 



shown in Fig. 
wavefunction IILehOi 



additional diagram contributes with 2 coUinear gluons projecting onto the meson 
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where is the amplitude for the production of a meson with hehcity Am 
tudinal photon, and where Hn 



hjq are the initial and final nucleon helicities. 



in Eq. (217), the standard kinematic function A(a;, y, z) is defined by 



A(a::, y, z) = + + — 2xy — 2xz — 2yz . 



: by a longi- 
Furthermore 



(218) 



which gives A(s, —Q^^mj^j) — 2mN\qL\, where Iq^l is the virtual photon momentum in the 
lab system. As the way to extract daL/dt from the measured fivefold electroproduction 
cross section is a matter of convention, all results for da^ / dt are given with the choice of 
the flux factor of Eq. (^). 



4.5.1 Hard electroproduction of vector mesons 

In this section, we consider the vector meson electroproduction processes 7£ + TV Vl + N 
at large on the nucleon A^, where Vl (/5° , pj, t^^L, ■•■) denotes the produced vector me- 
son with longitudinal polarization. The longitudinal polarization component of the vector 
meson is obtained experimentally through its decay angular distribution. p° electropro- 



duction data both at high W ( 30 <W< 140 G eV ) [ |AdlO0( | as well as at intermediate 
values oiW {^.S <W < 6.5 GeV ) [lAck99i |TytOl|| support the fact that s-channel helcity 
conservation (SCHC) holds to good accuracy. This implies that a longitudinally polarized 
vector meson originates from a longitudinally polarized photon. Therefore, when SCHC 
holds, one can access the leading order longitudinal photon cross section for vector meson 
production solely by measuring the decay angular distribution of the vector meson, i.e. 
without performing a Rosenbluth separation of longitudinal and transverse photon cross 
sections. 

In the valence region, the vector meson amplitude in leading order in Q is ob- 
tained from the hard scattering diagrams of Fig. ^ At small xb, and ui electro- 
prodution can also proceed through a two gluon exchange mechanism, as discussed in 
|Bro94, Fra96| , Rad96t] to which we refer for details. We restrict ourselves in the following 
discussion to the quark exchange contribution, and will mostly show results in the valence 
region. 

The resulting amplitudes M.^ for electroproduction of longitudinally polarized vector 



mesons were calculated in Rcfs. |Man98a, Vdh9S], yielding the following expression 



4 1 
9 Q 



dz 



X \av,n N{p) fiNip) + Bv,N N{p)ia^^'^^ N{p) \ , (219) 
L 2 lUN J 

where '^Vl (-2^) is the distribution amplitude (DA) for a longitudinally polarized vector 
meson, and where as is the strong coupling constant. 

For electroproduction on the proton, the amplitudes A and B in Eq. ( |2ig| ) are 
given by |VIan984 |Vdh98(| : 



^pIp 



j\xl.{e.H--e,H^) {-i 
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j dx-^ (e„ - ea 
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where = +2/3 (e^ = —1/3) are the u (d) quark charges respectively. 

For n cl cctroproduction on the proton, the ampUtudes A and B are given by 



| Man994 |Vdh99( | 



PL " 



J dx - E'^) 



X ~ ^ + ie X + ^ — ie 



x — ^ + ie X + S_ ~ ie 



(222) 
(223) 



For colp electroproduction on the proton, the amphtudes A and B are given by |[vlan98a, 
Vdh98| : 



A 



I^L P 



B.. 



ie X + S_ — ie 



(224) 
(225) 



One sees from Eqs. ( ^20[]225| ) that the ampUtudes for the hard electroproduction 
of longitudinally polarized p", and ui mesons involve different flavor combinations of 
the GPDs i/"(a;,^,t), H'^{x,^,t), and analogously for £:"(x,^,i), E'^{x,^,t). Therefore, 
measuring those different channels allows to make a flavor decomposition of the GPDs H'^ 
and E'i. 

In order to extract the GPDs from cross sections of hard meson electrop rodu ction, 
one needs to specify the meson distribution amplitude $y^(z) entering in Eq. (219). For 
the longitudinally polarized p'^ meson, an experimental estimate of the second Gegenbauer 
coefRcient of its distribution amplitude has been given for the first time in Ref. | CleOO |, 
and was found to be 03''^ = —0.1 ± 0.2 at a scale (/i^) = 21.2 GeV^. Unfortunately the 
precision in the determination of ai'^ is still too low to discriminate between different 



model predictions for this quantity (QCD sum rules: 0.18 ±0.1 [gall9§, 0.08 ± 0.02 
| Bak98 and instanton model: —0.14 | Pol99a all at a scale around /i « 1 GeV ). However 
all these analyses point to relatively small values for the second Gegenbauer coefficient 
and favor therefore a DA for the longitudinally polarized p meson that is rather close to 
its asymptotic form. 

In view of these findings, we use in all calculations shown below, the asymptotic DA 
for longitudinally polarized vector mesons : 



= /y 6z(l-z) 



(226) 



with fp w 0.216 GeV and f^, ~ 0.195 GeV, determined from the electromagnetic decay 

V e+e-. 

One sees from Eq. ( |219| ) that the leading order electroproduction amplitude for lon- 
gitudinally polarized vector mesons is of order 1/Q, in contrast to the leading order DVCS 
amplitude which is constant in Q. This difference is due to the additional gluon propagator 
in the hard scattering amplitude for the leading order meson electroproduction amplitude 
(see Fig. |3|). This also leads to the dependence of the leading order amplitude on the 
strong coupling constant as{^ Q^)- At large scales Q^, the running coupling constant 
is given by its expression from perturbative analyses. However, the average virtuality of 
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the exchanged gluon in the leading order meson electroproduction amphtudes can be con- 
siderably less than the external Q^, which is therefore not the "optimal" choice for the 
renormalization scale. This puts some caveats on the applicability of the leading order 
meson electroproduction amplitude to analyse absolute cross sections at accessible scales. 
The corrections to the leading order vector meson electroproduction amplitude (both in 
powers of l/Q and in a^) is still an open question to be addressed in future work. 

We next compare the results for electropr odu ction cross sections to the available 
data. Because the leading order amplitude of Eq. ( ^19[ ) is of order 1/Q, it predicts a l/Q^ 
scaling behavior for the cross section daL/dt. By measuring the behavior of the meson 
electroproduction cross section, one may study how fast one approaches the scaling regime 
predicted by perturbative QCD. In particular, the measurement of hard electroproduction 
reactions in the region « 1 — 20 GcV^, which is accessible experimentally, may pro- 
vide information on the importance of power corrections to the leading order amplitudes. 
One source of power corrections is evident from the structure of the matrix element of 
Eq. (^) which defines the GPDs, where the quarks are taken at zero transverse separation. 
This amounts to neglect, at leading order, the quark's transverse momentum compared 
to its large longitudinal (-1- component) momentum. A first study of these corrections 
due to the quark's intrinsic transverse momentum, assuming a gaussian dependence, has 
been performed in Ref. | |Vdh99 |. At virtualities of a few GcV^, these corrections were 
found to lead to a sizeable reduction of the cross section |Vdh99 |, e.g. a factor 4 reduc- 
tion at xb ~ 0.3 and « 4 GeV^ ( see Ref. [Vdh99| for details ). A comparison of 
these calculations to existing p° electroproduction data in the few GeV^ range is shown 
in Fig. |3|. 



As the 

| Der94 



IlP ^ 

AdIOOt 



Bre98| , 

system ^ ), the calculations in Fig. 



p^, P reaction has mostly been measured at small values of xb 
( or equivalently large values of the cm. energy W of the j*p 
for electroproduction can reveal how the valence 
region is approached, with decreasing value of W, where one is sensitive to the quark 
GPDs. For the purpose of this discussion, we call the mechanism proceeding through the 
quark GPDs of Fig. ^, the quark exchange mechanism (QEM). Besides the QEM, p° elec- 
troproduction at large and small Xb proceeds predomina ntly th r ough a perturbative 
two-gluon exchange mechanism (PTGEM) as studied in Ref. | Bro94 , Fra96 1 . To compare 
to the data at intermediate Q^, the power corrections due to the parton's intrinsic trans- 
verse momentum dependence were implemented in both mechanisms using the estimate 
of Ref. |Vdh99| . The comparison of the calculations with the data in Fig. ^ show that 
the PTGEM explains well the fast increase of the cross section at high cm. energy (W), 
but substantially underestimates the data at the lower energies. This is where the QEM 
is expected to contribute since xb is then in the valence region. The results including the 
QEM describe the change of behavior of the data at lower W quite nicely. 

Recently, data have also been obtained by the HERMES Collaboration for up 
to 5 GeV^ and around W ^ 5 GeV lAirOCfl . The comparison of the calculations for pi 
with those data is shown in Fig. One sees from Fig. |3^ a clear dominance of the QEM 
in the intermediate W range as predicted in Refs. |Vdh98, Vdh99|. When including the 
model estimate for the power corrections, a fairly good agreement with these longitudinal 
p° electroproduction data |AirOO is obtained, as seen from Fig. 

A dedicated experiment to investigate the onset of the scaling behavior of p° elec- 
troproduction in the valence region (Q^ ~ 3.5 GeV^ and xb ^ 0.3) is also planned in the 

Q2 (1 - xb)/xb- 



is expressed in terms of xb as : W 
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Figure 33: Longitudinal forward differential cross section for p'l electroproduction. Calculations 
compare the quark exchange mechanism (dotted curves) with the two-gluon exchange mechanism 
(dashed curves) and the sum of both (full curves). The calculations include the corrections due 



to the intrinsic transverse momentum dependence for the quark exchange I Vdh99 1 as well as 
for the two-glu on exchange |Fra9(:]. The data a re from NMC (triangles) Arn94|| , E6 65 (solid 
circles) ||Ada97i, ZEUS 93 (open circles) |Der95| ] and ZEUS 95 (open squares) |Bre98| ]. Figure 
from Ref. [|Vdh9E|]. 



very near future |]Guid9 j | using the CLAS detector at JLab. 

Besides the p° electroproduction cross section, the ujl electro product ion cross section 
has also been measured recently by the HERMES Collaboration TytOl| at intermediate 
values of ~ 5 GeV. The measurement of the ratio of lol/p^ electroproduction cross 
sections at these lower values of W allows one to see the departure from the diffractive 
regime, at small xb, where this cross section ratio is given by lol ^ Pl = 1 ^ 9. In the 
valence region, around xb ~ 0.3, the quark exchange mechanism of Fig. |^ leads to a ratio 
wl : p° ~ 1 • 5 |Col97 Vdh99|| , which is nearly twice as large as in the diffractive region. 
This larger ratio of around 0.2 for the uji^j p\ cross sections at larger values oi xb is in 
good agreement with the recent HERMES measurements [TytOl], and also supports the 
dominance of the quark exchange process of Fig. ^. 

It is also instructive to compare the t-dependences for vector meson electroproduc- 
tion cross sections. In particular we found that a factorized ansatz in t for the GPD as in 
Eq. (149) gives a much less steep f-dependence around xb — 0.1 as the unfactorized (Regge 
motivated) form of Eq. (152). In particular, the unfactorized form of Eq. (152) leads, in 
the region ~ 0.1, to an exponential t dependence daj^ ~ in the small —t region. 
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Figure 34: Longitudinal cross section for p^, electroproduction. Calculations compare the quark 
exchange mechanism (dashed curves) with the two-gluon exchange mechanism (dotted curves) and 
the sum of both (full curves). The calculations include the corrections due to intrinsic transverse 
moment um dependence as calculated in Re f. Vdh99|. The d ata are from E665 (open circles) 
|Ada97| and HERMES (sohd circles) ]AirOC|]. Figure from Ref. ]AirOO[. 



When using the value a' = 0.8 GeV^^ for the Regge slope parameter in Eq. (152), we 
found for the t-slope of the (twist-2) cross section daL for /9° electroproduction the value : 
BpO ~ 6.5 GeV~^. This value is remarkably close to the experimentally extracted slope 
parame ters for electroproduction on the proton which has been measured at HERMES 
| TytOl| in the kinematic range 0.7 < < 5.0 GeV^ and i.O < W < 6.0 GeV, and which 



was found to be |Tyt01| 



BpO = 7.08 ± 0.14 (stat.) tom (sys*-) GeV"^ . (227) 



This hints that a Regge type (unfactorized) form as in Eq. (|152|) might be a reasonable 
starting point for a more realistic parametrization of the t-dependence of GPDs. The t- 



slope of the vector meson electroproduction cross sections was also addressed in Ref. [ Fra96 



where it was shown that higher-twist contributions are able to explain the change of the 
t-slope with Q^. In particular the value of BpO decreases with increasing Q^. 

Although the data for the vector meson electroproduction channels at intermediate 
values of xb discussed above point towards the dominance of the quark exchange mecha- 
nism, it would be too premature at the present stage to try to extract quark GPDs from 
these vector meson electroproduction cross sections. To reach this goal, one first needs 
to get a better theoretical control over the power (higher- twist) corrections, which is an 
important topic for future work. 

Besides the cross section cti,, the second observable which involves only longitudinal 
amplitudes and which is therefore a leading order observable for hard exclusive meson elec- 
troproduction, is the transverse spin asymmetry, Av^n (TSA) for a proton target polarized 



89 



perpendicular to the reaction plane (or the equivalent recoil polarization observable) . For 
definiteness let us consider the following (azimuthal) asymmetry : 

\S±\ 

where /3 is the angle between the transverse proton spin S± and the reaction plane (spanned 
by the virtual photon and the produced meson). 

For the electroproduction of longitudinally polarized vector mesons, induced by a lon- 
gitudinal virtual photon, we found this transverse spin asymmetry to be given by ( compare 



with the corresponding Eq. (24C) for the pion ) 



^ _ lMAB*)/mN 

TT \Ami-e)~\B\^ {e+t/{4:m%))~Re{AB*)2^^' ^ ' 

which is proportional to the modulus |Aj^| of the perpendicular component of the mo- 
mentum transfer A of Eq. ( |186| ), and where A and B are given by Eqs. ( |22^p25| ) for the 



different vector meson channels. One sees that the transverse spin asymmetry is propor- 
tional to the imaginary part of the interference of the amplitudes A and B, which contain 
the GPDs H and E respectively. Therefore, it depends linearly on the GPD E. Note that 
in contrast, both in the DVCS cross sections and SSA as well as in the longitudinal cross 
sections for (longitudinally polarized) vector mesons, the GPD E only enters besides a 
large contribution of the GPD H . In order to increase the sensitivity to the GPD E in 
those observables, one needs to increase the value of the momentum transfer A, because 
the function E is kinematically suppressed in the amplitudes when A ^ 0. Since the value 
of A should remain small however in comparison with the hard scale Q in order not to be 
totally dominated by higher twist effects, the DVCS cross sections, SSA and vector meson 
cross sections give only a limited handle in the extraction of the GPD E. In this context. 



the transverse spin asymmetry of Eq. (229) provides a unique observable to extract the 
GPD E. 

Besides, also the theoretical uncertainties and open questions which we discussed 
for the meson electroproduction cross sections largely disappear for the transverse spin 
asymmetry. Indeed, because the transverse spin asymmetry involves a ratio of cross sec- 
tions, the dependence on the distribution amplitude and the strong coupling constant in 



Eq. (21£) drops out in this ratio. This also suggest that the transverse spin asymmetry is 



less sensitive to pre-asymptotic effects and that the leading order expression of Eq. (22£) 
is already accurate at accessible values of (in the range of a few GeV^). 

Due to its linear dependence on the GPD E, the transverse spin asymmetry for lon- 
gitudinally polarized vector mesons opens up the perspective to extract from it the total 
angular momentum contributions J" and J'' of the u— and d-quarks to the proton spin. 



Indeed, in the three-component parametrization for i?" and E'^ presented in Eqs. (154 
174), the values of J" and J'^ enter as free parameters (see Figs. |o|, |l]). We therefore 



investigate in the following the dependence of the transverse spin asymmetry for p^, ujl 
and electroproduction on J" and J''. Due to the different u- and d-quark content of the 
vector mesons, the asymmetries for the p^, ujl and channels are sensitive to different 



combinations of J" and J'^. From Eqs. (220-225), one finds that the p^ production is 
mainly sensitive to the combination 2 J" + J'*, ujl to the combination 2 J" — J'^, and 
to the isovector combination — J'^. 
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Figure 35: Transverse spin asymmetry for the 7|,p — > p5,P reaction sXt — -0.5 GeV^ and = 5 
GeV '^ . T he estimates are given using the three-component model for and E"^ as described in 
Sec. i.2.3, where J was fixed to the value J =0. The curves show the sensitivity to the value 



of J" as indicated on the curves, which enters in the model for the GPD S", shown in Fig. |2C|. 
For the GPDs ff", H"^ and for the valence part of E^ , E"^, the values byai ~ &sca ~ 1 have been 
used. 
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Figure 36: Same as Fig. ^ but for the jlp ^ lulP reaction. 
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In Figs. ^ and [3^ we firstly show the xb dependence of the transverse spin asym- 
metry for p° and lul production at = 5 GeV^ and —t = 0.5 GeV^. In these figures, 
we fixed the value of J'' = 0, which is close to the valence estimate of Eq. (162), for which 
the VM part of E"^ vanishes (see Fig. ^T]). For this value of J'', we then demonstrate in 
Figs. the sensitivity of the p° and lol transverse spin asymmetries to different values 

of J". Note that the value J" « 0.3 would correspond to the valence estimate of Eq. (161), 
as is also seen from the vanishing of the VM part of E" in Fig. ^ for this value 



From Figs. ^ 36, one sees that the transverse spin asymmetries for pi and ul elec- 
troproduction display a pronounced sensitivity to J" around ccs « 0.4, where asymmetries 
are predicted in the -15 % to -30 % range according to the value of J". This transverse 
spin asymmetry for p° and at large xb gets enhanced by the (isoscalar) D-term, which 
cont ribut es to both processes. In particular, it gives a large real part to the amplitude B in 
Eq. ( |229| ) which is then multiplied by the imaginary part of the amplitude A, proportional 
to iJ(^,^,t). At these larger values of Xb a. precise extraction of the values of J„ and Jd 
from the p° and ujl asymmetries requires a good knowledge of the D-term which can be 
measured e.g. through the charge asymmetry for DVCS as discussed in Sec. 1.3.4 . At 
smaller values of a;^, where the D-term contribution shrinks away, the sensitivity to J" 
and J'^ arises through both the real and imaginary parts oi B, which depend on through 
the VM contribution to i?" (giving a purely real part to B) and the valence contribution 
to i?" (giving predominantly an imaginary contribution to B). 

For comparison, we show in Fig. ^ the transverse spin asymmetry for at the 
smaller value of —t = 0.25 GeV^ (the asymmetry of Eq. ( ^29[ ) is proportional to 
and = 2.5 GeV^. Those values are at present already accessible at HERMES. It will 
therefore be very interesting to provide a first measurement of this asymmetry in the near 
future for a transversely polarized target. 
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Figure 37: Same as Fig. BSl but for t 



-0.25 GeV^ and 



2.5 GeV^ 



Although we showed the transverse spin asymmetries in Figs. 35 3^, for the value 
J'^ — Q, one can vary this value and exploit the different dependence of the p^, (sensitive to 
2J'"' + J'^) and ojl (sensitive to 2 J" — J'*) transverse spin asymmetries to extract information 
on the different flavor contributions. 



92 



Yl + P ^ Pi^ + n 



< 

z 
a! 

00 

Z 
< 

H 



0.6 



0.4 



0.2 



-0.2 



-0.4 



-0.6 











/' = 0.3 


= 2.5 CeV^ 






-t = 0.25 GeV^ 






...^/=-0.1 






/ = o.o^ 






/ = 0.1 /I 


, , , , 


.... 1 .... 1 . 





0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Figure 38: Transverse spin asymmetry for the "flp ^ p^n reaction at t — -0.25 GeV^ and — 
2.5 GeV^ . The estimates are given using the three-component model for and E"^, where J" 
was fixed to the value J" — 0.3. The curves show the sensitivity to the value of J as indicated 
on the curves, which enters in the model for the GPD E'', shown in Fig. plj 
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Figure 39: Transverse spin asymmetry for the j^p — > pjn reaction at = 5 GeV^ for different 
values of t, as indicated on the curves (in GeV^). The estimates are given using the three- 
component model for i?" and E'', where J" and J"* were fixed as indicated on the figure. 
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One is directly sensitive to the difference J" — J'^ by measuring the transverse spin 
asymmetry for p~j^ production. The corresponding asymmetry is shown in Fig. ^ for a 
value of J" = 0.3 and for different values of J'^. One observes that this asymmetry can 
be very large over a wide range in xb- Because for pj, the isoscalar D-term contribution 
is absent, this large value for the asymmetry mainly originates from the product of the 
imaginary part of the amplitude A ( proportional to (C, t) ) with the real part of the 
amplitude B ( mainly proportional to the VM part of i?" — E'^ which contains the sensi- 
tivity to J" - J'' ). 

In Fig. |3^, we also show the t-dependence of the transverse spin asymmetry. One 
observes how by increasing the value of —t, the kinematical accessible region over which 
one can access this asymmetry gets larger. 

Summarizing this section, we have identified the transverse spin asymmetry for lon- 
gitudinally polarized vector mesons as a unique observable which depends linearly on the 
GPD E. These asymmetries for p^, ojl and display a large sensitivity to the total 
angular momentum contributions J" and J'^ of the u- and d-quarks to the proton spin. As 
these asymmetries represent a ratio of cross sections, several theoretical uncertainties drop 
out. Also from the experimental point of view, such ratios are much easier to measure 
than absolute cross sections. A first measurement of these asymmetries will therefore be 
of high interest in the quest of our understanding of the quark contributions to the proton 
spin. 

4.5.2 Hard electroproduction of pions 

In this section, we discuss the hard pion electroproduction processes 7|, + p — > ti"'" + n and 
7l + P ^ 7r° -hp. 

The pion electroproduction amplitude Ai^ for a longitudinal virtual photon is also 
obtained in leading order in Q from the hard scattering diagrams of Fig. This leads to 
the amplitude : 



Mi 



4 1 
9 Q 



A • n 



X <^ A^N N{p ) N{p) + B^N N{p ) 75- N{p) 

IrriN 



(230) 



where ^niz) is the pion distribution amplitude, and where as is the strong coupling con- 
stant. 



For 7r°p, the amplitudes A and B in Eq. ( |230[) are given by [|Man98at |Vdh98 | 
. _ . . ^ , „„ , ^ 1 

^ttO p — 



B-^Op 



j\.l={e.E--e.E^) [-1 



(231) 
(232) 



+ ie X + ^ — ie 

For TT"*" n electroproduction on the proton, the amplitudes A and B are given by 



[ Man99b| , |Fra99a | 
A j- 



dxin^ - H") \ — — + — — 



(233) 
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dx - E'^ 



ed 



(234) 



In the amplitudes for pfon electroproductfon of Eq. (230), the pfon distribution am- 
phtude <I>7r enters. Recent data |Gro98| for the 7r"7*7 transition form factor up to Q'^ = 9 
GeV^ are in agreement |^ with the asymptotic form of the distribution amphtude, given 
by: 

$,(z) = V2/, 6z(l-z), (235) 

with /tt = 0.0924 GeV fr om th e pion weak decay. 

As discussed in Sec. 2.4.2, the function E^ — E^ contains a strong pion pole singularity 



in the limit t — s- to^. The pion pole gives a dominant contribution to the hard 7r+ elec- 
troproduction amplitude ■A4^+ in the valence region, which can be worked out explicitely. 
Using the pion pole formula of Eq. ( ^l|) for the induced pseudoscalar form factor hA{t), 
and by using the PCAC relation qaI fir = g-wNN /it^Ni where Qt^nn is the ttNN coupling 
constant, one obtains for the pion pole part of the amplitude ■ 



ML 



9-itNN 

-t + mi 



Nip HNip) 



(236) 



where i^V represents the pion electromagnetic form factor. When using an asymptotic 
distribution amplitude for the pion, the leading order pion pole amplitude is obtained by 
using in Eq. (236) the asymptotic pion form factor F^", which is given by : 



167ras 



(237) 



Besides the tt hard electroproduction, one can also study the r] and rf hard electro- 
production. For the r/ electroproduction, when neglecting the effects of the QCD axial 



anomaly, the corresponding amplitudes A and B as in Eq. (230) are given by |Man98a 



A. 



dx . — ( 6^ 



i/" + i?" - 2 es H' 



1 



1 



Br, 



1 / 



dx ^{eu E" + ed E"" - 2e, E' 
-1 V6 ^ 



X — ^ + ie X + — ie 
1 1 



X ^ + ie X + — ie 



, (238) 
(239) 



The effects of the QCD axial anomaly have been studied in Ref. Eid9S|] . 

We show in Fig. ^ the leading order cross sections for various pseudoscalar meson 
production channels. As discussed before for the absolute cross sections for vector meson 
electroproduction, also for pseudoscalar meson production the (higher-twist) power cor- 
rections to the cross section have not yet been worked out systematicall y. Ther efore, these 
leading order results could receive sizeable corrections (see e.g. Ref. | Vdh99[ for a first 
estimate of such corrections). To be on the safe side, we therefore show in Fig. ^ the 
leading order predictions at a relatively large scale (Q^ = 10 GeV^), where we use for the 
running coupling constant as its expression from perturbative analyis at a scale = 10 
GeV2. 



^^See however Ref. |Bak01| for a recent discussion of the deviations from the asymptotic distribution 
amphtude for the pion. 
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By comparing in Fig. 40, the tt+tt, and n^p channels, one immediately sees the promi- 



nent contribution of the charged pion pole to the tt"*" cross section. The cross section for 
the TT*^ channel, where this pion pole is absent, is more than a decade lower for xb > 0.15. 
Going to smaller values of xb, on the other hand, one finds an increasing ratio tt" : tt"*", due 



to the amplitudes A of Eqs. (231,233). This was also noticed in Ref. |Man99t|. One can 



also understand the ratio 7r° : 7r+ at smaller xb by comparing the amplitude of Eqs. (231) 



and (233). One sees that in the tt^p channel, the GPDs and H'^ enter with the same 
sign, whereas in the 7r+n channel, they enter with opposite signs. These GPDs and H'^ 
are constructed starting from the polarized quark distributions Au and Ac? as discussed 
in Sec. 4.2.4. Because Au and Ad have opposite signs, this leads to a partial cancellation 
for the channel compared to the tt^ti channel, yielding a ratio tt" : 7r+ « 1 : 5 for the 



non-pole part of the 7r+ channel |Vdh99|. 

For the rjp channel, using the ansatz for the GPDs H"^ discussed in Sec. 4.2.4 in 
Eqs. (238,239), a slightly smaller cross section was found as for tt", i.e. 77 : tt" « 4 : 5 
I Vdh99|, comparable with the ratio 77 : tt'^ « 2 : 3 found in Ref. [Man98a]. This ratio is also 



very sensitive to the x-dependence of the polarized u- and d-quark distributions, giving a 
possibility to cross-check various parametrizations for those polarized quark distributions. 
Besides, when including the effects of the axial anomaly and SU (3) breaking, markedly 
different results were obtained in Ref. | |Eid99| |, i.e. tt" : 77 : 77' = 27 : 1.6 : 1.4, which will 
be interesting to check experimentally. 



YlP— > TtN, n A , KY 

= JO GeV\ t = -0.3 GeV^ 




Figure 40: Leading order predictions for the iiN , ttA and KY longitudinal electroproduction 
cross sections aXt — -0.3 GeV^, as function oi xb (plotted up to the xb value for which t — tmin)- 
Results for pion and charged kaon channels are given using an asy mptotic distribution amplitude. 
For K^T,^ , predictions are shown using the Chernvak-Z hitnitsk v |Cher84| distribution amplitude 
with antisymmetric part : rj"^ = 0.25. Figure from Ref. | FraOOa ] . 



As we discussed already for the hard electroproduction of longitudinally polarized 
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vector mesons there is, besides the cross section ctl, a second relevant observable, which 
involves only longitudinal amplitudes and which is of leading order. This is the transverse 
spin asymmetry for a proton target polarized perpendicular to the reaction plane, see 



Eq. (22? 



For the hard electroproduction of ttN final states, the transverse spin asymmetry 



Attn is given by | Fra99a 
2|A| I 



N 



(240) 



where A and B are defined as in Eqs. (231-234) 



>H 0.7 

Pi 

W 0.6 



* + + 

Yl p — > 71 n, 71 A 



< 



0.5 
0.4 



« 0.3 



0.2 
0.1 

-0.1 
-0.2 









71^ n 








■V;;^ - -.,-0.5 






-0.1 \ 


-0.3 




-0.2 J 


-0.35 / 


/o.s 








7t A 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Figure 41: Transverse spin asymmetry Attn for the longitudinal electroproductio n of n'^ n and 
7r'''A°, at different values of t (indicated on the curves in GeV^). Figure from Ref. [FraOOa|. 



In Fig. ^ (upper plot), we display the transverse spin asymmetry A^^+n foi' 7r^?i 
electroproduction at several values of — i, as a function oi xb- Fo r the GPDs if " and 77'*, 
we use the phenomenological parametrization as discussed in Sec. 4.2.4. It is obvious from 
Fig. |4^ that the predicted -4^+„ is large and comparable with the res ults of [ Fra99a | in 
which the GPDs computed in the chiral quark soliton model |PcnOOa| were used. As a 
consequence, the investigation of these processes can provide unique tests of the soliton 
type approach to baryon structure. 

The transverse spin asymmetries are likely to be less sensitive to higher twist effects 
and to NLO (in as) corrections. The latter corrections have been calculated recently in 
Ref. |Bel01d| and were found to be sizeable as well as strongly scheme dependent for ab- 
solute cross sections, but nearly drop out of the transverse spin asymmetries. Hence these 
asymmetries can be explored already at presently accessible scales at HERMES and JLab, 
using a transversely polarized target. 

The HERMES collaboration reported already preliminary a rather large asymme- 
try for 7r+ electroproduction in case of a longitudinally polarized target along the beam 
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direction |Tho01 . As was noted in Ref. |PolOO|, such a measurement includes a small 
target polarization component transverse to the virtual photon direction Sj_, leading to 
the following form for the polarized cross section ap : 



C7p sin <i> 



L<^LT 



(241) 



where cr^ receives only contributions from the longitudinal amplitude which dominates in 
the hard regime for the exclusive channel. The cross section cr^-p contains the interference 
of longitudinal and transverse amplitudes, which is suppressed by a power of the hard 
scale (1/Q) relative to cr^. Note, however that the transverse spin component S± drops 
as 1/Q relative to the longitudinal spin component for a target polarized along the 
beam direction, so that in this case both terms in Eq. (241) are of the same order in 1/Q. 
Therefore, to extract the transverse spin asymmetry unambiguously, an experiment using 
a transversely polarized target is really needed. 

A measurement of At^n can also provide an important help in the extraction of the 
pion electromagnetic form factor. For xb < 0.15 where one reaches values tmin ^ 2m^, we 
find that the pion pole constitutes about 70 % to the longitudinal cross section. Measure- 
ment of the transverse spin asymmetry, which is an interference between the amplitudes 
A and B in Eq. 23C , would help to constrain the non-pole term, and in this way help to 
get a more reliable extraction of the pion form factor. 



4.5.3 Hard electroproduction of strangeness 



In this section, wc turn to strangeness hard electroproduction reactions + p K + Y . 
In leading order, the kaon electroproduction amplitude is given by an expression analogous 
as Eq. ( ^ . 

For the charged kaon electroproduction channels, the amplitudes A and B are given 
by ||Fra99a| |FraOOa[ : 



A 



K+Y 



dx HP- 



B 



K+Y — 



-1 
1 



dxEP- 



X 



(242) 
(243) 



where eg is the s-quark charge, and Y is the produced hyperon {Y = A, E"). To estimate 
the GPD HP^^, we use the SU{3) relations as discussed in Sec. IJ. The SPD EP^^ 
contains a charged kaon pole contribution, given by : 



B 



pole 
K+Y 



-3/(20 Vk fK 9KNY i2mN)/{-t + m\) , 



(244) 



where Jk — 159 MeV is the kaon decay constant, and where gKNY are the KNY coupling 
constants. In line with our use of 5[/(3) relations to estimate H, we use also SU{i) 
predictions for the coupling constants : gKNA/ ~ —3.75 and qknt,/ ~ 1.09, which 
are compatible with those obtained from a Regge fit to high energy kaon photoproduction 
data [Guid97]. In Eq. (^44|) , ryf^ is defined as : 



vk 



dc 



C '^k (0 



1 



(245) 



98 



where (C) is the symmetric (in part of the distribution ampHtude (DA) of the kaon, 
which runs over the range -1 to +1 in the notation of Eq. (245). One has Vk — ^ ^^r an 
asymptoti c DA ( $"(0 = 3/4(1 - C^) ), and r;|- = 7/5 for the Chernyak-Zhitnitsky (CZ) 

we have also defined 



kaon DA |Cher84]. Furthermore, in Eq. ( p45D we have also defined - for further use - 77^ 
for the antisymmetric part ^^iO of the kaon DA. The latter is due to SU{3) / symmetry 
breaking effects. 

For the i^^S"*" electroproduction, the expressions for A and B, all owing f or both a 
symmetric and antisymmetric component in the kaon DA, are given by jFraOOal : 




.(246) 



In contrast to tt'^ electroproduction, K'^ electroproduction can contain a pole contribution, 
which is given by : 



B 



pole 



fK9KNT. (2™Af) 



(247) 



and which vanishes when the kaon DA is symmetric (i.e. when 'q'^ = 0). Therefore, the 
pole contribution to Bi(q^+, provides a direct measure of the antisymmetric component of 
the kaon DA. 

In Fig. we show besides the leading order predictions for the pion electroproduc- 
tion cross sections also the leading order predictions for kaon hard electroproduction cross 
sections at = 10 GeV^. The charged pion and kaon channels obtain a large contribution 
in the range xb ^0.1 from the pion (kaon) pole. This largely determines the ratio between 
these channels at larger xb- For values of —t in the range 0.1 ^ 0.5 GeV^, this yields 
TT+n : A'+A « 7 : 1 ^ 1.8 : 1, using an asymptotic DA for both vr and K. The kaon DA is 
not well known however, and the results with a CZ kaon DA yield cross sections, for 
the pole contribution, larger by a factor (7/5)^ ~ 3.8. The ratio iC+A : K^YP at large xb 
is determined from the ratio of the couplings : Qkni^/ Qknt. ~ 12. For the Ar°I]+ channel, 
the pole contribution is absent if 77^ = (as for 7r°). In this case, the ratio ir^p : Ar°E+ 
is determined by the amplitude A and is very sensitive to the input valence quark distri- 
bution into H. For Auy « — Adv expected in the large Nc limit, 7r° : « 1 : 3, while 
for Auy « -2Mv prefered by the global fit to DIS of Ref. gca9|, 71° : AT*' « 3 : 1. The 
sensitivity of this ratio to the polarized quark distributions might be interesting to provide 
cross-checks on such global fits from DIS In Fig. 1^, we show the results for by 
using the polarized distributions of | Lca98| as input for H (as discussed in Sec. 1.2.4 ). 

Besides the contributi on of the amplitude A, Ar^E^ electroproduction has also a pole 
contribution, given by Eq. (247), which is nonzero if ry^ 7^ 0. In the estimates shown here, 
we include the pole contribution of Eq. (247), and use the CZ kaon DA with = 0.25. 
The resulting AT^ pole contribution provides a sizeable enhancement of the K^Yi'^ cross 
section (it gives roughly half the value of the amplitude A at the largest xb)- 

After the cross section, we next discuss the transverse spin asmmetry Aky for the 



strangeness channels, which is given by |FraOOa| 
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Figure 42: Transverse spin asymmetry, for K^A, K'^T,'^ and if^E^ longitudinal electroproduc- 
tion for different values of t (indicated on the curves in GeV^). For K+A and K+T,°, the thick 
(thin) lines are the predictions with the asymptotic (Chernyak-Zhitnitsky) kaon distribution am- 
plitude respectively. For Tf^E"*", the thick (thin) lines are the predictions with Chernyak-Zhitnitsky 
type kaon DA, with antisymmetric part : rj'^ = 0.25 (0.1). 



where A and B are as defined before. It is interesting to note that in the case of hyperon 
production, the same transverse spin asymmetry can be measured on an unpolarized target 
by measuring the polarization of the recoiling hyperon through its decay angular distribu- 
tion. The transverse spin asymmetries Ak+a^ -^k+i:" and Akos+ are shown in Fig. A 
comparison with Fig. ^ for the ir^n channel shows that the transverse spin asymmetries 
for the charged kaon channels are as large as for 7r+n. One also sees that Ak+ho ~ Ak+a- 
For if ° production, the sensitivity to the SU (3) / symmetry breaking effects in the kaon 
DA is illustrated (lower panel of Fig. by plotting Ako-s+ for two values of 77^. Because 
-^K"s+ is directly proportional to 77^, it provides a very sensitive observable to extract the 
if" form factor. 

As we discussed before, the hard electroproduction of strangeness also gives an unique 
access to the second class currents (e.g. weak electricity), for a first estimate see Sec. 21. 

To summarize this section, we have shown that the yields for hard exclusive electro- 
production reactions of decuplet and octet baryons are similar. Strange and nonstrange 
channels can be comparable and in some cases the strange channels can even dominate (de- 
pending on the distribution amplitude and polarized parton distributions), in contrast to 
low-energy strangeness production. Large transverse spin asymmetries are also predicted 
for these reactions. 
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4.5.4 Hard electroproduction of ttA final states 



In this section, we discuss the hard electroproduction of vrA final states : 7^ + p ^ tt + A. 
In Sec. 2.7, we introduced the GPDs for the iV ^ A transition, and have seen 



how they can be expressed in terms of ^ GPDs through Eq. (|93D, using large 
relations. The large Nc relations make definite predictions for the ratios between the 
different cross sections for charged pion production as [ FraOOa : cr^ " : cr^ ^ ■ 



7 n — >7r p 



1 : 0.5 : 1.25 : 0.8, see also Fig. ml For the production 



of neutral pseudoscalar mesons these estimates give : a 



7*p^r)(?/)A+ 



< a 



and 



0.1 al P. 



In the — > A DVCS process considered above, both unpolarized and polarized 
skewed quark distributions enter, as shown in Eq. ( |192| ). The longitudinal electroproduc- 
tion of pions in the reaction ^*j^N ^ ttA at large allows to select out the polarized 
quark distributions in the A'^ ^ A transition. 

For charged pions, the leading order amplitudes for longitudinal ttA electroproduction 
on the proton are given by [FraOOa] 



{j*p vr+A") 



4 1 

= — ie 

9 Q 



dz 



1 

^ 2 



dx 



ed 



-V3/2 
-1/2 



1 



A • n 



(47ras) 



AT 



rriM 

+ C4(a:,^,t)^(A.Pn0-A/3)W(p) . (249) 

For neutral pions, the leading order amplitude for longitudinal it A electroproduction 
on the proton is given by 



(7*p-> TT^A 







'le 



4 1 
9 Q 

+1 



dz 



dx 



1 eu + ed 
V2 2 
1 

: + ^ — ie 



X {^^ip) \C\{x,tt)n0 + C2{x,^,t)Ap 



1 



{Anas) 
A • n 



+ C3ix,^,t) {/^np-fiAf}) 

rriN 

+ C4{x,tt)^{A-Pni3~Ap) }N{p) . (250) 



The charged pion production amplitudes (Eq. (249)) receive an important contribu- 
tion from the pion pole contained in the function C2{x, t), which can again be evaluated 
analytically. Using the asymptotic distribution amplitude for the pion, and the pion-pole 
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part of hA(t) given by Eq. (p]), the pion-pole contribution to the vr A++ amplitude of 
Eq. ( |249| ) is given by 

■^^ A++ - pole) = -^e Q (q2) ^ llUA ^^(p') A , (251) 



where Ft, (Q ) represents the pion electromagnetic form factor as given by Eq. (237). 
In Eq. ( p51| ), we have introduced the TriVA coupling Jt^na by using the "flavor-skewed" 
Goldberger-Treiman relation 

gA „ \/2 /ttATA 

"T = 2— , (252 

which is equivalent to the ordinary Goldberger-Treiman relation through Jt^na = "i/ ^/^firNN ■ 
Using the phenomenological value /^TVAf/^"" ~ 0.08, this yields fj^NA ~ 2.13. 

The leading order ttA electroproduction cross sections are also shown in Fig. |4^. One 
sees that those ttA channels which contain a charged pion pole exchange give again large 
cross sections, e.g. the 7r~A++ cross sections are comparable with the tt+tt, cross sections. 
Comparing the 7r^A++ and tt+A" channels, their ratio is around 2.25 in our calculations, 
close to the ratio 3 which one would obtain if there was only charged pion pole exchange. 
We found the non-pole contribution to 7r^A++ to be roughly about one-tenth of the total 
7r^A++ cross section. Also the small 7r°A+ cross section, which receives no contribution 
from pion exchange confirms the smallness of the non-pole contributions to the forward 
ttA cross sections in the valence region. 

Similar as for the electroproduction of t:N final states, one can also define a transverse 
spin asymmetry A-kA for ttA electroproduction, which is given by [FraOOa| : 

Xa^-^^"^^^*^^^"^^"^\ with, (253) 

D.,A = \A\^m%{l-£)'' + \B\^e[t' ~2t{m\ + ml) + {m\~mlf] 
+ Re(AB*) 2Cm?^ [C^ - C (3m^ + m%) - t - 171% + m%] . 
For 7r+A° electroproduction , A and B are given by : 

A.+AO = fdxC^! ^f— + . I , (254) 

where the functions Ci, C2 are as defined in Eq. (po|). Eq. ( |255| ) implies that the pion 
pole contribution to B^^+^o is given by ^^+^0 = V3/4-B^+^. We use an asymptotic pion 
distribution amplitude (DA), for which 5^+,'^ is given by : 

^ -3/(20 gA {2mMf/{-t + ml) . (256) 

For the electroproduction of tt^ A++ one has - up to a global isospin factor - analogous 
expressions as Eqs. (254, 25£), by making the replacement <-*■ ea- 

In Fig. ^ (lower plot), we also show the transverse spin asymmetry A.„+/^a for tt+A 
production at several values of — as a function of xb- We use the large Nc relations 
Eqs. ( p3|) for the — > A GPDs. Fig. ^ shows that .4^+ a" has the opposite sign compared 
with A-^ +n- Although the magnitude of ^1,^+^0 is smaller than yl^+„, as anticipated in 
| Fra99a |, it is still sizeable. Similarly, we find At^- /^++ ~ O.5^7r+Ao- 
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5 CONCLUSIONS, PERSPECTIVES AND KEY EX- 
PERIMENTS 



In this work, we outlined in detail the structure of generalized parton distributions (GPDs), 
which are universal non-perturbative objects entering the description of hard exclusive 
electroproduction processes. We focussed mostly on the physics of the internal quark- 
gluon structure of hadrons encoded in the GPDs. For this we presented calculations of 
the GPDs within the context of the chiral quark-soliton model. Guided by this physics 
we constructed a parametrization of the GPDs H , E, H and E, all depending on the 
variable a:, the skewedness f , and the momentum transfer to the target nucleon, t. The 
parameters entering this parametrization can be related in a rather general way to such 
(not yet measured) quantities as the contribution of the nucleon spin carried by the quark 
total angular momentum (J", J'^, etc.), qq components of the nucleon wave function (in 
particular the D-term and the "vector meson" ( VM) part of the GPD E) , the strength of 
the skewedness effects in the GPDs {b^^ai and bgea), the weak electricity GPD in nucleon 
to hyperon transitions (Ewe), flavor SU{3) breaking effects, and others. 

Below we give a list of promising observables which we studied in this work and 
indicate their sensitivity to the different physics aspects of the GPDs. We also discuss for 
each of them the present experimental status and perspectives. We concentrate on deeply 
virtual Compton scattering (DVCS) and hard meson electroproduction (HMP). 



DVCS : single spin asymmetry (SSA) ( see Sec. 1^ ) 



The DVCS SSA, using a polarized lepton beam, accesses (in leading order) directly 
(i.e. without convolution) the GPDs along the line x = e.g. _ff'?(^,^,i). We found 
that the DVCS SSA displays : 

— a sensitivity to the strength of the skewedness effects in GPDs (parametrized 
by byai and bsea) 

— a modest sensitivity to the total angular momentum contributions of the u- 
quark (J") and c?-quark (J**) to the proton spin 

— an important sensitivity to the t-dependence of the GPDs. In particular in 
kinematics where the Bethe-Heitler (BH) contribution dominates, the SSA is 
basically given by a ratio of the imaginary part of the DVCS amplitude to the 
BH amplitude. Therefore it depends linearly on the GPDs and therefore the 
parametrization of the i-dependence directly influences the size of the SSA. 

On the experimental side, the first data for the DVCS SSA have been obtained very 
recently, and several projects both by the HERMES Collaboration and at JLab are 
already planned in the near future. 



First measurements of the SSA were performed by H ERMES j|AhOl| (27 GeV 



beam energy), and at Jlab by the CLAS Collaboration jStepOl, GuidOlj (around 
4.2 GeV beam energy) and display already a clear "sin $" structure of the DVCS 
SSA, indicating that the higher twist effects are modest for this observable at 
the presently accessible scales. 
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Further measurements of the DVCS SSA are planned at JLab (at 6 - 11 GeV 



beam energy), both in HaU-A |BerOO|, and by the CLAS CoUaboration jBurOlb 
BurOlal. 



At HERMES a dedicated large acceptance recoil detector project has been pro- 
posed |Kai01 . This will allo w a su bstantial improvement over existing mea- 



surements of the DVCS SSA | Air01 |, in particular regarding the exclusivity of 
the reaction by detecting the recoiling proton. 



DVCS : cross section and charge asymmetry ( see Sec. 4.3.4 ) 



One accesses the real part of the DVCS amplitude by detecting the DVCS cross 
sections or the DVCS charge asymmetry (i.e. the difference between the reactions 
induced by an electron and a positron). The real part of the DVCS amplitude 
contains (in a convolution) the information on the GPDs away from the line x = ^. 
In particular we have demonstrated the foUowings points. 

^ The DVCS charge asymmetry displays an important sensitivity to the D-term, 
which is required by the polynomiality condition of GPDs. The D-term leads 
to a DVCS charge asymmetry with opposite sign as compared to the charge 
asymmetry when this D-term would be absent. 

— It was discussed how this D-term is a consequence of the spontaneous breaking 
of the chiral symmetry and an estimate of the D-term in the chiral quark-soliton 
model has been given. 

— More generally, the GPDs in the region — ^ < x < ^ contain totally new infor- 
mation on nucleon structure compared to the forward quark distribution. In 
particular, in this region the GPDs behave like a meson distribution amplitude 
and one accesses qq (mesonic) components in the proton state. 

— We have shown that the real part of the DVCS amplitude displays quite some 
sensitivity to the twist-3 effects, which were estimated in the Wandzura-Wilczek 
approximation. 

On the experimental side, different facilities are complementary in the measurement 
of the DVCS charge asymmetry and DVCS cross sections. 

— The measurement of the DVCS charge asymmetry provides a unique opportu- 
nity for HERMES where both e+ and e~ beams are available. A dedicated 
measur ement of the DVCS charge asymmetry is planned in the near future 
|lKaiOl| . 

— The measurement of the DVCS cross section is favorable in a regime where the 
DVCS dominates over the Bethe-Heitlcr process. This is the case at high beam 
energy. An experiment to measure the DVCS cross sectio n at CO MPASS (in 
the 100 - 200 GeV beam energy range) has been proposed [ l'Ho99 . 



At small x b, first dat a for D VCS cross sections have been obtained by the 
HI I Adl01| and ZEUS [3au99| Collaborations, and can map out the behavior 
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of GPDs in the small xb region. These first DVCS data at small xb are in 



agreement with theoretical predictions [Fra98b, Fra99a , implying that higher- 
twist effects are rather small for > 4 GeV'^. 

HMP : Transverse spin asymmetries for longitudinal vector meson production 



( see Sec. 4.5.1 



We showed that in the electroproduction of longitudinally polarized vector mesons 
on a transversely polarized proton target, the GPD E enters linearly in the spin 
asymmetry. Therefore, this transverse spin asymmetry provides a unique observable 
to extract the GPD E and the quark total angular momentum contribution J'' to 
the proton spin. 

— We found large asymmetries for longitudinally polarized vector meson electro- 
production, which display a clear sensitivity to the quark total angular momen- 
tum contribution J" and J'' to the proton spin. Specifically, p*^ electroproduc- 
tion is mainly sensitive to the combination 2 J" -|- J'', p"*" electroproduction to 
the combination J" — J'' and ui electroproduction to the combination 2 J" — J''. 

— To control the model dependence in the extraction of the values J" and J"^, 
the measurement of the transverse spin asymmetries for p^, j and p J electro- 
production over a large kinematic range may be helpful as an additional cross 
check on the modeling of the functions E"^ and E"^. 

— Because the transverse spin asymmetry involves a ratio of cross sections, it is 
expected to be a clean observable already at accessible values of (in the few 
GeV^ range), where the precise form of the t-dependence of the GPDs, the NLO 
corrections and higher twist effects drop out to a large extent. 

These transverse spin asymmetries for longitudinally polarized vector meson electro- 
production have not yet been addressed experimentally until now. At HERMES, an 



experimental program with a transversely polarized target is just starting |Now01|. 
In view of their large sensitivity to J" and J'^, we strongly encourage dedicated 
experiments, at different facilities, to access these transverse spin asymmetries. 



HMP : Transverse spin asymmetries for pseudoscalar meson production 



( see Sees. 4.5.2, 15.3, 4.5.4 ) 



The electroproduction of charged pseudoscalar mesons 7r+ , on a transversely 
polarized proton target also leads to large asymmetries due to the pion (kaon) 
pole contribution to the GPD E. 

Given the knowledge of this pion (kaon) pole contributions, this asymmetry is 
mainly sensitive to the GPD H{£,, t), i.e. along the line x — With the help 
of the flavor SU{3) relations for the N ^ Y GPDs one accesses in this way 
polarized strange quark and antiquark distributions in the nucleon. 
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— Transverse asymmetries for kaon production can teach us about SU{S) sym- 
metry breaking effects. In particular, the K'^ production gives an access to the 
antisymmetric part of the kaon distribution ampUtude which is due to the flavor 
SU{3) breaking effects. Also the transverse spin asymmetry for kaon produc- 
tion is sensitive to the effects of second class currents, especially to the weak 
electricity. 

Besides these observables which have been discussed above, other observables can be 
valuable in the study of the onset of the scaling regime, such as e.g. cross section measure- 
ments for meson electroproduction ( for exp erimental plans at JLab, see Ref. lGuid98| , 



and projects at COMPASS, see Ref. Poc99 ) and double spin asymmetries (DVCS with 
polarized beam and target). 



For a more detailed extraction of the physics content of the GPDs from observables 
further theoretical work is definitely needed. Let us just give an (incomplete) list of 
theoretical problems which, in our opinion, need further clarification. 

• Calculations of NLO perturbative corrections to hard meson production amplitudes. 
For DVCS this work is largely done. Apart from a more accurate access to the GPDs, 
a NLO analysis of the data would allow to perform a global analysis of various data 
on hard reactions from small to large xb- 

• A major open theoretical issue is to get a better understanding of the power (higher 
twist) corrections to meson production which are expected to be sizeable. For DVCS 
observables the immediate task is to quantify the contributions of the target mass 
corrections and the genuine twist-3 operators related to quark-gluon correlations in 
the nucleon. 

• The analysis of the higher twist effects can be useful in identifying observables for 
which higher twist effects tend to cancel, such as expected in transverse spin asym- 
metries for meson production. 

• A further effort is required in modeling the t— dependence of the GPDs. In this 
respect the emerging chiral theory of hard processes can be of a big help. 

• Studying hard exclusive processes such as DVCS in the presence of an additional soft 
pion. 

• Developing the theory and phenomenology of a new resonance spectroscopy with 
light-cone probes can be complementary to traditional low energy probes and bring 
us a breakthrough in the understanding of hadrons as bound states of quark and 
gluons. 

• Establishing links between GPDs and other nucleon structure quantities, such as 
generalized form factors, polarizabilities, etc. 

• Studying the transition from exclusive — s- semi-exclusive — > semi-inclusive — > inclusive 
processes. 
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Summarizing, one may say that the study of hard exclusive processes has generated 
a lot of activity and excitement over the last few years, bringing together several fields 
investigating low and high energy strong interaction phenomena. One may expect new 
discoveries in the near future both on the theoretical and experimental side. We hope that 
the underlying work can stimulate such further efforts. 
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